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ANNALEN DER PHYSIK 


6. FOLGE ~ BAND 11, HEFT 2-3 x 1952 


Die Miesche Theorie der Beugung 
durch dielektrische Kugeln mit absorbierendem Kern 
und ihre Bedeutung für Probleme der interstellaren Materie 
und des atmosphärischen Aerosols 


Von A. Giittler 
(Mit 10 Abbildungen) 


Inhaltsiibersicht 


Im ersten Teil der vorliegenden Untersuchung wird die Miesche Beugungs- 
theorie auf kugelférmige dielektrische Kérner mit absorbierendem Kern ange- 
wandt. Die Amplituden des gebeugten Lichtes werden exakt angegeben. Dann 
werden fiir die optischen Wirkungsquerschnitte sehr kleiner Teilchen Reihen- 
entwicklungen angegeben und schlieBlich aus der physikalischen Bedeutung der 
Rayleighschen Streuung einfache Näherungsformeln fiir die a a 
schnitte gewonnen. 

Der zweite Teil enthält die numerische Durchführung der Theorie für kleine 
Wassertropfen mit Eisenkern im sichtbaren Spektralgebiet als Beispiel. Dabei 
wird die Brauchbarkeit der Näherungsformeln nachgewiesen. Aus den Berech- 
nungen dieses Teils werden einige allgemeine Schlüsse über die optischen Eigen- 
schaften von Trübungen des Aerosols gezogen. 

In einem den Anwendungen gewidmeten Schlußteil wird zunächst kurz ge- 
zeigt, daß bei den interstellaren Körnern eine Schalenstruktur mit absorbierendem 
Kern sehr wahrscheinlich ist. Ein Versuch, die interstellare Verfärbung mittels 
kleiner Teilchen der genannten Art zu interpretieren, führt dagegen wegen deren 
zu kleiner Albedo zu keinem befriedigenden Ergebnis. Anschließend wird in einem 
Schlußabschnitt über die opaleszente Trübung der Atmosphäre der optische Unter- 
schied zwischen absorbierenden und dielektrischen Kondensationskernen unter- 
sucht. Während diese i. a. die optischen Eigenschaften des Aerosols kaum ändern, 
erhöhen absorbierende Kerne‘ die Streuung und ganz besonders die Extinktion 
der Wassertröpfchen bedeutend. 


I. Die Miesche Theorie der Beugung durch dielektrische Kugeln 
mit absorbierendem Kern 


1. Berechnung der Wirkungsquersehnitte 


Als Modell betrachten wir eine Kugel vom Radius R, aus absorbierendem 
Material (Brechungsindex m, = n- (1—i x)) mit einer umgebenden dielektrischen 
Kugelschale mit dem Außenradius R (reeller Brechungsindex m,) (s. Abb. 1). 
Die Berechnung der Wirkungsquerschnitte eines solchen Korns für Streuung, 
wahre Absorption und die gesamte Extinktion läßt sich nun in ähnlicher Weise 


Ann. Physik. 6. Folge, Bd. 11 


65 
7 
113 | 
131 # | 
145 
155 
161 
166 = 
169 
175 
rher 
des 
ing, 
der- 
üge 
den- 
‚12, 
rek- 
ind 
Ver- 
len. y 
tel- 
igl. 
5 


66 Annalen der Physik. 6. Folge. Band 11, 1952 


wie bei der Mieschen Theorie der Beugung an homogenen Kugeln!) durchführen, 
die u.a. bei Born?) ausführlich dargestellt wird. Dort findet man weitere Literatur- 
hinweise für diese Theorie, und auch die hier benutzten Bezeichnungen schließen 
sich zunächst weitgehend dieser Darstellung an. Einen sehr instruktiven Über- 
blick über alle bisher vorliegenden derartigen Untersuchungen gibt überdies 

van de Hulst®) (vgl. insbesondere 1. c. 8.13). Über das Wesentliche des Mie- 
: schen Verfahrens diirfen hier 
wohl einige sehr kurze Bemer- 
kungen geniigen. 


falle eine ebene Welle (e) auf 
das Korn. Die dadurch erregten 
Feldstärkewellen in dessen 
Kern (I) und Hülle (II) sowie 
die gebeugte Welle (b) im 
Außenraum (I/II) ergeben sich 
Abb. 1. Modell eines dielektrischen Teilchens mit nme bekanntlich en 
Harn Maxwellschen Gleichungen zu- 
sammen mit den passenden 
Randbedingungen. Der durch (e) geforderte periodische Ansatz (E, 9 m ei®t) 
führt nun auf zwei Wellenpotentiale J7’ und //’, die elektrischen bzw. magne- 
tischen Multipolschwingungen zugeordnet werden können und der Differential- 
gleichung 


AI + IT =0 (1) 
genügen. Aus diesen lassen sich die Feldstärke-Amplituden berechnen gemäß 
(r IT) ow? 
Oreo ce r-sind 
__1 io, 1 
62 (r IT’) 
_ 1 If) 1 ‚@(r II’) 
i... 1 (r IT’) 1 (r IT”) 


Die Lösungen müssen eindeutige Ortsfunktionen sein und außerdem die fol- 
genden Randbedingungen erfüllen: 

1. Beschränkte Amplituden auch im. Zentrum des Korns. 

2. Lösung im Außenraum I/II additiv zusammengesetzt aus (e) und (b), webei 
(e) als einfallende ebene Welle vorgegeben und (b) im Unendlichen eine auslaufende 
Kugelwelle (Ausstrahlungsbedingung): 

3. Stetiger Übergang der Tangentialkomponenten’ der Feldstärken an den 
Grenzflächen I/II (r = R,) und II/III (r = R). 


1) G. Mie, Ann. Physik (4) 25, 377 (1908). 
2) M. Born, Optik, Berlin 1933, S. 2 
3) H. C. van de Hulst, Rech. ‘aan. Obs. Utrecht 11/1 (1946). 


Aus dem Außenraum (III): 
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Mit dem Produktansatz 
II (r, p, 8) = F (r) -® (9) (0) (3) 
gelingt nun in bekannter Weise die Separation der Variablen, und wenn man noch 
einführt, so erhält man als Eigenfunktionen 
k 
— (0) Pl (cos 8). ©) 


Dabei ist = 0,1, 2,3,...und ~=—l,...,+ 1, und die y, und x, bedeuten 
die Bicchti-Resselechan bzw: -Neumannschen Funktionen 


¥s +0 für z 0 (6) 


(2) = > 00 


Die erste Randbedingung läßt im Kern (J) nur die y, zu, während in den beiden 
anderen Räumen JJ und III passende Linearkombinationen der y, und x, zu be- 
stimmen sind. Für die gebeugte Welle (b) sind diese bis auf einen noch zu berech- 
nenden konstanten Amplitudenfaktor die Funktionen 


entsprechend der zweiten Randbedingung. 
Entwickelt man die einfallende, auf die Intensitätseinheit normierte ebene 
Welle (e) nach den Eigenfunktionen, so erhält man 


r- IT, g- Pi (cos #) sin p, 
womit im Außenraum 
IT, = IT, + H,. (9) 


Wegen (8) lassen die letzten beiden Randbedingungen in allen drei Räumen nur 
= 1 zu, so daß ” die Lösung in der Form 


SA (2) 
r- Ih = 4 -y, (kr) - P, (cos®)cosp 


r- Ty = = Ai’ - ; (ky 7) Pi (cos 8) 


r- IT, = (Bi vi (kar) + (ker) - Pi) (c08 8) cos 
(10) 


or (By (har) + (kar) PP? (c08 8) sin 


= (kat) Pi? (cos 8) cosp 


r- Di’ +f, (kar) - P (cos 8) sing 
schreiben läßt. 
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Aus der dritten Gruppe der Randbeaingungen ergeben sich die Gleichungen 


ki R, IT; (R,) = ks R, -IE (R,) 
Ry IR’ (Ry) = Ry (Ro) 


ky R-II, (R) = ky- R-II; (R) 


R-IE (R) — 


ö 73 | | 
ar Ih) II) 
ö ” 
1) = (r-H; Re 
ut é 
R-II (R) 


wobei die beiderseits gleichen Argumente y und # weggelassen wurden. Aus (8), 
(9), (10) und (11) ergeben sich nun für die Koeffizienten A, B, C, D die folgenden 
Bedingungsgleichungen: 


yı(kı R)- A; 
Be) gy villa 
k, ky 


(ky R) xi (ky R) 
ky ks 


Yyı (k, Ro) (ky Ro) 


(ke Ro) - 
(ke R) - By+ (ke R) (ks R) - D, 


By — (ke Ro) €; 
Xi (ke Ro) 
ky 


21+1 
y, (ks R) 


.D; — ql- ‚2i+1 wi (ks R) 


ok 
(ky Ro) - AY— yi (ke Ro) By — (ke Ro) Ci’ 0 


Aus diesen acht linear unabhängigen Gleichungen können die Amplituden der 
gebeugten Welle sofort bestimmt werden. Es ergibt sich 


mit 
i (k, Ro) —k,- (k, Ry) hi (k, Ro) 0 | 
| (kı Ro) (ka Ro) — (ke Ro) 0 
| 0 ky (k, R) ks (ks R) (13) 
Vi (ke R) 
kıyılkı Ro) — (he Ro) — ke - (he Ro) 0 
(kı Ro) — (k, Ro) — x1 (ke Ro) 0 
| 0 “ky (ky R) ky (ks R) 
0 y; (k, R) Xi (k, R) Ci (ks R) | 


| 
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Yı (kı Ro) —yı (k, Ro) — Xı (ke Ro) 0 
- (k, Ro) —k,- (k, Ro) —k,- X (k, Ro) 0 
0 yı(k, R) Xı (ky R) yı(k, R) 
0 ky 177 (k, R) ky - (k,R) ky-yi(ks R)| (13) 
vi (ki Be) Be) — (ka Bo) 0 
ky - vi (k, Ro) (k, Ro) —k,- Ui (k, Ro) 0 
0 yı (k, R) Xı (ke R) &, (ks R) 
0 k, - (k, R) k, (k,R) k,- Ci (ks R) |, 


womit wir bereits zu den von van de Hulst®) gebrauchten Bezeichnungen über- 
gegangen sind, wenn wir noch 
2nR_ 


7 


(14) 
setzen. Der von Greenstein eingeführte Extinktionsfaktor, d. h. das Verhältnis 
des Extinktionsquerschnitts zum geometrischen Querschnitt des Teilchens, ist 
dann 


ah +4, + (214 1)- Re [a, + (15) 


Und auch die Gleichungen für die entsprechend definierten Wirkungsfaktoren 
für die gesamte Streuung, die Vorwärtsstreuung, die wahre Absorption und den 
Strahlungsdruck nehmen in der von van de Hulst*) angegebenen Schreibweise 
eine sehr einfache Form an. Der Streufaktor wird nämlich 


= 1) - (var + bf) 1) (ja,|? + (16) 
und der 
(ay + af + by + OF by 4.1) 


4 + 2 2l 


Und hieraus erhält man den Absorptionsfaktor A und den Wirkungsfaktor P/c 
für den Strahlungsdruck gemäß 
A=E-—-S 


P=E—F. 


(18) 
(19) 


Die auf die Volumeneinheit des von den Teilchen erfüllten Raumes bezogenen 
Werte der Extinktion usw. erhält man aus diesen Wirkungsfaktoren durch Multipli- 
kation mit N -2 R?, wenn N die Teilchendichte bedeutet. 
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Zur Auflésung der Determinanten in (13) dienen die folgenden Abkiirzungen: 


= xi (ke R) - ky yp, (ky R)—k, Xı (ke R) yi (ky R) 
[x Hi (ke R) - ky (ky R)— ky (he R) - (ky R) 
yı (k, R)- kyy, (ky R)— ky p, (k, R) - yj (ks R) 
= yi (ke R) - ky (hs R)— kz, (kg R) - (ks R) 
Ixvlin = ke xi (ka R) - y, (ks R)— (hy R) - key yi (hy R) 
[x =k, (k, R) - &,(k, R)— (ke R) ky Ci (ky R) (20) 
[yylir = kayı (ky R) - (ks R)—y; (ky R) ey yt (ky R) 
= ke Yi (ke R) - (ks R) —y, (ke R) - ky Ci (k, R) 
= xi (ke Ro) kıy, (kı Ro) — k, Xı (he Ro) - yi (kı Ro) 
[yyY, 7 (ke Ro) - ky (ki Ro) — ke Yı (ke Ro) i (kı Ro) 
[x R, = he x1 (he Ro) - (hi Ro) — Xı (he R,) kı yi (k, Ro) 
[py]: Ry = Ie yi (he Ro) (hi Ro) (he Ro) kı yi (kı Ro). 
Man erhält dann nach kurzer elementarer Rechnung 


Man überzeugt sich leicht, daß diese Ausdrücke für R, — 0 (homogene dielek- 
trische Kugeln) in die entsprechenden Determinantenquotienten der Mieschen 
Theorie der Beugung an homogenen Kugeln übergehen. Dasselbe läßt sich auch 
für R,— R (homogene absorbierende Kugeln) leicht zeigen, wobei man zweck- 
mäßig von der Determinantendarstellung (13) ausgeht. 

Das allgemeine Verhalten der Amplituden für beliebige Teilchengrößen zu 
diskutieren, ist eine Aufgabe, die durch die in den Amplituden auftretenden Bessel- 
Funktionen beliebiger komplexer Argumente des vierten Quadranten sehr er- 
schwert wird, und auf deren Behandlung wir hier wegen des dazu erforderlichen 
Aufwandes an Rechnung leider verzichten müssen. Man würde sich bei einer 
derartigen Untersuchung vor allem auch mit der Erscheinung der optischen Re- 
sonanz befassen müssen, die jedesmal dann eintritt, wenn für einen speziellen Wert 
von x in einer der Amplituden (13) die Nennerdeterminante sehr klein wird ohne 
ein entsprechendes Verhalten des Zählers, so daß diese Amplitude einen extrem 
hohen Wert annimmt. Es wäre dann nachzuweisen, daß dabei ein völliges Ver- 
schwinden des Nenners in keinem Falle möglich ist, da man sonst einen negativen 
Absorptionsfaktor erhalten würde. 

Aber wenn man auch diese Fragen übergeht bzw. ihre Beantwortung in der 
angegebenen Weise vorwegnimmt, so bleibt es noch wünschenswert, etwas über 
die Konvergenz d r Reihen (15) bis (19) für die Wirkungsfaktoren zu erfahren. In 
Anbetracht der mit einer vollständigen Durchführung der Konvergenzbeweise 
verbundenen recht umfangreichen Rechnung wollen wir uns aber auch hier mit 
dem Hinweis begnügen, daß die Amplituden a, und b, für » = 21-+ 1 — oo ganz 
wie im Fall homogener Teilchen von den Resonanzschwankungen abgesehen wie 
(z/v) verschwinden. Man übersieht dies schnell, wenn man in die in (13) vor- 
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en: kommenden Determinanten die von Debye‘) angegebenen asymptotischen Ent- 
wicklungen der Zylinderfunktionen fiir groBe Indizes und beliebige positiv-reelle 
bzw. komplexe Argunente des vierten Quadranten einsetzt. Aus diesem starken 
Verschwinden folgt dann aber sofort auch die Konvergenz der Reihen fiir alle 
Teilchengrößen. 


2. Entwicklung der Wirkungsfaktoren für kleine Teilchen 

Nach Einführung des Radienverhältnisses 

20) om > (22) 
lassen sich die Amplituden a, und b, nach Potenzen von R entwickeln, wofür sich 
am besten wieder die Determinantendarstellung (13) eignet. Man erkennt sofort, 
daß die Entwicklungen der vier in diesen Gleichungen vorkommenden Deter- 
minanten jedesmal mit einem Gliede (k, R)?!*" beginnen. Um die Entwicklungen 
der Wirkungsfaktoren (15) bis (19) bis zu den Gliedern vierter Ordnung zu erhalten, 
reichen daher die den Reihen aus: 


21) (23) 

Ya(2) =240@) Yo (z) = 

3 
ale) +0@) ne) +0). 
len Mittels der ‘ites Bezeichnungen gelingt nun yi gleichzeitige Behandlung 
ich _ der elektrischen und magnetischen Amplituden 
ok- und aj’ =), (24) 
al | = kılks = m 1=1 
92 = ka/ks = m; he (25) 
er g=h =m 
en 1 1 1/%3 1 
1er 9 = hy = = m,. 
te- Die folgenden Gleichungen gelten nunmehr in gleicher Form für die elektrischen 
art und magnetischen Schwingungen, so daB die diese unterscheidenden Striche an den 
ne oben eingeführten Bezeichnungen vorerst wegbleiben können. 
m Ist nun noch 9, = y, bzw. x,, so lassen sich alle in (13) vorkommenden Deter- 
minanten in folgender Form schreiben: 
en 
R) ys (leg R) 0 
1 

In 0 (k, R) (k, R) (ks R) em 
| | , 
nz = fi (R) - 
ie 4) P. Debye, Jahresber. Dt. Math.-Ver. 1933; — s. a. E. Jahnke u. F. Emde, 


Funktionentafeln, 2. Aufl., Leipzig-u. Berlin 1933, 8. 204. 


. 
T- 
. 


= 19- 
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wenn ein kürzbarer Faktor fi (R) herausgezogen wird. Setzt man hierin nun die 
Entwicklungen (6) ein, so erhält man für 1 = 1 mit 


1k 
h= R 


( 


-( 
2h, -(1— (ky (i— Ay 0 
( 


q 
0 2h, -(1— 2") 17 


Bei der Entwicklung dieser Determinante brauchen alle Glieder von höherer als 
dritter Ordnung nicht mehr berücksichtigt zu werden. Mit (14) erhält man daher 
in dieser Weise 


Die Koeffizienten & und 7 bedeuten dabei die folgenden Polynome in q: 


& = 20h, {Gihs + 29h) — q} 
& = 2h, -{— 2m} - (gh, + 
— [my + + 5 my 
+ [mi - (91 4g — 2 + m - (29: 92 h;)] - q°} (28a) 
{(gı hs + 29,h,) + 2 (gh, — - q*} 
{m; - (91 + 29, hy) 
+ [my (gr he + 4.921) + 5 my - 2 92h,)] - 
— 10m; - (9: h2— - 
+ 2 [mi (gi + me - (2 —G2h,)] - q°}. 


Für / - 2 erhält man in ähnlicher Weise mit 


h= 75 ke 
0 
=| 1 +. (29) 
Is 92:9? ER 
| 0 3h, — 2h, q? 


und hieraus weiter wie oben 


mit (30) 
E, = 3h, {(2gih, + 2 - | 
=— {(2gihe +39 h,) +3 (gi he— 9: hı) 


(21) 


mit 


mit 


wo 


un 


fa 


I 

weg 

Dab 

und 

= 

Da 

sch 

= 

so 

jer 

x 


(21) 


4. Güttler: Die Miesche Theorie der Beugung durch dielektrische Kugeln 73 


Hieraus berechnet man nun die elektrischen und magpetischen Amplituden 
wegen (13), (24) und (26): 


[ye] _ _ a, [yi] (31) 


Dabei ist nach (28) und (23) 


= {A + at + 

mit 

(32) 


und 

[xa] = (1 + &ı 2°) + =) — (m+ x*) 

mit 


X= et m) 


Damit erhält man nach (31) für die Amplituden der elektrischen und der magneti- 
schen Dipolschwingung 


wobei jeweils (33) 


In derselben Weise erhält man aus (30) und (23) zunächst 


De [ya] ==. ¥, mit ¥, 
und (34) 
3 
Ws [x2] X, mit X, = 2 ne, 


so daB nach (31) fiir die elektrische und die magnetische Quadrupolamplitude 
jeweils 


mit (35) 


Bevor diese Amplituden in die Ausdrücke (15) bis (19) für die Wirkungs- 
faktoren eingesetzt werden, müssen die Polynome &(q), n(g) und danach ¥(q), 
X(q) für die elektrischen und magnetischen Oszillationen getrennt berechnet 


S 

r 

x Pı Pr xX, » | 
| | | 
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werden. Mit (25) findet man aus (28a) und (30): 
& = (20 m} + 40 mi) — (20 m! — 20 mi) - q° 
&=— (4 mim} + 8 mi) 
— (2m; + 18 my m; — 20 m;) - g? 
+ (10 m; m; — 10 m;) - q* 
+ — 2 m3) - 


71 = — (10 my m; — 20 m3) + (2 mı — 2 mi) - g? 
= (mi m; + 2 ms) 

+ (mi mz + 9 m} m; — 10 ms) - q? 

— (10 m} m; — 10 ms) - q? 

+ (2 mj m; — 2 ms) - 95 
& = (6 mi + 9 m3) — (6 mi — 6 mi) - g5 (36) 
72 = — (2 mi m; + 3 m3) — (3 mim; — 3 m}) - 45 
60 m, m; 
& 12m, m} 

— (10.m} m; — 10 m, m}) - q? 

— (2 m} m; — 2 m, ms) 
= — 30 m, m; 
n= 3m, m; 


+ (5 mı m; — 5 m, m3) - q? 
— (2 m} m; — 2 m, m}) - 95 
& = 15m, m; 
"a 5 m, m}. 
Und hieraus wieder findet man nach (32) und (34) die Polynome 


Wo + Ws Xi = + 


Pi = Wy + Wiz + + Wis i (37) 
Xi = Vo + + + 


= Wig + Xz = + 


wobei — jetzt wieder getrennt fiir die elektrischen und magnetischen Amplituden — 


w), = — 20 (m? — 1) - (m? + 2 m?) 
= — 20 (m? m? + 2 mi — m? — 2 m?) 
Wi, =— 20 (2m} + 1) - (m?— mi) 
= — 20 (2 m? m? — 2 m! +- m? — m?) 
vi, = 10 (m2 +'2) - (m? + 2 m?) (38) 
= 10 (m? m3 + 2m$ + 2m? + 4m!) 
= 20 (m}— 1) - (m} — m}) 
= 20 (mm? — m$—m? + m2) 
Wyo = 2 (m}— 1) - (m? + 2m?) 
= 2 (m?m$ + 2m — m? — 2 m2) 


Ww 
. 
4 
u 
[2 


(36) 


37) 


— 


32) 
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Wi, = 2 (m3 — 1): (mt + 9 m? mi — 10 md) 

= 2 mt m? + 9 m? m! — 10 m$ — mt — 9 m? m? + 10 mt) 
w), = — 2 (10 mi — 9 m? — 1) - (m? — m?) 

= — 2 (10 m? m$ — 10 m8 — 9 m? m? + 9 mi — m? + m?) 
©, = 2 (2m? + 1) (mt— mi) 

= 2 (2 mi m? — 2 m + mi — ms) 
Bf) =— (mt + 9m2— 10) - (m? + 2m}) 

= — m? m} — 2 m$ — 9 m? m? — 18 m! + 10 m? -+ 20 m? 
= — + 2) - (mi + 9 mı m? — 10 m;) 

= — mi m; — 9 mi m; + 10m; — 2 mt — 18m! mi + 20m: 
Bis = 10 (m; — 1) - (m — m;) 

= 10 (mi m; — m; — mı + m,) 
— 2 (m; — 1) - (m; — m;) 
— 2 (mi m; — m; — m; + ms) 
Way = — (m — 1) - (2m; + 3m,) 

— 2m, m; — 3m; + 2m; + 3m; 
Whs = — 3 (3m; + 2) - (m — my) 

= — 3mm; + 3m — 2m; + 2m, 
un = (2m; + 3) - (2mı + 3m;) 

= 4 + 6 ms + 6 mi + 9m 
6 (m5 — 1) - (mi— mi) 

= 6 (mim; — ms — m} + mi) 


I 


Vis 


wio = 0 wi; — 0 
Vio = W m, m; V3 = 0 
Wyo = — 6 m,m; - (m;— 1) = — 6 (m, my — m, m3) 
Wi = u; = 0 
Wis = — 6 m,m, - (m) —m,) = — 6 (m? m, — m, m2) 
55 = — 15 my mi — 1) = — 15 (m, mi — m, mi) 
= — 15 m, m; (m} — m3) — — 15 — m, m3) 
Wyo = ws = 0 
= 5m, = 0. 
Da nach diesen Gleichungen 
1 = 90m, m; X, = 5m, m, 
wird 
p; -0 p, = 0 


1 2 2 
+ mini) <q] 


und fiir reelles m, daher 
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(38) 


} 
| 
= 
| (39) 
u 
| (40) 
Im py” = Immi. (41) 
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Damit ergibt sich aus (33) und (35) fiir die elektrischen und magnetischen Dipol- 
und Quadrupolamplituden mit den Pg ei nach van de Hulst 


a, = —t-p,- x +.” 

6, = 


Wie 6, enthalten auch alle Amplituden der höheren Multipolschwingungen keine 
Glieder bis zur 6. Ordnung mehr. Vergleicht man die Gleichungen (42) mit den 


entsprechenden Entwicklungen für homogene absorbierende Teilchen®), die von. 


Schoenberg und Jung)®) angegeben wurden, so erkennt man, daß die bei den 
einzelnen Amplituden vorkommenden Glieder, soweit sie hier berechnet wurden, 
genau die gleichen geblieben sind. Lediglich die Werte der nicht verschwindenden 
Koeffizienten haben sich geändert. Aber auch die Beziehung zwischen dem ersten 
und dem dritten nicht verschwindenden Koeffizienten in a, ist, wie (42) zeigt, 
erhalten geblieben. Daher sehen auch die sich hieraus ergebenden Entwicklungen 
der Wirkungsfaktoren ebenso aus wie bei homogenen absorbierenden Teilchen. 
So erhält man nach (15) für den Extinktionsfaktor wieder 

mit 

= 6Im py (43) 

= — 6 Im pi — 6 Im + 10 Im 

e = 6 Re pi’ = 6 (Re 91)? — 6 (Im pi)? 
und nach (16) fiir den Streufaktor 

mit 

8, ~ 6 91? = 6 (Re pi)? + 6 (Im 21)? . (44) 

8, — 12 Re (pi = —28,- Red 

8, =— 12 9° 

worin 6 die gleiche Bedeutung hat wie in (33). Ähnlich wird nach (17) der Wir- 

 kungsfaktor für die Vorwärtsstreuung 


mit 
fa = 2Re (Pı 4*) 
fa = — 2 Re (p; -4*) (45) 
Ja = — 2 Im (pi - 4*), 

wobei 


A=3p,—2 py’. 
Fiir den Absorptionsfaktor ergibt sich dann nach (18) 


| 

mit 
a, = & Ag = | (46) 
a, = &— 8 = — 12 (Im 9)? 


5) E. Schoenberg u. B. Jung, Astronom. Nachr. 258, 261 (1934) = Veröff. Sternw. 
Breslau 7, Nr. 6. 
6) E. Schoenberg u. B. Jung, Mitt. Univ.-Sternw. Breslau 4, 61 (1937). 
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und für den Wirkungsfaktor P/h des Strahlungsdrucks gilt nach (19) 


mit (47) 
Da bei der Herleitung von (33) die Entwicklung von (1 + z)-! benutzt wurde, 
ist der Konvergenzbereich aller dieser Reihen beschränkt. Wie das im folgenden 
Teil ausgeführte Beispiel zeigt, beginnt die Konvergenz daher schon schlecht zu 
werden, sobald die Absolutbeträge der Argumente x, m, und m, der in 
(13) entwickelten Riccatischen Zylinderfunktionen nicht mehr hinreichend weit 
unter 1 liegen. 


3. Sphärische Albedo und Asymmetriefaktor 
Das gleiche gilt für die Entwicklungen der Ausdrücke für die sphärische Albedo 
und den Asymmetriefaktor, die zur Vervollständigung hier noch mitgeteilt werden 
sollen. Unter Benutzung der aus dem vorhergehenden Abschnitt ersichtlichen 
Entwicklungen der Wirkungsfaktoren für Streuung, Vorwärtsstreuung und Ex- 
tinktion findet man sofort für die sphärische Albedo 


8 
mit 
“= 3 e 
und für den Asymmetriefaktor 
F 
mit 
_ fe _ fe—% _ fs (49) 
_h,a 


4. Einfache Näherungsformeln für die Wirkungsfaktoren 
kleiner zusammengesetzter Teilchen 

Wenn|m|-~<1 und außerdem die Brechungsindizes. sich nur wenig von 
1 unterscheiden, ist das Volumen des Teilchens so klein, daß hierin die Felder 
der einfallenden Welle nahezu homogen sind und daher alle streuenden Oscillatoren 
des Teilchens in Phase schwingen. In diesem Bereich der Rayleighschen Nä- 
herung kann man sich alsdann alle diese Oscillatoren im Zentrum des Teilchens 
zusammengezogen denken und kommt dann zu den folgenden einfachen Näherungs- 
formeln: 

Es seien Z,»x und x! die Extinktionsfaktoren homogener absor- 
bierender bzw. nicht absorbierender Teilchen. Dann gilt unter den genannten 
Voraussetzungen für den Extinktionsfaktor eines zusammengesetzten Teilchens 
offenbar 


E(R, a) — B, a) + B((1— 9?) R,A), 


_ 
ne 
en 
On: 
en 
en 
en 
st, 
en 
n. 
r- 
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woraus hier sofort 


E(R, A) = - E,(R, A) + (1— 9°)? - A) (50) 
und ebenso mit $, ~ x* und 8, ~ x* als den entsprechenden Streufaktoren 
+ (1— 8, (RB, A). (51) 


Auch die sphärische Albedo y = S/E läßt sich auf diesem Wege sehr einfach be- 
rechnen. Ferner erhält man ebenso noch die folgenden Näherungsformeln für die 
Wirkungsfaktoren der Absorption x und der Vorwärtsstreuung F ~ 
nach (45) und (46): 


A=q-A, (52) 
(53) 


während fiir den Wirkungsfaktor des Strahlungsdrucks nach (47) einfach P=E 
zu setzen ist. 

In allen Fällen ist die von einem sehr kleinen Kern hervorgerufene Anderung 
der Wirkungsfaktoren proportional zu g*, während andererseits die Änderung der 
Wirkungsfaktoren absorbierender oder auch nicht absorbierender Kerne mit 
einem dünnen rein dielektrischen Überzug hiernach proportional zu 1— gq, d.h. 
der relativen Dicke desselben ist. 

Die Brauchbarkeit dieser Formeln wird im folgenden Teil der Untersuchung 
geprüft werden, wobei es sich zeigen wird, daß sie in den meisten Fällen, wo es 
sich um kleine Teilchen handelt, einen völlig ausreichenden Ersatz für die recht 
mühsame Auswertung der obigen Reihenentwicklungen darstellen und auch die 
Voraussetzurg über die Brechurgsindizes keineswegs so einschneidend sind, wie 
es zunächst scheint. 

Schließlich sei hier noch eine Folgerung aus der vorstehenden Betrachtung 
hervorgehoben, die für viele Anwendungen der Theorie zu beachten ist: Im Gültig- 
keitsbereich der vorstehenden Näherungsformeln besteht kein optischer Unter- 
schied zwischen einer Mischung homogener Teilchen aus verschiedener 
Materie und gleichartigen, aus den gleichen Komponenten schalenförmig 
aufgebauten Teilchen derselben Größe, wenn die Mengenverhältnisse der Kom- 
ponenten in beiden Fällen gleich sind. 


II. Wassertrépfchen mit Eisenkern 
1. Dia Brechungsindizes 


Die Auswertung der exakten Formeln (21) wird leider durch das Fehlen von 
Tafeln der Zylinderfunktiopen für beliebige komplexe Argumente so erschwert, 
daß der Verfasser bei der Erprobung der Theorie an einem Beispiel darauf ver- 
zichten und sich auf die Benutzung der Entwicklungen für kleine Teilchen be- 
schränken mußte. Für den absorbierenden Kern wurden in diesem Beispiel die 
von W. Meier’) im sichtbaren Spektralgebiet gemessenen Brechungsindizes des 
Eisens (s. Abb. 2) angenommen, fiir die dielektrische Hiille die des Wassers bei 


7) W. Moser; Untersuchungen iiber Dispersion und Absorption bei Metallen fiir das: 
sichtbare und ultraviolette Spektrum, Diss. Göttingen 1909, Leipzig 1909, s. a. Ann. Physik 
(4) 81, 1017 (1910). 
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25° C entsprechend den Messungen von Tilton und Taylor’) (s. Abb, 2). 
Diese spezielle Wahl wird besonders angesichts der in diesem Zusammenhang 
bereits in einer früheren Untersuchung des Verfassers!) beschriebenen Unzu- 
verlissigkeit der genannten Brechungsindizes des Eisens allein durch die sich dabei 
ergebende Möglichkeit eines Vergleichs der Ergebnisse mit bereits vorliegenden 
Berechnungen für q—0 (Wassertrépfchen) und g=1 (kleine Eisenkugeln) 
gerechtfertigt, eine Möglichkeit, die sich als sehr nützlich erweisen wird. Aus 
diesem Grunde wurde auch auf eine strenge Übereinstimmung der benutzten 


18 r n 
14 
12 

A 
10 
sa ae. 400 500 600 700 mu 

A> 


Abb. 2. Spektraler Verlauf der Brechungsindizes: m, = n- (1 — ix) für galvanisch zer- 
stäubtes Eisen nach W. Meier’) mit dem hier benutzten geglätteten Verlauf nach (54) 
(— — —), mg: Brechungsindex des Wassers bei 25° C-nach Tilton und Taylor®) 


Brechungsindizes mit den gemessenen Werten verzichtet und ihr spektraler Ver- 
lauf zugunsten einer erleichterten Erkennbarkeit von Rechenfehlern etwas ge- 
glättet. So konnte für den Eisenkern x konstant gesetzt und der spektrale Verlauf 
von n(A) sehr gut durch eine Gerade approximiert werden (s. Abbb. 2),welche die 
Werte für den von Meier angegebenen Verlauf bei 400, 508 und 650 mu exakt 
wiedergibt. 

Auf diese Weise ergaben sich für den Brechungsindex des Eisens zwei sehr ein- 
fache Näherungsformeln, welche mit einer geringen Extrapolation für den ganzen 
Spektralbereich zwischen 400 und 700 mu brauchbar sind, der hier behandelt 
werden soll: 


n=a:(A+ 4) 
mit = 1,91 u! Ay = 0,215 u (54) 
& x= 1,10. 


Die entsprechenden Ausgangswerte für die nachfolgende Berechnung, in die nach 
(38) nur die Quadrate der Brechungsindizes in Produkten und einfachen Potenzen 
eingehen, wurden demgemäß aus 

\ m? = — (0,765 + 8,02 2) - (Au) + 0,215)? 
berechnet (s. Tabelle 1). F 


8) L, W. Tilton u. J. K. Taylor, J. Res. Bur. Stand. 20, 420 (1938); s. a. J. d’Ans 
u. E. Lax, Taschenbuch fiir Chemiker u. Physiker, Berlin 1943, S. 1189. 
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Die Brechungsindizes m, des Wassers wurden der gleichfalls in Abb. 2 einge- 


tragenen Kurve entnommen. Die hier beriicksichtigte Dispersion ist allerdings ~ 


so klein, daß man sie sehr wohl auch ganz vernachlässigen kann, wie es bei ähn- 
lichen Berechnungen häufig geschieht. 

In der nachfolgenden Tabelle 1, welche die für die Berechnung der Koeffi- 
zienten v und w nach (38) erforderlichen m-Produkte enthält, sind die unabhängigen 
Ausgangswerte der m; und m; durch Kursivschrift hervorgehoben. 


Tabelle 1 
Für (38) erforderliche Produkte der’ Brechungsindizes 

A [mu] 400 450 | 500 550 | 600 650 700 

z — 0,28934 — 0,33830 — 0,39109 — 0,44770 — 0,50813 |\— 0,57239 — 0,6408 
— 3.0334 — 3,5466i — 4,1000 |— 4,6935i — 5,3271i — 6,0071i 6.7145: 

22° |— 0,52157 |— 0,60627 0,69805 — 0,79645 — 0,90143 |— 1,01371 |— 1,133% 
mm |_ 5,4679: |— 6,3559: — 7,3181i |— 8,34975 — 9,4502: |—10,6385: |—11,8807i 

2.4. 0,94018 |— 1,08650 |— 1,24595 |— 1,41689 |— 1,59913 |— 1,79527 |— 2.00519 
mim; |_ 985655 |\—11,3905i —13,06215 |—14,8542i —16,7647: —18,8408: |—21,021% 

— 9,1176 |—12,4642 |—16,6573 |—21,8285 |—28,1197 |--35,6816 |—44,6749 
+ 1,7554: |+ 2,39975 |+ 3,2069: |+ 4,20255 + 5,4137: |+ 6,8696: |+ 8,6010i 

{16,4353 |—22,4371 —29,7316 |—88,8329 |—49,8843 |—63,1921 |—79,0478 
|4 3,1642¢ |+ 4,2005¢ + 5,72415 |+ 7,4763: 9,6040% |-+12.16605 |+15,2187 
m; 1,8026 1,7921 | 1,7849 | 1,7790 | 1.7740 1,7710 1,76% 
m; 3,2494 3.2116 3.1859 | 3.1648 3,1471 3,1364 3,1308 
mi | 5,8573 5, 7555 5,6865 5,6303 5,5829 5,5546 5,53% 


Die Werte diese Tabelle wurden ebenso wie die folgenden Zahlenwerte mit der 
Maschine berechnet, um die erforderliche innere Genauigkeit zu erreichen. 


2. Die Koeffizienten des Extinktionsfaktors 


Aus den m-Produkten der Tabelle 1 wurden nun nach (36) zunächst die Koeffi- 
zienten v und w und dann damit die Yund X für die Radienverhältnisse g —0 
(homogenes Wassertrépfchen), 0,2, 0,4, 0,6, 0,8 und 1 (homogenes Eisenkügelchen) 
berechnet. Das Ergebnis zeigt die folgende Tabelle 2. 


Tabelle 2 
Die Koeffizienten Y’ und X’ nach (37) 
A[mp]| 400 450 500 550 | 600 650 700 
q 
00 53226. |— 51,421 |— 49,900 |— 48,459 |— 47,057 |— 45,791 |— 44,600 
|-+ 48,6911 |+ 56,186: |+ 64,3624 |+ 73,12bi |+ 82,4631 + 92,6294 |+103,323i 
0.2 |— 51,685 |— 49,859 |— 48,309 |— 46,835 |— 45,397 |— 44,088 |— 42,850 
50,9274 + 58,7874 |+ 67,360: 76,548: |+ 86,3401 |+ 96,994: |+108,200i 
04 40,895 |— 38,921 |— 37,172 |— 35,468 |— 33,772 |— 32,167 30,5 
66,5725 |+ 76,997: |+ 88,345: |+100,508% |-+113,4754 |+127,5525 |+142,338i 
06 11,608 |— 9,232 |— 6,942 |— 4,614 |— 2,219 0,190 2,652 
|4109,038i |+126,423i |+145,303i |+165,543i |+187,1265 |4+210,496i \+234,980i 
08 45,424 48.584 51,925 55,470 | 59,225 | 63,200 | 67,405 
© |+222,673: |4+256,222i |+292,189i +330,553i |+372,018i |+415,397 
10 139,450 | 143,902 | 148,977 | 154,527 | 160,526 | 167,083 | 174,160 
|+328,076: |+381,356: |+439,088i |+500,9845 -+567,0143 +638,310: |+712,843% 


0,2 
0,4 
0,6 
0,8 
1,0 


0,0 
0,2 
0,4 
0,6 
0,8 
1,0 


4 [m 
— 
0,0 
0,2 
0,4 
0,6 
0,8 
1,0 
q 
0,0 
0,2 
0,4 
0,6 
0,8 
| as 
00 
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81 
Tabelle 2 (Fortsetzung) 
Afmu]| 400 | 460 500 | 550 | 600 | 650 | 700 
x! 

0.0 126,089 | 123,088 | 120,311 | 117,588 | 114,725 | 111,984 | 109,250 
|-115,8474 |—194,492¢ |—155,1825 —177,3675 _—201,044i |—226,526i 253,098: 
7 125,820 | 122,817 | 120,038 | 117,261 | 114,442 | 111,695 | 108,953 
 |_115,736¢ |—134,942¢ |—155,6975 |\—177,952i |—201,7045 253,924: 
4 123,940 | 120,918 | 118,125 | 115,318 | 112,464 | 109,671 | 106,876 
118,463: |—138,0885 |—159,3015 |—182,0475 |—206,322i |—259,711¢ 
os 118,836 | 115,766 | 112,933 | 110,045 | 107,094 | 104,179 | 101,239 
|_425,864¢ |\—146,628% |—169,0845 |—193,162i |—218,856i |—246,534¢ 275,416: 
as 108,896 | 105,731 | 102,822 99,776 96,637 93,483 90,263 
1 140,2775 |—163,259% |—188,135i |—214,807¢ |—243,265¢ 273,952: |—305,999% 
in 92,509 89,189 86,153 82,847 79,397 75,849 72,166 
-164,0385 |—190,6785 |—219,5445 |—250,492: |—283,507: |—319,155¢ |—356,421¢ 
q 

09 14,917 14,357 13,897 13,467 | 13,054 | 12,689 12,351 
13,6464 |— 15,688i 17,924i |— 20,3217 |— 22,8754 — 25,668% — 28,6144 
03 12,419 11,652 10,913 10,136 ' 9,300 | 8,419 7,462 
% 15,9474 |— 18,308¢ |— 20,879: |— 23,6314 |\— 26,554i — 29,7434 33,0968 
5,947 4,290 2,396 0,188 |— 2,390 8,874 
* |_ 19,7405 |— 22,7305 |— 25,5975 28,8214 32,2055 — 35,864: 39,6615 
06 1.6899. |— 11,562 |— 17,267 |— 24,185 |— 32,491 — 42,411 — 54,162 
19,8025 22,449: |— 24,849: |— 27,4554 30,0475 — 32,718: 35,291: 
08 37,267 |— 50,774 |— 67,620 |— 88,287 |—113,318 —143,413 |—179,230 
91555 |— 8,189 — 6,516 — 4,1965 — 0,970: 
10 {109,427 |—144,775 |—189,099 |—243.671 |—309,950 |—389,779 484,903 
“14 18,9854 |+ 25,803: |+ 34,3444 |+ 44,858: |+ 57,6244 |+ 72,9963 |+ 91,3124 

x 
00 31410. |— 30,318 |— 29,403 |— 28,540 |— 27,708 |— 26,950 — 26,246 
28,7345 |+ 33,1285 |+ 37,9257 |+ 43,0675 |+ 48,546: |+ 54,5174 |+ 60,8038 
0a 24787 |— 23,098 |— 21,477 |— 19,746 |— 17,859 |— 15,842 |— 13,639 
1+'35,4065 40,812: |+ 46,6927 |+ 52,976: |+ 59,651 + 66,8935 + 74,4796 
04 „6046 |— 2,720 1,057 5,442 10,539 16,416 23,198 
58,2125 |-+ 61,3205 |+ 70,0817 |+ 79,393: |+ 89,2295 99,8195 |+110,8116 
06 23,358 29,744 37,412 46,580 57,462 70,295 85,331 
78,7415 |+ 90,705: |+103,548: |+117,1275 |+131,3895 |+146,628i |+162,296: 
08 63,694 75,132 89,192 | 106,213 | 126,598 | 150,860 | 179,476 
14.408,369¢ |+124,727¢ |+142,1645 |+160,484: |+179,595i |+199,836i |+220,415i 
10 117,466 | 137,244 | 161,589 | 191,373 | 227,212 | 270,077 | 320,823 
|4-138,142¢ |+158,709¢ |+180,417i |+203,014i |+226,345¢ |+250,7415 |4+275,123% 

00 | 3876. |— 3,723 |— 3,589 |— 3,460 |— 3,333 |— 3,214 |~ 3,099 
4,8695 5,6195 |+ 6,4865 |+ 7,3125 |+ 8,2465 |+ 9,263¢ |+ 10,3326 
0a — 3871. |— 3,718 |— 3,584 |— 3,455 |— 3,327 |— 3,208 |— 3,093 
|4 |+ 5,627¢ 6.4467 |+ 7,324¢ 8,259 |+ 9,2776 |-+ 10,3481 
04 | 3717 |— 3,562 |— 3,425 |— 3,293 |— 3,162 |— 3,038 |— 2,918 
4 5.0995 |+ 5,8865 |+ 6,7455 |+ 7,665: |+ 8.646: |+ 9,7135 |+ 10,885: 
2671 |— 23501 |— 2,345 |— 2,190 |— 2,083 1,881 |— 1,729. 
6,616 |+ 7,6585 |+ 8,7815 + 9,9905 |+ 11,279% |+ 12,679 |+ 14,148: 
08 1,202 1,427 1,655 1,893 2,143 2,402 2,672 
12,232: |4+ 14,1915 |+ 16,317: |+ 18,5975 |+ 21,027¢ |+ 23,658: |+ 26,412; 
10 11,621 11,992 12,415 12,877 | 183,377 13,924 14,513 
297.3404 |+ 31,780: |+ 36,5915 |+ 41,749: |+ 47,2515 |+ 53,193: |+ 59.4047 


Ann. Physik. 6. Folge, Bd. 11 


nge- 
ings 
ihn- 
effi- 
igen 
700 
‚64048 
145% | 
‚133% 
‚8807; 
‚00519 
‚021% 
‚6749 
‚6010: | 
‚0478 
‚218% | 
‚7694 | 
‚1308 | 
‚3% | 
der 
sf fi- 
en) 
700 
4,600 
3,3231 
2,850 
3,2005 
),599 
2,3331 
2,652 
1,9805 
7,405 
3,3974 
1,160 
2,8431 | 


82 


Annalen der Physik. 6. Folge. Band 11. 1952 


Tabelle 2 (Fortsetzung) 


Alma] | 400 450 | So | 550 | 600 | 650 700 
X 
00 | 31897.) 30,94 | 30,041 29,128 | 28,194 | 27,268 | 26,388 
40,0721 — 46,704i |— 53,8731 — 61,5605 — 69,7631 — 78,5961 — 87,810i 
02 31,894 | 30,940 | 30,037 29,125 & - 28,191 | 27,265 | 26,330 
40,0775 46,7091 — 53,8791 — 61,5671 |— 69,7714 — 78,6051 — 87,820i 
04 31,794 | 30,840 | 29,936 29,021 | 28,085 | 27,157 26,219 
40,2214 -— 46,876: |— 54,0701 — 61,7851 — 70,017 — 78,8811 — 88,127 
06 31,114 | 30,153 © 29,244 28,319 | 27,370 | 26,426 | 25,468 
41,208% |— 48,0141 — 55,3741 — 63,2665 |— 71,6871 — 80,7571 — 90,220 
08 28,596 | 27,611 | 26,683 | 25,718 24,721 3 23,716 | 22,688 
44,8581 |— 52,2271 — 60,2001 — 68,748 77,8701 — 87,025 — 97,6% 
10 21,823 | 20,774 19,793 18,721 17,596 | 16,428 | 15,208 
054.679: |— 63,5591 — 73,181¢ — 83,4971 — 94,502: —106,3855 —118,80% 


Aus diesen Werten ergeben sich nun die in der folgenden Tabelle 3 zusammen- 
gestellten Hilfsgrößen und Amplitudenkoeffizienten: 


Tabelle 3 
Werte der Hilfsgröße 6 und der Amplitudenkoeffizienten p 
Almu]| 400 450 50 | 550 | 600 650 700 
q Re 6’ nach (33) 
0,0 0,03115 0,03320 0,03411 0,03509 | 0,03593 | 0,03644 0,03671 
0,2 0,03003 0,02824 0,02661 0.025022 | 0,02360 | 0,02217 | 0,02085. 
0,4 0,01581 —0,00548 —0,02482 —0,04213 |—0,05699 —0,06964 07991 
0,6 —0,06521 —0,11242 —0,15633 —0,19257 |—0,22222 —0,24680 —0,26603 
0,8 —0,17075 —0,23956 —0,29665 —0,34320 |—0,37992 —0,41044 10, 43305 
1,0 —0,26147 —0,34052 —0,40095 —0,44756 |—0,48300 —0,51072 Lo; 53048 
Im 6’ nach (33) 
0,0 0 0 . 0 0 0 0 0 
0,2 0.01781 0,01863 0,01869 0,01832 0,01765 0,01678 0,01583 
0,4 0,13261 0.13540 0,13424 0,13016 0,12417 | 0,11799 | 0,10969 
0,6 0,32873 0,32195 | 0,30914 0,29226 0,27334  0,25384 | 0,23500 
0,8 0,48481 0,46405 | 0,43654 0,40530 0,37406 0,34338 | 0,31509 
1,0 0,60029 0,55533 0,50725 0,46093 0,41776 0,37821 | 0,34336 
—6 Re pi nach (33) 
0,0 | 0,84426 | 0,83553 0,82951 | 0,82456 0,82035 0,81782 0,81647 
0,2. | 0,84837 0,84440 0,842.76 | ¢ ),84170 0,84082 0.84109 0,84201 
0,4 | 0,88145 | 0,91058 0, 93893: | 0,96415 | 0,98548 | 1,00447 | 1,02027 
0,6 | 1,00812 | 1,12349 1,22644 1,31458 | 1,38768 1,44836 1,49701 
08 | 1,39204 | 1,65024 1,86504 | 2,04037 | 2,18016 | 2,29164 2,37838 
1,0 | 2,30734 2.70269 3,00473 3,23778 3,41502 3,55062 3,65240 
6 Im pi = e, nach (33) s. Tabelle 4 
—6Im pi nach (33) 

0,0 0 o | 0 0 0 0 
0,2 0,01423 0,01490 | 0,01497 0,01471 | 0,01421 0,01355 0,01283 
0,4 | 0,11322 0,12456 | 0,13174 0,13481 0,13433 | 0,13134 0,12663 
0,6 | 0,38144 0,44728 | 0,49436 0,51922 0,52569 | 0,51908 0,50331 
0,8 | 0,96997 1,16440 1,28029 1,32946 1,32953 | 1,29724 1,24400 
1,0 | 2,17323 2,44535 . 2,54106 2,52383 .| 2,43558 2,30858 2,16479 
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Tabelle 3 (Fortsetzung) 
[mp] | 400 450 50 | 550 | 600 | 650 700 - 
— 6 Im pj’ nach (41) 
0,0 | 0 | 0 0 IR. 0 0 0 
02 | 0,00129 0,00151 0,00175 | 0,00200 0,00227  0,00256 0,00286 
0,4 0,00414 0,00484 0,00560 0,00641 0,00727  0,00820 0,00917 
0,6 0,03145 0,03677 0,04251 0,04866 | 0,05523 0,06228 0,06962 
0,8 0,13253 0,15500 0,17913 0,20506 | 0,23274 0,26245 0,29336 
1,0 0,40445 0,47289 0,54667 0,62580 | 0,71028 | 0,80094 0,89527 
10 Im p% nach (35) : 
0,0 0 0 | 0 0 | 0 0 0 
0,2 0,00010 0,00010 | 0,00009 0,00009 | 0,00008 0,00007 0,00007 
0,4 0,00320 0,00309 | 0,00292 0,00272 | 0,00251 , 0,00231 0,00212 
0,6 | 0,02396 | 0,02296 | 0,02158 | 0,02003 | °0,01845 | 0,01691  0,01550 
0,8 0,09510 0,08908 0,08226 | 0,07529 | 0,06858 | 0,06233 0,05677 
1,0 0,23697 0,21207 | 0,18940 | 0,16905 0.15119 0,13549 0,12211 


Aus diesen Werten ergeben sich schließlich nach (43) die in der folgenden Ta- 
‘belle 4 verzeichneten Koeffizienten des Extinktionsfaktors E: 


®) Koeffizienten nach Schoenberg und Jung‘), umgerechnet nach (55). 


Tabelle 4 
Die Koeffizienten des Extinktionsfaktors E nach (43) 
A{mu}| 400 450 500 550 600 | 650 700 
q e 
0,0 0 0 0 0 0 | 0 0 
0,2 0,02914 | 0,02956 0,02920 0,02826 0,02693 | 0,02539 0,02379 
0,4 0,23177 | 0,23365 0,22956 0,22109 0,20986 | 0,19713 0,18422 
06 | 0,76737 | 0,76110 0,73703 0,70115 0,65876 | 0,61358 0,56954 
0,8 1,72827 | 1,66392 1,57130 1,46420 1,35299 | 1,24339 1.14217 
1,0 3,01438- | 2,77355 2,53626 2,30463 2,08882 | 1,89084 1,71678 
1,0°) 3,00 2,70 2,54 2,30 2.18 1,92 _ 
e3 
00 | 0 0 0 0 0 | 0 0 
02 | 0,01563 0,01651 0,01681 0,01680 0,01656 . 0,01619 0,01576 
04 | 0,12057 0,13249 0,14025 0,14394 0,14412 0,14185 | 0,13792 
0.6 | 0,43684 0,50700 0,55846 0,58792 0,59937 0,59827 | 0,58843 
08 | 1,19761-| 1,40844 1,54168 1,60981 1,63086 1,62203 | 1,59413 
10 | 2,81466 | 3,13031 3,27712 3,31868 3,29705 3,24502 | 3,18217 
1,0°) 2,78 | 3,18 3,38 3,48 | 3,42 3,28 | — 
&4 
0,0 0,11880 0,11635 0,11468 0,11332 | 0,11216 0,11147 0,11110 
0,2 0,11981 0,11869 | 0,11823 0,11794 | 0,11771 0,11780 0,11807 
0,4 0,12054 0,12909 | 0,13815 0,14678 | 0,15452 0,16168 0,16784 
0,6 0,07124 0,11383 0,16015 0,20609 | 0,24862 0,28688 0,31944 
0,8  —0,17486 —0,00756 0,16823 0,33654 0,48709 0,61760 | 0,72536 
10 (—0,62711 —0,06468 | 0,43264 0,86198 | 1,21654 | 1,50527 1,73209 
1,0°) —0,60 0,06 0,48 0,78 1,08 1,36 — 
es 
1,0°) |—1,02 |\—0.36 0,04 | 0,28 0,52 0,74 | _ 
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Zum Vergleich sind in Tabelle 4 für g = 1 (homogene Eisenteilchen) außer den 
hier berechneten Koeffizienten auch die sich aus den Berechnungen von Schoen- 


berg und Jung®) ergebenden Werte verzeichnet. Wie man aus dem Zusammen- 


hang zwischen dem Extinktionsfaktor # und der von diesen Verfassern angegebenen 
interstellaren Extinktion pro kpe (Kxpec) 

d 
0,462 [1] 
für den Teilchendurchmesser d sofort findet, sind die Koeffizienten e,, e,, e, der 
Reihe nach einfach das Doppelte der von den genannten Verfassern berechneten 


Werte —A,, 4 A,, A, und ebenso kann der gleichen Untersuchung nach 


és = 2-(—2A4,) fiir g = 1 auch noch der Koeffizient des hier nicht mehr berück- 
sichtigten Glieles 5. Ordnung entnommen werden. Die Abweichungen der hier 
berechneten Koeffizienten von ihren von Schoenberg und Jung angegebenen 
Werte. sind i.a. besonders in Anbetracht ihrer sehr starken Abhängigkeit vom 
Brechungsindex, wie sie in der spektralen Variation zum Ausdruck kommt, zwar 
nicht erheblich, aber doch nicht ganz zu übergehen. Sie werden sich zum größten 


Kipe (55) 


Teil dadurch erklären lassen, daß die genannten Verfasser nicht bereits den von. 


W. Meier angegebenen spektralen Verlauf des Brechungsindex geglättet haben, 
wie es hier geschah, sondern erst die damit berechneten Koeffizienten. Da deren 
ursprüngliche Werte leider nicht angegeben wurden, ist ein strenger Vergleich 


natürlich nicht möglich, und es hätte auch nichts genützt, von genau den gleichen ' 


Werten des Brechungsindex auszugehen. 


3. Die Werte des Extinktionsfaktors 


Bevor man die in der nachfolgenden Tabelle 5 zusammengestellten Ender- 
gebnisse vergleicht, betrachte man die Abb. 3a—c, in denen der mit den Koeffi- 
zienten der Tabelle 4 nach (43) berechnete Extinktionsfaktor E für die Wellen- 
längen 400, 550 und 700 mu und verschiedene Radienverhältnisse g als Funktion 
von x dargestellt ist. Die Pfeile an den Kurven verweisen auf die Abszissenwerte £, 
bei denen bereits einer der Werte m,x und m, -q x die Zahl 1 überschreitet. 
Für g = 1 (homogene Eisenteilchen) konnte auch der Verlauf der Extinktions- 
kurve eingetragen werden, wie er sich nach 


E =e. ImnZ(—a,+ ,) ” 


(56) 
aus den von Schalen!) a entsprechend unseren Gleichungen (21) be- 
rechneten Werten Im --- durch Interpolation gewinnen läßt. Wie man sieht, 
bleiben für g = 1 die hier berechneten Näherungen bis zu x-Werten in der Nähe 
der in Abb. 3a—c gekennzeichneten £ und bei A = 700 mu sogar noch erheblich 
weiter durchaus brauchbar. Ist nun g <1, so zeigt eine Durchsicht der Rech- 
nung, daß davon nur eine Verbesserung der Konvergenz zu erwarten ist. So werden 
z.B. die bei der Herleitung der Gleichung (33) als klein gegen 1 behandelten 
Glieder pj| -x*, pı|-x° und 2%, -x°, von denen die ersten beiden mit den in den 
Tabellen 3 und 4 angegebenen Werten für g = 1 bei x = 1,2, der letzte erst bei 
x — 1,5 und darüber den Wert 1 erreicht, mit abnehmendem g schnell kleiner. 


8) C. Schalén, Nova Acta Reg. Soc. Sci. Upsaliensis (4) 10, Nr. 1 RM“ Medd. 
Astronom. Obs. Uppsala Nr. 64. 
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Die gefährliche Konvergenzschranke rückt also offenbar weiter fort. Hiernach 
darf man erwarten, daß auch für kleinere q die berechneten Näherungen jedenfalls 
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jeweils für den in den Diagrammen gekennzeichneten Bereich m, 2,9: m, x <1 
brauchbar sein werden. 

In Abb. 3a, b sind weiterhin noch die Ergebnisse eingetragen, die man durch 
Hinzunahme des nächsten Gliedes e, - x° in der Entwicklung des Extinktions- 
faktors für g = 1 erhält, wenn man für e, die in Tabelle 4 angegebenen auf die 
Berechnungen von Schoenberg und Jung zurückgehenden Werte benutzt. 
Offensichtlich hätte die Berechnung dieses Gliedes, die bereits für g = 1 wesentlich 
mühevoller ist als die der vorhergehenden Glieder, keine merkliche Verbesserung 
gebracht. 

Es folgt nun in Tabelle 5 eine Gesamtübersicht über alle hier berechneten 
Extinktionsfaktoren. Wo einer der Werte m, - x oder q - |m,| - x die Zahl 1 über- 
schreitet, wurde der entsprechende Extinktionsfaktor in Klammern eingeschlossen. 


Iga 
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(Schoen Schoenberg u. Jung (x) a= 75 my 
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Abb. 4. Spektrale Variation des Extinktionsfaktors für den Teilchendurchmesser 
75 mu und verschiedene Radienverhältnisse q 


Für g = 1 werden auch noch die in der beschriebenen Weise durch Hinzunahme 
des Gliedes 5. Ordnung nach Schoenberg und Jung yerbesserten Werte an- 
gegeben sowie auch die oben erwähnten interpolierten Werte nach den Berech- 
nungen von Schalén ohne Entwicklung der Amplituden. 
Zur Erleichterung der Diskussion dieser Ergebnisse dienen die Figuren 4, 5a, 
5b und 6, die zwar immer nur einen Ausschnitt aus der Tabelle 5 darstellen, 
aber doch alles Wesentliche deutlich erkennen lassen. Zunächst zeigt die Abb. 4 
für einen Teilchendurchmesser d = 75 mu |g E über.lg A. Man erkennt schon aus 
dieser Abbildung, daß mit zunehmendem Kernanteil der Übergang vom A- 
Gesetz der Rayleigh-Streuung zu den Verfärbungsgesetzen absorbierender 
Teilchen vorwiegend schon bei sehr kleinen Kernen stattfindet. Im übrigen zeigen 
diese Verfärbungsgesetze von der charakteristischen Rotverfärbung durch kleine 
Teilchen abgesehen keine bemerkenswerten Grundzüge, die nicht entscheidend 
von den Brechungsindizes abhängen, so daß sich ein weiteres Eingehen auf die 
Verfärbungskurven nach Art der Abb. 4 nicht lohnt. 
Anders die Abhängigkeit des Extinktionsfaktors vom Radienverhiltnis q, 
die für die Wellenlängen 400 und 700 mu und die Teilchendurchmesser 10, 25, 
50, 75 und 100 mu in den Abbildungen 5a und 5b gezeigt wird. Man erkennt 
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einerseits, wie bei konstantem g der Extinktionsfaktor in dem hier untersuchten 
Wellenlängen- und Größenbereich ganz roh proportional dem Durchmesser an- 


Tabelle 5 
Extinktionsfaktoren E für Wassertröpfehen mit Eisenkern 
[mu] | | 400 | 450 500 550 600 | 650 700 
q d= 10 
0,0 | 10-8-| 4,520 | 2,769 | 1,787 1,206 0,843 0,608 6,451 
0,2 10-3. | 2,300 | 2,071 1,840 1,618 1,413 1,229 1,069 
0,4 | 10-®-| 1,827 | 1,636 | 1,446 1,266 1.101 0,954 | 0,828 
0,6 10-2.) 6,048 | 5,331 | 4,645 4,016 3,458 | 2,972 | 2,562 
0,8 | 10-2-| 1,363 | 1,166 0,991 0,839 0,711 | 0,603 | 0,514 
1,0 | 10-1-| 2,381 | 1,947 | 1,602 1,323 1,098 0,918 0,773 
1,0%)| 10-1.) 2,381 1,947 | 1,602 1,323 1,098 | 0,918 _ 
d = 25 mu 
00 | 10-4 1,766 | 1,080 | 0,698 | 0,471 | 0,329 | 0,238 0,176 
0,2 | 10=. 6,018 | 5,356 | 4,724 | 4,183 | 3,697 | 3,121 | 2,701 
0,4 10. 4,729 | 4,160 | 3,669 | 3,205 | 2,784 | 2,410 2,089 
0,6 | 10-1.) 1,541 | 1,356 | 1,180 | 1.019 | 0,876 0,753 0,648 
0,8 | 10-1-| 3,482 | 2,979 | 2,529 | 2,139 | 1,809 1,532 1,305 
1,0 "| 10-2.| 6,122 | 5,007 | 4,114 | 3,896 | 2,812 2,345 1,974 
1,04) 10-1.| 6,119 | 5,006 | 4,114 | 3,396 | 2,812 2,345 ma 
d = 50 mu 
0,0 | 10-%-| 2,825 1,728 | 1,117 | 0,754 0,527 | 0,380 0,282 
02 |10-8-| 1,524 1,278 | 1,085 | 0,926 0,790 | 0,677 | 0,582 
0,4 | 10.) 1,012 | 0,891 0,778 | 0,675 | 0,583 | 0,502 0,433 
0,6 |101.| 3,295 | 2,889 | 2,504 2,153 | 1,844 1,577 1,353 
0,8 10. 7,471 | 6,406 | 5,431 | 4,579 | 3,858 3,255 2,762 
1,0 | 1,339 | 1,100 | 0,903 | 0,741 | 0,612 | 0,508 | 0,426 
1,02)| | 1,830 | 1,098 | 0,903 | 0,742 0,612 0,508 _ 
1.092), 130 | 1,08 | 088 | 0,75 08 
d = 75 mu 
0,0 rs 1,430 | 0,874 | 0,565 | 0,382 0,267 | 0,193 0,143 
0,2 | 10-%-| 3,478 | 2,677 | 2,185 1,740 | 1,438 | 1,201 1,012 
04 104.) 1,757 | 1,511 | 1,297 | 1,110 | 0,948 0,810 | 0,694 
0,6 | 101.) 5,499 | 4,798 | 4,137 3,535 | 3,009 | 2,559 | 2,183 
0,8 1,242 | 1,03 | 0,910 | 0,765 | 0,642 | 0,539 | 0,455 
1,0 (2,275) | 1,897 | 1,559 | 1,277 | 1,049 | 0,866 | 0,721 
1,011) (2,203) | 1,883 | 1,560 1,281 |. 1,054 0,870 | — 
1,022), 2,04 1,77 | 1,49 1,27 1,05 =e _ 
d = 100 mu 
0,0 |10-2.| 4,520 | 2,769 | 1,787 1,206 0,843 0,608 0,451 
02 | 10.) 7,604 | 5,444 | 4,094 3,183 2,533 2,053 | 1,689 
0,4 | 101.) 2,863 | 2,368 | 2,006 | 1,687 1,422 | 1,201 1,020 
0,6 10-1.) 8,414 | 7,309 | 6,016 | 5,320 4,497 3,798 3.218 
0,8 (1,871) | (1,639) 1300) | 1,172 0,979 0,818 | 0,686 
1,0 (3,493) | (2,986) | (2,474) | (2.027) | (1,658) | (1,362) | (1,128) 
1,022) (3,236) | (2,940) | (2,477) | (2,040) | (1,674) | (1, 377) _ 
1,022) 2,58 2,42 2,17 | 1,90 1,64 | se 


11) Durch Hinzunahme des Gliedes 5. Ordnung nach Schoenberg und Jung‘) ver- 
besserte Werte. 
12) Interpolierte Werte nach Berechnungen von Schal&n!®) ohne Entwicklung der 
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wachst, der Extinktionsquerschnitt also proportional zum Teilchenvolumen, wie 
zu erwarten, wenn d <4. Bemerkenswert ist jedoch die Abhängigkeit des Ex- 


+? T T T T 


4 4 b) 
0 02 04 06 08 1 


tors E von der relativen Größe des Kerns 
(Radienverhältnis q) für verschiedene Teilchendurchmesser d und die Wellenlängen 400 
und 760 mu. Die punktierten Kurven ( -- - )stellen den Verlauf nach der Näherungsglei- 

chung (50) dar 
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tinktionsfaktors von q: Es zeigt sich, daB der Ubergang der Extinktion von den 
Werten fiir g = 0 (homogenes Dielektrikum) zu denen für g = 1 (homogenes Ab- 
sorbens) für sehr kleine Teilchen im stärksten Maße bei kleinen Kernen (q < 0,2) 
stattfindet. Für größere Teilchen dagegen erstreckt sich dieser Übergang immer 
gleichmäßiger über den ganzen Bereich von g = 0,1 bis g= 1. Im Bereich der 
kleinsten Kerne (q <0,1) steigt mit wachsendem Kern E zunächst nur sehr langsam 
an, nämlich proportional g*, wie man leicht nachrechnet. Solche Kerne machen 
sich daher, wie zu erwarten, in derExtinktion praktisch überhaupt nicht bemerkbar. 
Dieser Bereich der ,,getarnten“ 
Kerne reicht bei einem Teilchen- 
durchmesser von 100 mu bis zu 
einem Kerndurchmesser von etwa 
10 mu, bei einem Teilchendurch- 
messer von 10 mu bietet er je- 
doch mit einem Kerndurchmesser 
unter 5Ä nur wenigen Atomen 
Raum. 

Die in Abb. 5a und 5b ein- 
getragenen punktierten Kurven 
wurden nach der Näherungs- 
formel (50) berechnet. Wie man 
sieht, ist die Annäherung im allge- 
meinen durchaus brauchbar, ob- 
wohl sich die Brechurgsindizes , 


bereits erheblich von 1 unter- 2 wa Te 7 
scheiden. Da bei A = 700 mu (yl 


dieser Unterschied für m, erheb- Abb. 6. Abhängigkeit der Wirkungsquerschnitte 


lich größer ist als bei A = 400 mu, 
ist die schlechtere Annäherung 
bei jener Wellenlänge durchaus 
verständlich. 

Schließlich zeigt mit Rück- 
sicht auf die meteorologischen 
‘Anwendungen Abb. 6 noch die 


vom Tropfendurchmesser d bei konstantem Kern 
für die Wellenlängen 400 ( ) und 700 mu 
(——-—). a: Extinktionsquerschnitt eines homo- 
genen Wassertrépfchens, 6: ds. für ein Wasser- 
tröpfchen mit Eisenkern von 10 mp Durchmesser, 
b’: zugehöriger Streuquerschnitt, c: Extinktions- 
querschnitt beieinem Eisenkern von 25 mu Durch- 
messer. Die oberen Kurven stellen das Verhältnis 
der Extinktionsquerschnitte für die beiden Wellen- 


Variation des Extinktionsquer- längen dar 

schnitts mit zunehmender rein 

dielektrischer Hülle bei konstantem absorbierenden Kern von 10 bzw. 25 mu 
Durchmesser. Wenn man von kleinen Nuancen absieht, deren Feststellung wegen 
ihrer starken Abhängigkeit von den Brechungsindizes praktisch ohne Wert ist, so 
ändert sich der Extinktionsquerschnitt bei zunehmender Aniagerung einer dielek- 
trischen Hülle an den absorbierenden Kern zunächst in einem sehr weiten Bereich 
nur außerordentlich wenig. Man könnte dieses Gebiet als den Bereich der ge- 
tarnten Hüllen bezeichnen. Erst wenn die Hülle so groß geworden ist, daß ein 
ebenso dickes homogenes dielektrisches Teilchen fast den gleichen Extinktions- 
querschnitt erreicht, beginnt auch der Extinktionsquerschnitt des Teilchens mit 
Kern merklich zuzunehmen und folgt schließlich nach einem augenscheinlich 
recht engen Übergangsgebiet dem wesentlich steileren Verlauf für das homogene 
dielektrische Teilchen. Dies ist jetzt der Bereich der getarnten Kerne. Nach den 
Überlegungen im letzten Abschnitt des ersten Teils ist dieses Verhalten physika- 


4 
wie 


90 


lisch leicht zu verstehen. Nach den gleichen Uberlegungen wird es aber auf kleine 
Teilchen beschränkt bleiben. Übrigens kommt das gleiche Verhalten in dem 
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ul 
in Abb. 6 miteingezeichneten Verlauf des Extinktionsverhiltnisses für A — 400 

und 700 mu wegen der oben iibergangenen, hier aber weit mehr ins Gewicht fallen- ur 

den Feinheiten nicht ganz so deutlich zum Ausdruck. B 

di 

4. Die Streuung 

Für die Wellenlängen 400 und 700 mu, d.h. die Kernbrechungsindizes m, = di 

1,174—1,291 2 und 1,747—1,921 7, wurden nun auch noch mit den gleichen Werten si 

des Radienverhältnisses q der Streufaktor S und die sphärische Albedo y berechnet. de 

Die folgende Tabelle zeigt jeweils die ersten drei nicht verschwindenden Koeffi- A 

Tabelle 6 ps 

Die Koeffizienten des Streufaktors und der sphärischeu Albedo nach (44) bzw. (48) 2 

| 6 % | Ya Ys | Ye 

2 = 400 mu a 

0,0 | 0,11880  —0,00740 0 _ d 

02 | 0,12010 —0,00721 —0,00117 4,1215 —0,7838 | —4,1518 ti 

0,4 | 0,13845 | —0,00438 —0,01068 | 0,5974 | —0,5822 | —0,5662 h 

0,6 0,26586 0,03467 —0,06800 | 0,3465 —0,5241 —0,1815 d 
0,8 0,82078 0,28030 —0,47284 | 0,4749 —0,5308 —0,1724 

1,0 2,60172 1,36053 —2,61418 | 0,8631 —0,4842 —0,6592 n 

A = 700 mu | 

0,0 | 0,11110 | —0,00816 0 | 

0,2 | 0,11826 —0,00493 —0,00094 | 4,9713 | —0,8700 —5,0028 fi 

0,4 0,17915 0,02863 —0,01100 | 0,9725 | —0,5932 —0,9708 t 

0,6 0,42757 0,22749 | —0,08117 | 0,7507 | —0,6337 —0,7034 f 
0,8 1,16021 | 1,00485 —0,44172 1,0158 —0,5159 —1,0218 

1,0 | 2,71456 |  2,88003 —1,55344 1,5812  —0,1760 —1,9138 

-] T s 

t 

f 


75 a2 [7 06 08 1 
Abb. 7. Streufaktor S für verschiedene Werte des Radienverhältnisses g bei 400 ( —) 
und 700 mu (— — — ) und verschiedenen Teilchendurchmessern d. Die punktierten Kurven 


(+++) stellen den Verlauf nach der Näherungsgleichung (51) dar 
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zienten dieser Größen, wie sie sich nach (44) und (48) aus den Werten der Tabellen 3 
und 4 ergeben. 

Die Reihe (48) für y konvergiert schlechter als die Reihen (43) und (44) für £ 
und S. Daher wurde hier die 


Berechnung „ auf Teilchen- ya (Fe) 
durchmesser d bis 50 mu be- N 

schränkt. Das Ergebnis zeigen N a a 
die Abb. 7 und 8. Während 7 AO a 

sich die Streuung mit wachsen- N u 5 } 50 


40% von d nur sehr wenig 
ändert, findet geradeindiesem 97 RN 


4 
dem Kerndurchmesser bis zu ls WKS. 


? 400 
Gebiet der stärkste Abfall der 700 
Albedo statt. Der kleine 400 
sehr großen Kernen rührt von + m mu 
dem stärker zunehmenden Bei- 
trag des Kerns zur Streuung “a 
her, hinter dem die Abnahme 0 02 04 06 08 1 
des Beitrags der Hülle nun- 
mehr zuriickbleibt. Abb. 8. Abhängigkeit der sphärischen Albedo y bei 


Die physikalische Erkli- 400 und ping 
rung des Verhaltens des Streu- messern d vom ienverhältnis g. O: berechnete, 


ten von q ergibt sich wieder aus der Überlegung, die zur Näherungsformel (51) 
führte. Der nach dieser Gleichung berechnete Verlauf wurde auch in Abb. 7 wieder 
punktiert eingezeichnet. Auch hier wird die Näherung (51) wieder in den meisten 
Fällen ausreichen. Was die Unter- ER 
schiede in der Güte der Annäherung 5 
bei beiden Wellenlängen anbetrifft, 
so gilt hier sichtlich das gleiche wie 
für die Extinktion (vgl. Abb.5a und b). 
Schließlich interessiert für die 
meteorologischen Anwendungen auch 
hier wieder die Änderung von Streu- 
querschnitt und Albedo bei zuneh- 
mender Wasseranlagerung an absor- 
bierenden Kernen gegebener Größe. 
Die erste ist für einen Kerndurch- 
messer von 10 mu in Abb. 6 mitein- a 4 1 f f 
gezeichnet. Sie ähnelt im wesent- 4 2 j 4 5 6 


J 
lichen dem Verhalten der Extinktion, Yo 


- Abb. 9. Sphärische Albedo y für Wasser- 
wenn auch das Hinbiogen der Streu trépfchen mit einem Eisenkern von 10 mu 


T TT 


„kurvein die eines homogenen Wasser- Durchmesser als Funktion des Tropfen- 


tröpfchens (S = E) bereits früher er- durchmessers d bzw. der durch die Kon- 
folgt. Die aus diesen Kurven in ohne densation bewirkten Teilchenvergrößerung 


1/q. Die für 400 ( ) und 700 mu (— ——) 
weiteres Weise graphisch seneichneten Kurven ergeben sich aus den 
zu entnehmenden Werte von Igy= Kurven der Abb. 6, die eingetragenen Ein- 


lg (SQ)—Ig (EQ) zeigt die Abb. 9 zelwerte wurden nach (48) berechnet 
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als Funktion der durch Wasseradsorption erfolgten Vergrößerung g-' des 
Teilchens — wieder für die Wellenlängen 400 und 700 mu. Die in dieser Ab- 
bildung eingetragenen Einzelwerte wurden aus (48) mit den Koeffizienten der 
Tabelle 6 berechnet. Die Übereinstimmung der auf beiden Wegen erhaltenen y- 
Werte ist durchaus befriedigend. 


III. Schlußbemerkungen 
über Anwendungen in Astrophysik und Meteorologie . 


1. Die interstellare Verfärbung des Sternenlichtes 


Die kontinuierliche Verfärbung des Sternenlichtes beim Durchgang durch - 
interstellare Materiewolken wurde bekanntlich erstmalig von Schoenberg und 
Jung*) wesentlich allein auf Streuung und Absorption durch Körner von beugender 
Größe zurückgeführt. Schoenberg und Jung')*), Schalen!) #2), Green- 
stein“) und van de Hulst!) haben versucht, unter verschiedenen Annahmen 
über die Materie und Größe der interstellaren Körner die beobachtete Verfärbung 
möglichst gut darzustellen, worüber von Schoenberg"*) und zuletzt von van 
de Hulst**) zusammenfassend berichtet wurde. Die sehr verschiedenen Ansichten 
der genannten Autoren über die Materie der Körner wurden kürzlich vom Ver- 
fasser‘”) kritisch gegeneinander abgewogen, der hierbei zu dem Ergebnis kam, 
daß zwar nicht metallische, wohl aber stark absorbierende Substanzen am Aufbau 
der Körner wahrscheinlich ganz wesentlich beteiligt sind. Allen genannten Unter- 
suchungen gemeinsam ist jedoch die Annahme einer homogenen Zusammensetzung 
der Körner, die zum Teil nur als eine versuchsweise Vereinfachung gedacht war, 
zum Teil aber auch nach den von Lindblad**), van de Hulst") und D. ter 
Haar?) entwickelten Vorstellungen über die Entstehung der Körner durch Konden- 
sation des interstellaren Gases folgendermaßen begründet wurde: 

1. Die Temperatur des Gases liegt außerordentlich hoch (um 10* ° K), die 
der bereits gebildeten Körner dagegen infolge ihrer Temperaturstrahlung so niedrig 
(unter 20° K), daß praktisch jedes ankommende kondensierbare Gasatom fest- 
gefroren wird und eine Oberflichenwanderung, die den Aufbau eines Kristall- 
gitters ermöglichen könnte, nicht stattfindet. 

2. Wegen der genannten hohen Gastemperatur ist die interstellare Materie 
von dem der Korntemperatur entsprechenden Kondensationsgleichgewicht weit 
entfernt. 

Auf diese Bedingungen für die Bildung homogener Körner wird der Verfasser 
in Kürze an anderer Stelle zurückkommen. Sie scheinen im allgemeinen durchaus 
nicht erfüllt zu sein, und man muß daher eine fraktionierte Kondensation an- 
nehmen, die zu einem schalenförmigen Aufbau der Körner aus Schichten ver- 
schiedener Zusammensetzung führen muß. Dabei reicht allein schon die Tat- 


18) C. Schalen, Nova Acta Reg. Soc. Sci. Upsaliensis (4) 11, Nr. 6 (1939) = Uppsala 
Astronom. Obs. Ann. 1, Nr. 2. 

14) J. L. Greenstein, Circ. Harvard Coll. Obs. Nr. 422 (1937). : 

15) H. C. van de Hulst, Rech. Astronom. Obs. Utrecht 11/2 (1949). 

16) E. Schoenberg u. H. Lambrecht, Erg. exakt. Naturwiss. 19, 1 (1940). 

1) A. Güttler, Z. Astrophysik 81, 1 (1952) = Mitt. Sternw. München I Nr. 5. 

18) B. Lindblad, Monthly Notices Roy. Astronom. Soc. 95, 20 (1934); Nature 
(London) 135, 133 (1925). 

12) H.C. van de Hulst, Nederl. Tijdschr. v. Natuurkunde 10, 251 (1943). 

20) D. ter Haar, Bull. Astronom. Inst. Netherl. 10, 1 (1943). 
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sache, daß unter den hieran vermutlich beteiligten Verbindungen die stark ab- 
sorbierenden wie die Oxyde des Eisens weit weniger flüchtig sind als die meisten 
der nicht absorbierenden (z. B. H,O, NH,, CH, neben dem viel weniger häufigen 
schwer flüchtigen SiO,), für eine Bildung abso bieiende: Ke.ne hin. 

Auch Cernuschi?!) glaubt aus anderen Gründen die bisherigen Vorstellungen 
über die Bildung der Kö.ner bestreiten zu müssen. Er veimutet, daß sich diese 
gegenwärtig in der Wiederauflösung befinden und ihre Bildung in einem ver- 
dichteten Frühstadium des Milchstraßensystems oder aus noch hinreichend dichter 
abgestoBener Sternmaterie (Nebelhüllen der Supernoven und Noven) stattge- 
funden haben muß. Jedoch wurden Cernuschis Einwände, in eister Linie das 
nach einer Theorie von Devonshire‘?) hergeleitete Verschwinden des für den 
Ausfrierprozeß maßgeblichen Akkomodationskoeffizienten für den Zusammenstoß 
der kondensierbaren Gasatome mit den Körnern bei sehr tiefen Korntempe:aturen, 
von ter Haar) und van de Hulst") unter Hinweis auf die Ergebnisse von 
Laboratoriumsversuchen von Langmuir’), von J. K. Roberts®) sowie von 
Keesom und G. Schmidt*) 2’) zurückgewiesen. Trotzdem scheint Cernuschis 
Gedanke einer Bildung der Kérner durch Kondensation eines wesentlich dichteren 
Gases, als es im gegenwirtigen interstellaren Raum angenommen werden muB, 
richtig zu sein, wie der Verfasser in Kiirze darzulegen beabsichtigt. 


Die Diskussion interstellarer Körner mit schalenförmiger Struktur und ab- 
sorbierendem Kern wird auch noch durch die gleichfalls beim Durchgang du ch 
interstellare Materiewolken von Hiltner‘), Hall, Markowitz und Mike- 
sell2%) 3°) 3!) beobachtete Polarisation des Sternenlichtes nahegelegt. Der Ve such 
einer Theorie dieser Erscheinung führte nämlich Spitzer und Tukey*) zu der 
Annahme von mit dielektrischer Materie überkrusteten gestreckten ferromagneti- 
schen Kernen, deren Entstehung sie nach einer Theorie von Oort und van de 
Hulst**) auf einen Ausschmelzvorgang bei Zusammenstößen interstellarer Wolken 
zurückführten. Als geeignete Ferromagnetika kommen nach diesen Verfassern 
gerade stark absorbierende Verbindungen in Betracht wie der hochopake Magnetit 
(Fe,0,), y-Hämatit (Fe,0,) und Magnesium-Eisen-Spinell (MgFe,0,). 

Leider führt nun der in dieser Weise nahegelegte Versuch, die vorliegende 
Theorie der Beugung durch dielektrische Kugeln mit absorbierendem Kein auf 
die interstellaren Körner anzuwenden, zu keinem befriedigenden Ergebnis. Zwar 
ist es leicht möglich, die beobachtete Verfärbung des Sternenlichtes beim Durch- 
gang durch interstellare Materie mittels der mit den hier benutzten Brechungs- 


21) F. Cernuschi, Astrophys. J. 105, 241 (1947) = Harvard Reprint Nr. 292. 
22) A. F. Devonshire, Proc. Roy. Soc. London (A) 158, 269 (1937). 
23) D. ter Haar, Astrophys. J. 106, 484 (1947). 
24) I. Langmuir, Phys. Rev. 8, 149 (1916). 
%) J. K. Roberts, Proc. Roy. Soc. London (A) 129, 146 (1930); 135, 192 (1932); 
142, 518 (1933); 152, 445 (1935). 
26) W. H. Keesom u. G. Schmidt, Physica 3, 590, 1085 (1936); 4, 828 (1937). 
27) G. Schmidt, Thesis Leiden 1938. 
” ph A. Hiltner, Astrophys. J. 109, 471 (1949); Science (Lancester Pa.) 109, 
165 (1949). 
29) J.S. Hall, Astronom. J. 54, 39, 187 (1948); Science (Lancester Pa.) 109, 166 (1949). 
30) J. S. Hallu. A. H. Mikesell, Publ. U. S. Naval Obs. Washington 17/1 (1950). 
81) W. Markowitz u. J. S. Hall, Astronom. J. 55, 175 (1950). 
3°) L. Spitzer jr. u. J. W. Tukey, Science (Lancester Pa.) 109, 461 (1949); Astro- 
phys. J. 114, 187 (1951). 
38) H. Oort u. H.C van de Hulst. Bull. Astronom. Inst. Netherl. 10, 187 (1946). 
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indizes berechneten Teilchen sehr gut wiederzugeben. Man kann dies nach einer 
in ihren Grundlagen auf van de Hulst*) zurückgehenden schnell fördernden 
graphischen Methode sofort zeigen, die in der schon genannten früheren Unter- 
suchung des Verfassers!”) bereits beschrieben und angewandt wurde, so daß wir 
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Abb. 10. Extinktionsfaktor E von Wassertröpf- 
chen mit Eisenkern bei verschiedenen Teil- 
chendurchmessern d (10 bzw. 25 mu) und Radien- mit Wassertrépfchen mit Eisen- 
verhaltnissen q als Funktion der zur gleichen kern einheitlicher Größe (d = 
Wellenlänge gehörenden scheinbaren Helligkeits- 10 mu, q = 0,2 bzw. 0,4) als Mo- 
differenz zwischen £ und e Persei nach Messun- dell-Rauchteilchen: Scheinbare 

gen von Stebbins und Whitford Helligkeitsdifferenz £—: Persei als 
Funktion von A-!. Große Kreise: 
Messungen von Stebbins und 


Abb. 11. Darstellung der inter- 
stellaren Verfärbung von £ Persei 


uns hier sehr kurz fassen können. Diese 
Methode beruht darauf; daß für den rich- 
tigen Extinktionsfaktor H(A“) zwischen 
diesem und dem beobachteten scheinbaren 
Helligkeitsunte:schied zwischen einem 
verfärbten und einem unverfärbten Stern von 
gleichem Spektraltyp als notwendige Be- 
dingung eine lineare Beziehung 


Whitford®°)3) bei der Messung 
bei 0,5 a”! von Whitford®®) ist 
der für diese von van de Hulst") 
abgeschätzte Fehler eingetragen. 
Kleine Kreise: ‘Zwischen den 
vorigen Werten interpolierte Wer- 
te. Kreuze: Aus Abb. 10 gefundene 
Werte. Die gestrichelten Teile der 
Kurven wurden willkürlich ergänzt 


Am (4) = Am(0) + - (57) 


gelten muß. Die im wesentlichen nur auf dem Entfernungsunterschied der beiden 
Sterne beruhende Konstante Am(0) sowie der von Zahl und Querschnitt der ab- 
sorbierenden und streuenden interstellaren Teilchen abhängige Koeffizient c er- 
geben sich aus der graphischen Da stellung Z (Am), sobald diese linear erscheint. 
In Abb. 10 ist dieses Verfahren für den veıfäi bten Stern £ Persei und den unver- 
färbten Vergleichsstern ¢ Persei, beide frühe B-Sterne, durchgeführt. Der Ein- 
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fluß eines kleinen spektralen Unterschiedes der beiden Sterne, der von L. Divan*) 
untersucht wurde, kann in diesem Zusammenhang übergangen werden. Von den 
erfolglos durchgeführten Hypothesen wurden nur die den biauchbaren benach- 
baıten eingezeichnet. Wie man aus dieser Abbildung ersieht, erhält man für sehr 
kleine Teilchen (d < 10 mu) und Radienverhältnisse q zwischen 0,2 und 0,4 schon 
unter der Annahme von Teilchen einheitlicher Größe brauchbare Darstellungen, 
und Abb. 11 zeigt, wie sich die Beobachtungen von Stebbins und Whitford®) 
im sichtbaren Spektralgebiet mit solchen Teilchen sehr gut erklären lassen, aber 
auch die Beobachtungen der gleichen Verfasser ®) %) im Infrarot einer derartigen 
Annahme nicht widersprechen. Für die spezifische Rötung 
Am (2,28 u!) — Am (0) 

Am (2,35 — Am (2,10 (6) 
erhält man für d = 10 mu und q = 0,2 bzw. 0,4 die Werte 7,0 und 6,8, also nicht 
viel weniger als den von Greenstein und Henyey*’) berechneten Mittelwert 
8,1 + 0,4. Die Werte der interstellaren Rauchdichte 9, erhält man mittels der 
aus einer Formel von van de Hulst*) leicht herzuleitenden Beziehung 

0, = 16-108 -d- —1)- (59) 


pg 8 em }’ 


worin die Dichten für Hülle und Kern s = 1,0 gem” und s, = 7,9 gem” und der 
photographische Extinktionskoeffizient apg = 1,2 mag/kpe nach McCuskey*) 
im Feld ZF 7 (gal. Länge: 133°) zu setzen und der photographische 
Extinktionsfaktor Hyg = E(2,28 uw!) Abb. 10 und 11 zu entnehmen ist. Es ergibt 
sich hier fiir q = 0,4 der durchaus plausible Wert 9, = 6 - 10% gem-*, während 
für q = 0,2 9, = 8- 105 gem~ allerdings wesentlich zu hoch erscheint, so daß 
ein Radienverhältnis in der Nähe von 0,4 den Vorzug verdient. 

Trotzdem kommen Teilchen der hier angegebenen Art für die Erklärung der 
interstellaren Verfärbung nicht in Betracht — nicht nur wegen der offensichtlichen 
chemischen Instabilität der Kombination Fe/H,O, sondern vor allem auch wegen 
ihrer viel zu geringen sphärischen Albedo. Aus Abb. 8 entnimmt man nämlich 
für A = 400 mu und d = 10 mu bei q = 0,2 den Wert y = 2-10-%, bei q = 0,4 
sogar nur 3 - 10*, also jedenfalls viel zu kleine Weıte, wenn man zum Vergleich 
die Bestimmung der photographischen Albedo der interstellaren Materie aus dem 
diffusen Licht der Milchstraße von Henyey und Greenstein®) heranzieht, 
welche 0,3 <y <0,8 fanden. Erst für x* = 0,1, d.h. dx 60 mu, würde man 
mit q = 0,2 bei den hier benutzten Brechungsindizes eine mit dem Ergebnis von 
Henyey und Greenstein vereinbare Albedo erhalten. Dann aber würde man 
zur Erklärung der interstellaren Verfärbung eine andere spektrale Dispersion der 
Brechungsindizes annehmen müssen. 

Durch die hier durchgeführte Berechnung wird erneut dargelegt, daß sehr 
kleine Teilchen für die optisch wirksamen interstellaren Körner kaum in Be- 
tracht kommen, solange man nicht die genannten und andere Albedountersuchun- 


34) L. Divan, C. R. Acad. Sci. Paris 233, 284 (1951). 
3) J. Stebbins u. A. E. Whitford, Astrophys. J. 98, 20 (1943); 102, 318 (1945). 
36) A. E. Whitford, Astrophys. J. 107, 102 (1948). 
37) J. L. Greenstein u. L. G. Henyey, Astrophys. J. 98, 327 (1941). 
3) S. W. McCuskey, Astrophys. J. 109, 414 (1949); s.a. B. J. Bok, Monthly Notes 
Astronom. Soc. South Africa 10, 51 (1951); Harvard-Coll. Bull. 920 (1951). 
3%) L. G. Henyey u. J. L. Greenstein, Astrophys. J. 98, 70 (1941). 


ler 

en 

PT- 

vir 

m 

’ 

sei 

n- . 

[0- 

ire 

als. 

se: 

nd 

ng 

ist 

15) 

on. 

en 

eT- 

ne 

ler 

zt 

7) 

en ” 

‚b- 

nt. 

in- 


96 Annalen der Physik. 6. Folge. Band 11. 1952 


gen in Anbetracht ihrer vomVerfasser bereitsfrüher auseinandergesetzten Schwierig- 
keit ganz in Zweifel ziehen will. Man wird die wirksamen Körner wohl eher unter 
den größeren Teilchen— jedoch mit anderen Brechungsindizes — und vielleicht 
sogar schon in dem Größenbereich suchen müssen, in den wir hier nicht eindringen 
konnten. Daiiiber hinaus — und das war ein wesentliches Ziel diese Exkurses — 
wird wieder gezeigt, wie wenig sich aus der beobachteten Verfärbung allein ent- 
nehmen läßt, wie wichtig dagegen weitere, zuverlässige Albedobeobachtungen sind. 


Nachdem der Verfasser schon in seiner früheren Untersuchung zu diesem 
Thema’) zeigte, daß es unmöglich ist, aus den gegenwärtig vorliegenden Beob- 
achtungsergebnissen zwingende Schlüsse über das spezielle Material homogen 
aufgebaut gedachter interstellarer Körner zu ziehen, wird diese Situation durch 
Einführung der Schalenstruktur zweifellos noch hoffnungsloser. Vorläufig wird 
man sich also damit begnügen müssen, auch diese Annahme als Arbeitshypothese 
zu benutzen. 


2. Die opaleszente Trübung der Atmosphäre 


Die hier behandelte Theorie zweischichtiger Teilchen würde für die Meteoro- 
logie von großem Interesse sein, da bekanntlich die Kondensation des Wasser- 
gehalts der Luft vorwiegend an Kondensationskernen stattfindet, den Aitken- 
kernen (Großionen bis gegen 0,2 u Durchmesser), Staubteilchen und hygroskopi- 
schen Tröpfchen (vgl. z. B. die von Moeller und Junge“) zusammengestellte 
tabellarische Übersicht). Jedoch fällt nur die Erscheinung der opaleszenten 
Trübung wenigstens zum großen Teil in das Gebiet der hier eingehender disku- 
tierten kleinen Teilchen, wobei hier die Kondensation durch die kleineren und 
mittleren Aitkenkerne eingeleitet wird. Bekanntlich zeigt sich die genännte Er- 
scheinung in einer Blauverfärbung des Landschaftsbildes, wie man sie z. B. be- 
sonders im Alpenvorland bei Föhn beobachtet. An diesem Phänomen ist neben der 
Rayleigh-Streuung der Luft jedenfalls beim antizyklonalen Föhn die Streuung 
durch die kleinen im Aerosol suspendierten Teilchen beteiligt, für die ebenfalls 
SA, während die größeren, stark, aber nur wenig selektiv streuenden fehlen. 
Dor Verfasser kann hier nicht über quantitative Beobachtungen der opaleszenten 
Trübung berichten, hofft aber, zu diesem Thema einen Beitrag beisteuern zu können 
durch einen Vergleich der optischen Eigenschaften derartiger Wassertrépfchen mit 
absorbierenden und nicht absorbierenden Kernen bei A = 500 mu. Als Materie 
des Kerns wurde dabei gewählt f 


I. Quarz, mittlerer Brechungsindex m = 1,55, 
II. Kalkspat, mittlerer Brechungsindex m = 1,58, 


III. Graphit, mittlerer Brechungsindex m = 3,02—1,76 7 nach Untersuchungen 
von Zakrzewski"), die mit einer Messung von v. Wartenberg**) gut überein- 
stimmen. (Es ist jedoch sehr wohl möglich, daß infolge eines durch die ausge- 
zeichnete Spaltbarkeit des Graphits nach der Basis verursachten Auswahleffekts 
dieser Wert nicht dem wahren Mittel über: alle Orientierungen entspricht.) 


40) F, Möller u. Ch. Junge in Landolt-Börnstein, Zahlenwerte und Funktionen, 
6. Aufl., Bd. 3, Berlin, Heidelberg, Göttingen 1952, S. 585. 

41) ©, Zakrzewski, Anz. Akad. Wiss. Krakau 1910 (A), 122. 

42) H.v. Wartenberg, Verh. dtsch. physik. Ges. 12, 105 (1910). 
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1V. Kohle, mittlerer Brechungsindex m = 1,94—0,66 2 nach Messungen von 
Senftleben und Benedict*) an reiner Bogenlampenkohle, deren Reflexions- 
vermögen mit Messungen von Aschkinaß“) an der abweichend zusammen- 
gesetzten Gaskohle annähernd übereinstimmte. 

Dabei wurden nun aber die Näherungsformeln (50) und (51) aes die nach 
den Feststellungen des vorhergehenden Teils dieser Untersuchung zweifellos fiir 
Berechnungen der vorliegenden Art ausreichen. Ferner wurde sinngemäß die Ent- 
wicklung der Wirkungsfaktoren Z und S für die homogenen Teilchen hinter dem 
ersten Gliede abgebrochen, so daß mit d als Teilchendurchmesser — hier 25 mu — 
und 


m? —1 
(60) 
und für absorbierende Stoffe 
E=—4ImM: Im M -d, (62) 


für nicht absorbierende dagegen E = S. Unter den gewählten Voraussetzungen 
erhält man dann die in der folgenden Tabelle 7 verzeichneten Werte dieser Wir- 
kungsfaktoren, denen zum Vergleich auch die sphärische Albedo y beigefügt wurde. 


Tabelle 7 
Wirkungsfaktoren der homogenen Teilchen bei 25 mu Durchmesser und 4 = 500 mu 
Kern | Stoff m s y 
— | Wasser 1,336 1 
I | Quarz | 1,55 | 6,00267 1 
- | Kelkepet | 1,58 | 0,00288 
III Graphit | 3,02—1,76 i | 0,00127 0,113 0,0112 
IV Kohle | 1,94—0,66 i 0,00056 0,138 | 0,0040 


Mit diesen Werten erhält man dann die in Abb. 12 dargestellte Abhängigkeit 
der Größen S, # und y vom Radienverhältnis q, die natürlich den in Abb. 5, 7 und 8 
gezeigten entsprechenden Zusammenhängen bei Eisenkernen ähnelt. Da nach 
Tabelle 7 Kalkspat sich hier nur außerordentlich wenig von Quarz unterscheidet, 
wurde der Kern II in der Abbildung weggelassen. Graphit und Kohle ergeben 
trotz der sehr verschiedenen Brechungsindizes einen merkwürdig ähnlichen Effekt. 
Dagegen ist der Unterschied der Extinktion bei relativ dicken absorbierenden und ~ 
nicht absorbierenden Kernen offensichtlich enorm. Kerne aus Graphit oder Kohle, 
die das halbe Teilchenvolumen einnehmen (q = 0,8), erhöhen die Extinktion auf 
das Tausendfache, während sich ein Kalkspat- oder Quarzkern bis zu dieser 
Größe noch so gut wie überhaupt nicht bemerkbar macht. Die Streuung 
zeigt ein ums anderes Verhalten und eine weit geringere Variation. 


is) H. Senftleben u. E. Benedict, Ann. Physik (4) 54, 65 (1917). 
44) E. Aschkinaß, Ann. Physik (4) 18, 373 (1905). 


Ann. Physik. 6. Folge, Bd. 11 


— 23 

ese 
TO- 
en- 
pi- 
Ite | 
ben 
nd 
be- 
ler 
ng 
lis 
on. 
en 
en 
rit 
rie 
on 
n- 
e- 
ts 
| 

7 


98 Annalen der Physik. 6. Folge. Band 11. 1952 


Es erübrigt sich, für diese absorbierenden Kerne nochmals die gleichen Be- 


trachtungen anzustellen, wie sie oben bei Eisenkernen an die Abb. 6 und 9 ange- 
knüpft wurden. Alle dortigen Fetsstellungen sind auch auf die hier besprochenen 
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Abb. 12. Die optischen Eigenschaften von Wasser- 
tröpfchen von 25 mu Durchmesser mit verschie- 
denen Kernen als Funktion des Radienverhält- 
i : Extinktionsfaktor, — — —: Streu- 


Kerne: Quarz (Q), Graphit (G) und Kohle (C) 


Kerne mutatis mutandis zu 
übertragen. Dagegen ist der 
Einfluß dielektrischer Kerne 
offensichtlich vergleichsweise ge- 
ring. Sobald die Wasserhüllen 
ungefähr das gleiche Volumen 
erreichen wie die Kerne, werden 
sich Teilchen mit dielektrischen 
Kernen, die bei der opaleszenten 
Trübung im «allgemeinen über- 
wiegen, von homogenen Wasser- 
tröpfchen nicht merklich unter- 
scheiden. Der auch in der Streu- 
ung noch beträchtliche Einfluß 
absorbierender Kerne stellt je- 
doch bei experimentellen Unter- 
suchungen der opaleszenten 
Trübung zweifellos eine gewisse 
Schwierigkeit dar, die auch nicht 
durch Messung der Verfärbung 
des Streulichts zu umgehen ist, 
da dieser in jedem Falle ein 2~- 
Gesetz zugrunde liegt. 


Den Herren Prof. Schoenberg und Prof. Bucerius ist der Verfasser für ihr 
dieser Untersuchung entgegengebrachtes Interesse und anregende Diskussionen 
zu Dank verpflichtet, den Herren Prof. R. Geiger und Dr. Roßmann für wert- 


volle Ratschläge zum letzten Abschnitt. 


München, Universitäts-Sternwarte. 


(Bei der Redaktion eingegangen am 16. Mai 1952.) 
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Theoretische Betrachtungen über die Beeinflussung 
der ferromagnetischen Koerzitivkraft durch Einschlüsse 
mit rotationselliptischer Form, für den Fall, daß deren 


Abmessungen klein gegen die Dicke der Blochwand sind 


Von Eberhard Schwabe 
(Mit 3 Abbildungen) 


Inhaltsübersicht 


Die Arbeit zeigt, daß für nichtkugelige Einschlüsse, deren Abmessungen klein 
gegen die Dicke der Blochwand sind, erhebliche Abweichungen gegenüber den 
bisher bekannten Abschätzungen für H, zu erwarten sind, sobald die Achsenab- 
messungen um mehr als 4% von der Kugelgestalt abweichen. Für den einfach 
zu behandelnden Fall, daß 180°-Wandverschiebungen über rotationselliptische 
Teilchen mit Hauptachsen der Länge L längs den Richtungen leichtester Magneti- 
sierbarkeit vorliegen, werden die zu erwartenden Abweichungen der bekannten 


Theorien als Funktion vom Dimensionsverhältnis < der Teilchen in einer Ab- 


bildung wiedergegeben. Danach liegen z. B. fiir 3 ay 3 wegen der Anisotropie 


der Streufeldenergie die zu erwartenden Koerzitivkräfte um etwa den Faktor 25 
höher als bisher angenommen wurde. 


Weiterhin sind Gleichungen zur Berechnung von H, für unmagnetische und 
magnetische Einschlüsse mit d < L< 6 bzw. L<d 5 unter Beriicksichtigung 
von 180°-Wandverschiebungen an 


Einleitung 


Seit der Entdeckung der Barkhausensprünge!) weiß man, daß die ur- 
sprüngliche Hypothese der nach ihrem Ergründer Weiß benannten Bezirks- 
struktur durchaus berechtigt ist. Diese Hypothese besagt, daß auch ein nicht- 
magnetisiertes Ferromagnetikum aus einzelnen, bezüglich der Magnetisierungs- 
richtung statistisch verteilten, bis zur Sättigung „spontan‘“ magnetisierten 
„Weißschen Bezirken‘ besteht. 

Nachdem Bloch?) die Berechnung der energetischen und geometrischen Ver- 
hältnisse der ,,Blochwinde“*) durchgeführt hatte, waren die Voraussetzungen 
für eine Theorie der Koerzitivkraft H, gegeben. 


1) H. Barkhausen, Z. Physik 20, 401 (1919). 

2) F. Bloch, Z. Physik 74, 295 (1932). 

3) Als Blochwände werden bekanntlich die Übergangsgebiete zwischen Weißschen 
Bezirken bezeichnet. 
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Die Wirkung innerer Spannungen auf H, beriicksichtigte R. Becker‘). 1942 
bestätigte Kersten®) mit einer anschaulichen Theorie, die bereits vorher von 
Gumlich, Köster, Gerlach, sowie Kußmann und Scharnow u.a. experi- 
mentell gefundene Tatsache einer Erhöhung der Koerzitivkraft durch Ausschei- 
dung von Fremdstoffen. Dabei zeigte er, daß nicht, wie von Kußmann und 
Scharnow®) angenommen, der Umweg über die eventuell mit der Ausscheidung 
verbundenen inneren Spannungen notwendig, sondern daß eine direkte Beein- 
flussung der Beweglichkeit einer Blochwand durch Änderung ihrer Oberflächen- 
energie an Fremdkörpern möglich ist. Die nach dieser Theorie gefundenen Werte 
stimmten in vielen Fällen mit den experimentellen Ergebnissen überein. 

Eine wesentliche Erweiterung der allgemeinen Theorie wurde durch die Er- 
kenntnisse von L. Néel’) gegeben. N&el erhob gegen die Kerstensche Theorie 
den Einwand, daß bei den Betrachtungen der Wandverschiebungen ein sehr 
wichtiger Term, die ,,Streufeldenergie“, nicht berücksichtigt wurde. Für Teilchen, 
deren Abmessungen groß gegen die Dicke der Blochwand sind (d > 6), wird 
dieser Einfluß durch die ebenfalls von N&el theoretisch vorausgesagten, später 
durch Williams und Bozorth®) experimentell bestätigten Zipfelmützen herab- 
gesetzt, liefert summarisch jedoch immer noch zu hohe Koerzitivkräfte. Erst 
durch die Einführung einer statistischen Verteilung der Störungen wurden die 
Ergebnisse mit dem Experiment größenordnungsmäßig wieder vergleichbar. 

Die Néelschen Ansätze sind sehr allgemein gefaßt und: beziehen sich sowohl 
auf Störungen durch innere Spannungen (im Sinne Beckers) als auch auf Sté- 
rungen durch Fremdkörpereinschlüsse, wie sie von M. Kersten unter Vernach- 
lässigung der Streufeldenergie Berücksichtigung fanden. 

Die vorliegende Arbeit beschränkt sich im Gegensatz zu N éel lediglich auf den 
Einfluß von Fremdkörpereinschlüssen. Dabei soll gezeigt werden, daß die geo- 
metrische Anisotropie von Einschlüssen wegen der richtungsabhängigen Streu- 
feldenergien gegebenenfalls recht erhebliche Abweichungen von den Ergebnissen 
der oben erwähnten Theorien zur Folge haben kann. Um den Kreis der Betrach- 
tungen nicht zu weit zu ziehen, werden in dieser Arbeit nur Teilchen mit der Form 
von Rotationsellipsoiden angenommen, bei denen die Abmessungen noch wesent- 
lich unter der Wanddicke 6 bleiben (d < 6). Weiterhin sollen sich die Abschätzungen 
lediglich auf 180°-Wandverschiebungen beschränken. 


Durchführung der Betrachtung 

Die Energiebilanz der 180°-Wand ist gegeben durch: 

2H,I,daV/ = 2H,I!,F dz = dEnax- (1) 
H, = Grenzfeldstärke 
I, = Sättigungsmagnetisierung 
dV = Volumenzunahme des zum äußeren Feld günstig liegenden Weiß- 

schen Bezirks bei einer Wandverschiebung um dx 
== Fläche der Wand (wegen d <6 gilt F = const.) 
dE nax = Maximaler Energiebedarf bei 180°-Wandverschiebung um dz. 


4) R. Becker u. W. Döring, „Ferromagnetismus‘“, J. Springer, Berlin 1939. 

5) M. Kersten, „Grundlagen einer Theorie der ferromagnetischen Hysterese“. 
Hirzel, Leipzig 1943. 

6) A. Kußmann u. B. Scharnow, Z. Physik 54, 1 (1929). 

7) L. Néel, Ann. Université Grenoble 22, 299 (1946). 

5) Williams Bozorth u. Shockley, Physic. Rev. 75, 155 (1949). 
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Aus Gl. (1) folgt der Ausdruck: 


dE 
(2) 
Bezeichnet man die Feldstärke, die zur Verschiebung der Wand in 2-Richtung 
über die ungiinstigste Stelle notwendig ist, als Grenzfeldstärke Ho, so ergibt sich 
aus Gl. (2) mit dem nach Becker*) angenommenen Wert: 


3 
H, N 
die Gleichung: 
dE 
3 i=) wax 
A,» (3} 


dx 
ob es exakt möglich ist, sämtliche Effekte zu erfassen, welche die Wandenergie- 
mulden über den Teilchen verflachen und über Bereiche von der Größenordnung 
der Wanddicke verschmieren. Jedoch kann man mit guter Näherung der weiteren 
Berechnung die folgende Gleichung zugrunde legen: 


dE 24E 
% \ (4) 


Dabei bedeutet AE die zu einer Wandverschiebung um = maximal benötigte 


Die genaue Berechnung des a) ist schwierig, und es erscheint fraglich, 
max 


2) 
Energie. Beachten wir nochmals die Bedingung d < 6, so ergibt sich praktisch 
für AE die Energiedifferenz eines Wandstückes mit Teilchenzahl r, und 
eines gleich großen unmittelbar angrenzenden mit Teilchenzahl n,. Kennt 
man die über mehrere Wände gemittelte Differenz Anmax = (Nz— N,)max, dann 
ist, wie man leicht einsieht: 


dE 
(Fe max N Anmax (5) 


AE, = Gesamtenergiedifferenz eines einzelnen Teilchens in- und auBerhalb 
der Wand. 
Die Koerzitivkraft folgt nach den Gin. (3) und (5) zu: 


3 AEB, Anmax 

(6) 
Der nächste Schritt wird die Ermittlung von AE sein. Es ist zweckmäßig, JE 
zu unterteilen in: 


H,x 


AE\wana = Energieänderung durch Volumenänderung der Wand. 
E\ streut. = entsprechende Änderung der Streufeldenergie. 


Es ergibt sich also: 
AE, N AE\ wana + AE, Streuf. - (7) 
Wir unterscheiden weiterhin: 


a) Geometrisch isotrope Teilchen, also Kugeln. 
b) Geometrisch anisotrope Teilchen, z. B. Nadeln oder Scheiben usw. 


: 
- 
1ei- 
ınd | 
Ing | 
in- 
en- 
rte 
Er- 
rie 
ehr 
en, 
ird 
ter 
ab- 
rst 
die 
ch- 
len 
eu- 
sen 
ch- 
rm 
gen 
| 
| 
iB- | 


102 Annalen der Physik. 6. Folge. Band 11. 1952 


a) Energieverhältnisse für kugelige unmagnetische Einzelteilchen mit d < 6 


In diesem Falle ist, wie N ée1*) feststellte, die Streufeldenergie sowohl innerhalb 
als auch außerhalb der Wand vorhanden und nur das Streufeld entsprechend den 
umliegenden Spins gedreht. Das Teilchen befindet sich jedoch wegen seines ge- 
ringen Durchmessers auch in der Wand praktisch in einem homogenen magneti- 
schen Fluß und wird somit nahezu die gleiche Streufeldenergie haben wie außer- 
halb. (Néel®) zeigte, daß z. B. bei Eisen mit kugeligen Einschlüssen vom Durch- 
messer d — 0,2ö die beiden Energieanteile AEjwana und AEistreur. bereits 
gleich sind, wobei sich letzterer aus den Quadrupolstreufeldern ergibt.) Für 
d <6 folgt aber, wie gesagt, daß die Differenz der Streufeldenergie in den beiden 
zu berücksichtigenden Stellungen der Blochwand — über und vor dem Teilchen — 
verschwindet. Also: 

AE | streut. = 0. (8) 


Der zweite in Frage kommende Term, die Einsparung von Wandenergie, liefert, 
falls die Wand über dem Teilchen liegt, bei unmagnetischen Einschlüssen den 
Beitrag: 


E\wana = (9) 
y = Flächendichte der Wandenergie, für Fe: Yıso & 1,7 = ; 
v = Volumen des Einzelteilchens. 
Außerhalb der Wand gilt ohnehin: 
E\wana = 0. (10) 
Also: 
Mit diesen Werten erhalten wir fiir die Koerzitivkraft nach der Gl. (6): 
(12) 


b) Energieverhältnisse für unmagnetisches Einzelteilchen mit der Gestalt 
eines Rotationsellipsoides, dessen Abmessungen klein gegen die Wanddicke sind 


Bei den ausgeprägt anisotropen Fremdkörpern ist zu bemerken, daß sie, soweit 
durch Ausscheidung im ungestörten Mutterkristall entstanden, womöglich keine 
statistische Verteilung ihrer Vorzugsrichtungen aufweisen. Dagegen sprächen die 
energetischen Verhältnisse bei der Entstehung der Ausscheidung. Es wird oft 
so sein, daß sich die langen Achsen von ausgeschiedenen Rotationsellipsoiden z. B. 
im kubischen Gitter nur längs dreier senkrecht aufeinander stehender Vorzugs- 
richtungen des Mutterkristalls bilden können. Anders liegen die Verhältnisse 
natürlich bei Schlackeeinschlüssen, die aus der Schmelze her eingefroren sind, 
und bei wenig von der Kugelgestalt abweichenden Fremdkörpeın. 

Wir werden hier nur den Fall der 180°-Wandverschiebungen über Rotations- 
ellipsoide durchrechnen, bei denen die Symmetrieachsen in den drei Richtungen 


%)'L. Néel, Cahier de Physique, Nr. 25—38 (Dez. 1944). 
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der leichtesten Magnetisierbarkeit liegen. Und zwar wird jeder einzelne der drei 
grundsätzlich möglichen, aus Abb. 1 zu ersehenden Fälle getrennt behandelt 
werden. 


Starkes schwaches starkes 
Streufeld Streufeld Streufeld 


Blochwand ochwand lochwand 


|| « 


Abb. 1. Die drei wichtigen Lagen der Einschlüsse bei Verschiebung einer 180°-Wand 


A 


Bezeichnen wir mit L die Länge der Symmetrieachse eines Einschlusses und 
mit d die Länge der dazu senkrechten, dann liefert der Fall A für AE\wana beim 
AbreiBen der Wand vom Einschluß ungefähr die Energiedifferenz : 


AE wand a (13) 


Für den Anteil der Streufeldenergie gilt, solange Z <0, sicher: 
AE streut. a = 0, (14) 


denn sowohl außer- als auch innerhalb der Wand wird das Teilchen unter der 
gemachten Voraussetzung senkrecht zu seiner Rotationsachse magnetisiert. 

Kommt dagegen L in die Größenordnung der Wanddicke 6, dann wird die 
Streufeldenergie innerhalb der Blochwand wegen des dort entstehenden Quadrupol- 
streufeldes geringer als außerhalb und Gl. (14) verliert ihre Gültigkeit. Doch 
können wir auch dann AKjstreur. a vernachlässigen, denn die entsprechenden 
Ausdrücke für AEjstreut.g bzw. AEjstreur.c Werden im Sinne der folgenden 
Betrachtungen bei der beabsichtigten Mittelwertsbildung aller 3 Größen wesent- 
lich stärker eingehen als ersteres. Im Falle B liegt nach Abb. 1 die Wand über 
dem Teilchen äußerst günstig. Durch die Einsparung an Wandenergie gilt: 


AE\ wana B at. (15) 


Auch die Streufeldenergie wird weitgehend vermieden, wenn das Teilchen in der 
Wandmitte liegt, und es gilt deshalb hier: 


N,— WN. 
AE | streut.B © I; v. (16) 


N, = Entmagnetisierungsfaktor quer zu L. 


N, = Entmagnetisierungsfaktor lings L. 
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Aus der Summation von Gl. (15) und (16) folgt für den in Abb. 1 gezeigten Fall 
länglicher Teilchen 


= (54 a (17) 
Wegen des bei Rotationsellipsoiden bestehenden Zusammenhanges: 
2N,+N,=42 (18) 


bzw. N,—N,=3N,—42 
vereinfacht sich Gl. (17) zu: 


Für den Fall C in Abb. 1 entsteht bei länglichen Rotationsellipsoiden die ent- 
scheidende Frage, wo die Wand energetisch günstiger liegt, vor oder über dem 
Teilchen. Was die Einsparung an Wandvolumen betrifft, so wird diese nach wie 
vor dafür sorgen, daß die Energie der Wand über dem Teilchen ein Minimum 
besitzt. Mit der Streufeldenergie sieht es jedoch in diesem Falle gerade anders aus. 
Sie tritt bei dem länglich angenommenen Teilchen in der Wand verstärkt auf 
und wird also die Wandstellung vor dem Teilchen begünstigen. Nach analoger 
Behandlung zu B ergibt sich: 


(20) 
Wir schätzen nun einmal für den Fall einer unmagnetischen langen Nadel mit 


I> dim «-Eisen die Werte von z und z J} ab, um danach zu entscheiden, welche 


Lage günstiger ist. (Für N,gilt:N,— F)- Dabei werden die folgenden Näherungs- 
werte für Eisen benutzt: 


erg 
cm? 


0,6 10-5 em 


> 


I, = 1700 (EME). 


> 


Man erhält: 


7=238.105 x I, = 0,9- 1074. 
Sehr ähnlich liegen die für das mit 
L<d. Nur sind hier die AE;streur. negativ, weil (3 N,—4r)<0. Da wir 
aber praktisch doch nur die Absolutbeträge von AE|3 bzw. AE,c benötigen, 
brauchen die Vorzeichen nicht beachtet zu werden. 

Die 180°-Wand wird also im Falle C, wenn L > d, vor einem Einschluß halt 


machen, weil in dieser Stellung die Energie um ungefähr den Faktor 55 gerin- 


ger ist als bei der um -; benachbarten Wandlage über dem Teilchen. 


Man könnte nun annehmen, daß der Term der Streufeldenergie bei regel- 
mäßiger Anordnung von Teilchen in Stellung B und C durch die verschiedenen 
Vorzeichen gerade kompensiert würde, denn z. B. für längliche Teilchen liegt die 
Wand in B über dem Fremdkörper günstig, während sie in € davor halt machen 
wird. (Bei abgeplatteten Teilchen würde es entsprechend umgekehrt sein.) Zur 
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Möglichkeit der Kompensation muß aber bemerkt werden, daß sie nur bei einer 
vollständig ebenen Blochwand besteht. Da vor derartig hohen Energiebergen 
jedoch sicher eine Durchbiegung der Wand auftritt, wird sich etwa eine Lage ein- 
stellen, wie sieschematisch aus Abb. 2zu ersehen 
ist. Dabei ist also Voraussetzung, daß die Streu- 
feldenergie der Einschlüsse die wesentlichen 
Beiträge zu bzw. AEjc liefert. Die 
Krümmung um das Teilchen B (Abb. 2) 
wird naturgemäß geringer sein als die um C, 
da die mit Energie benötigenden Oberflächen 
ladungen auf der Wand verbunden ist. Wir 
werden in unserer weiteren Betrachtung also 
sinngemäß lediglich die Absolutbeträge von 
bzw. überlagern. 


e) Einfluß der Teilchenverteilung 


Es müssen zwei Arten der Verteilung be- 
achtet werden, nämlich: 


Die nichtstatistische Verteilung. 
Die statistische Verteilung. 


Beide Fälle schließen die in Abb. 1 dargestellten 
drei grundsätzlichen Möglichkeiten der Lage 


von Einschlüssen sein und können nebenein- 
ander auftreten. Abb. 2. Lage der Blochwand vor 


bzw. über zwei in einer Ebene 
Die nichtstatistische Verteilung ist für Ein- orientierten Teilchen der Stellung 


schliisse, die durch Ausscheidung entstanden B und © 
sind, grundsätzlich möglich. Es liegt nämlich 
nahe, die Mosaikblock- oder auch die Lamellengrenzen als Orte größerer Keimdichte 
anzusprechen. Die örtlich höhere Keimzahl kann aber eine nichtstatistische Aus- 
scheidung der Teilchen zur Folge haben, die sich längs der gleichmäßig verteilten 
Ebenen ausbilden wird. Die beiden Strukturen sind voneinander zu unterscheiden. 
Durch die Lamellenstruktur wird das Material in Scheiben von sehr gleichmäßiger 
Dicke (0,5—1-10-* cm) unterteilt, während die meist röntgenographisch durch Linien- 
verbreiterung bzw. -Intensität nachgewiesene Mosaikstruktur in den Größen- 
ordnungen von 10~7—10-3 cm eine Unterteilung des Materials in Bléckchen be- 
wirkt. Die Meinungen über das Zustandekommen und das Wesen der beiden 
Strukturen gehen heute noch sehr weit auseinander. Sicher scheint jedoch, daß 
die Grenzflächen beider Phänomene, sei es durch Leer- oder durch Störstellen, 
eine erhöhte Energie besitzen, demzufolge die Keimbildung und damit die Aus- 
scheidung begünstigen. Es wird nur eine Frage der bisher unbekannten energe- 
tischen Verhältnisse sein, wie sich die Zahl der ausgeschiedenen Fremdpartikel 
auf das Korninnere, die Lamellen- und die Mosaikblockgrenzen verteilt. 

Wir bezeichnen im folgenden den relativen Teil der Ausscheidungen, der sich 
an den Lamellengrenzen abscheidet mit c, den an den Mosaikblockgrenzen mit e 
und schließlich den relativen Teil der statistisch im Korninneren ausgeschiedenen 
Partikel mit f. Also, wenn keine weiteren Gründe für nichtstatistische Ausschei- 


dung vorliegen: j= (21) 
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Zunächst soll der Teil der Partikel betrachtet werden, der sich an den Lamellen- 
grenzen abscheidet. Da sich die bekannten Lamellenabstände mit der Konzen- 
tration kaum ändern, wird es nicht schwer fallen, die Teilchenzahl auf einer La- 
mellenfläche zu berechnen. Die Gesamtfläche F, der Lamellengrenzen in einem 
würfelförmigen Kristalliten von der Kantenlänge S beträgt: 

83 


=. (22) 


7 — Lamellendicke (n ~ 10-‘cm). 
S = Kantenlinge des Kristallits. 


Demzufolge wird die Zahl n, der auf der Fläche F abgeschiedenen Teilchen gegeben 
sein durch: 
= F F N 2 
ny, =crnm=cpyF (23) 
ß = Teilchendichte. 
r = Gesamtzahl der Teilchen im Volumen 83. 
c = Zahlenfaktor <1, der den Bruchteil der auf La- 
mellengrenzen abgeschiedenen Teilchen angibt. 
F = Fläche der Blochwand ~ 8%. 


Wenn wir die für 7, wichtige maximale Änderung der Teilchenzahl bei einer Ver- 


schiebung der Wand um 3 berechnen wollen, dann genügt sicher: 


Atmax © (nx 10-4cm; s.o.) (24) 


(Ausdehnung der Lamellengrenzen klein gegen 4 angenommen) : 


Ganz ähnlich liegen die Verhältnisse bei den Mosaikblockgrenzen, nur muß 
man wegen der dreidimensionalen Teilung noch den Faktor } einführen. Also: 


© S* _ (25) 


e = Zahlenfaktor < 1, der den Bruchteil der auf Mosaikblockgrenzen abge- 
schiedenen Nadeln angibt. 
€ = Kantenlänge der Mosaikblöckchen. Bezeichnung sonst wie in Gl. (23). 


Dabei ist allerdings eine weitgehende Ausrichtung der Mosaikbléckchen in den 
einzelnen Ebenen angenommen. Für € sind jetzt nur Werte von über 10-5 cm 
sinnvoll. Denn Schwankungen in geringeren Abständen würden die Wand beim 
Fortschreiten nicht mehr stören, weil sie, ähnlich einem großen Wagenrad auf 
kleinem Kopfsteinpflaster, keine wesentlichen Potentialberge mehr zu überschrei- 
ten hätte. 

Nach Versuchen von Snoek”™) an reinem Eisen mit gelöstem Kohlenstoff bzw. 
Stickstoff wurde öfters diskutiert, daß die Ausscheidung vorzugsweise in den 
Blochwänden stattfinden kann. Dieser Standpunkt ist energetisch durchaus 
plausibel, gehört aber in das Gebiet der magnetischen Nachwirkungen; denn es 
steht wohl außer Zweifel, daß die Anreicherung der Fremdatome in den Bloch- 
wänden nur durch Diffusion und also mit einer Zeitfunktion möglich ist. Wir 


10) J. L. Snoek, Physica 6, 161 (1939). 
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wollen diese weitere Möglichkeit für nichtstatistische Ausscheidungen deshalb 
hier nicht beachten. 

Die statistische Verteilung wird mathematisch etwas einfacher charakterisiert 
als bei Néel, weil in der hier entwickelten Theorie nur die Beeinflussung der 
Koerzitivkraft durch Fremdkörper berücksichtigt wird. 

Gesucht ist also nach Gl. (6) ein vernünftiger Ausdruck für Anmax, die maxi- 
male Änderung der Teilchenzahl in benachbarten „Kästchen“ der Dicke 6. Wir 
teilen den quaderförmig angenommenen Weißschen Bezirk mit den Kanten- 


längen 8, S und Z in > gleich große Kästen (fiir d<6 ,,korrespondierende 


Länge“ = 6). Die Fläche der Wand F wird dabei gleich S®— der entsprechenden 
Seitenfläche des Kristalliten — gesetzt, was bei der praktisch immer vorhandenen 


Lamellierung der Weißschen Bezirke eine gute Näherung darstellt. 3 ist sicher eine 


eroße Zahl, und wir benutzen zur Abschätzung der Wahrscheinlichkeit W (n — n), 
mit der eine Abweichung (n —n) vom Mittelwert der Teilchenzahl in einem 
dieser kleinen Kästen vorkommt, die Gleichung: 
1 _(n-ny? 
W(n—n) = (26) 
V227 
n == mittlere Teilchenzahl pro Kasten. 


Nun wird aber nicht die Abweichung vom Mittelwert gesucht, sondern es 
interessiert uns vielmehr An, die Änderung der Teilchenzahl in benachbarten 
Kästen. Wir fragen also nach der Wahrscheinlichkeit w, mit der in zwei benach- 
barten Kästen die Abweichungen (n— x) bzw. (n—n-+ An) vom Mittelwert 
auftreten. Nach den üblichen Regeln der Wahrscheinlichkeitsrechnung ergibt sich: 


w [(n— n), (mn — n + An)] = w(n— n) wn— + An). (27) 
Also: 
n = n n 2 8 


Aus diesem Ausdruck wird n durch Summation über alle vorkommenden n bei 
vorgegebenen An eliminiert: 


An\* oo 
W (An) = zen n > dé. (29) 
n=0 -00 
Daraus folgt: 

An! 

Win)» 3 (30) 


Nun wird die Wand zum Erreichen von H, nur die Hälfte aller 3 Kästen 


durchlaufen. Die Zahl der zu beachtenden Fälle beträgt also =, = . In diesen rs 


Fällen darf es—so setzen wir fest— höchstens einmal lade ih An >Ammax- 
Dabei bedeutet Anmax jetzt den die Koerzitivkraft bestimmenden Mittelwert 
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von Anmax über mehrere Weißsche Bezirke. Höchstens einmal in 55 2 Fällen 
heißt aber, daß die Summe aller relativen Wahrscheinlichkeiten von Fällen, in 


denen An > Anmax ist, gerade > ergibt. Also: 


An! An? 
= Va _4*max 
— bad dA N 31 
2  2yan | : 
An = Anmax 
oder: 
— Z 
Atmax 2 (32) 
Wegen: 
n=S*6Bf (33) 


ergibt sich schlieBlich: 
= 28 66 fn” 


Diese Statistik erfaßt natürlich nur die nach den Voraussetzungen zu berück- 
sichtigenden Schwankungen der Teilchendichte und nicht etwa auch Schwankungen 
der geometrischen Vorzugsrichtungen anisotroper Teilchen. 


d) Zusammenstellung der einzelnen, für H, wesentlichen Anteile 


Es wäre jetzt entsprechend den Faktoren c,e und f (vgl. 8. 105) zu verteilen, 
und zwar einmal auf den Anteil, der auf den Mosaikblock und Lamellengrenzen aus- 
geschieden wird und zum anderen auf den statistisch im Korninneren eingelagerten 
Anteil. Es ergibt sich nach den jeweils angegebenen Gleichungen und für unmagne- 
tische Teilchen mit rotationselliptischer Form sowie ö > L > dbzw.ö>d> L"): 


| ge [n. 61. (6); (14); 
H, 0,5—1,|3N,—42| (19); (20); (24); (35%)) 
wenn: 
y \3N,—42| 72 
und: 
— [n. Gl (6) ; (13); (19) ; (20) 
yell n. G1(6); (13); (19); (20); 
H,~ (24); (25); (34)], 
2 Z 


11) Auch für längere Teilchen mit L > 6>d gelten die angegebenen Gleichungen 
näherungsweise noch. Bei abgeplatteten Teilchen mit d>6> LZ sind wegen des Auf- 
tretens von N éelschen „Zipfelmützen‘“ größere Abweichungen zu erwarten. i 

Ein re, muß eingeführt werden, da die Teilchenlage A (Abb. 1) für 


— nichts zu H, beiträgt. 
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wenn: 
| 
$>| 


Da « = Av (mit x Volumenanteil der Einschliisse) folgt entsprechend: 


en 
2 ı/övf 
— / —In - 
wenn: 
-4r 
5< I; 
bzw.: 
en | 
H, N 1557, ze‘ ’ (38) 
2 2 
Sy « 
wenn: 
y 3N,—A4n|. 
2 


Sinngemäß wird man in den Gln. (35—38) als letzten Faktor diejenige der drei 
eingeklammerten Größen wählen, die den maximalen Einfluß auf H, hat. 


e) Berücksichtigung des etwaigen Ferromagnetismus von Einschlüssen 


Für ferromagnetische Einschlüsse müssen die bisher gewonnenen Ergeb- 
nisse in sinnvoller Weise abgeändert werden. Dabei wird es von Bedeutung 
sein, ob während des Festhängens der Wand an den Teilchen die Spins im Ein- 
schluß gedreht werden oder nicht. Zur Drehung wird im allgemeinen Energie be- 
nötigt. Diese Energie kann aber nur von dem nach der Spindrehung im Ein- 
schluß vermeidbaren Streufeld geliefert werden. Wir unterscheiden deshalb: 


Y N 2 
b) (5), + (22—N,) (1; [für Lage nach Abb. 1 B] 


und: 


a) ( ö Einschluß > 


y 
b) > 3 Al? + (22 —,) — AN). 
Al, = Differenz der Sättigungsmagnetisierung von‘ Muttermetall und Einschluß. 


+ @a— N) 


(40) 


Gilt die Ungleichung (39), dann wird das Streufeld die Spins im ferromagne- 
tischen Einschluß ,,herumzwingen‘, und wir müssen die geänderten energetischen 
Verhältnisse berücksichtigen. Anderenfalls Gl. (40) wird das Streufeld nicht zur 
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Drehung der Spins im Einschluß genügen, und wir können praktisch unsere in 
den vorherigen Abschnitten für unmagnetische Partikel abgeleiteten Beziehungen 
unverändert benutzen. 

Wenn nur je eine Ungleichung in (39) und (40) erfüllt ist, dann bedeutet das, 
daß z. B. die Teilchen in den Lagen A und B ummagnetisiert werden, während für 
die restlichen Teilchen in Lage C das Streufeld zur Umorientierung der Spins nicht 
ausreicht. Dieser Umstand wird praktisch jedoch nur in extremen Fällen — also 
für Scheiben oder Nadeln auftreten — und man kann auch ihn durch Einführung 
eines geeigneten Faktors im Endergebnis berücksichtigen. 

Wenn (39) gilt, dann müssen weiter die Terme AE|wanaz bzw. AE|wanac 
dahingehend abgeändert werden, daß für die Volumeneinsparung in der Wand 
nicht mehr die Gln. (15) bzw. (18) gelten, sondern daß: 


denn der Energieaufwand gegen Austausch- und Anisotropiekräfte im Einschluß 

beträgt jetzt nicht mehr Null wie bei unmagnetischen Teilchen, sondern in erster 

Näherung (3 — Wandenergiedichte )) Die Streufeldenergieschwan- 
nschlu 

kung bei der Wandverschiebung charakterisieren wir nicht mehr durch Gin. (19) 

und (21), sondern durch: 


AE Streuf. . 


*2 | (42) 


Wir halten als Ergebnis dieses Abschnittes fest, daß bei ferromagnetischen Ein- 
schliissen (unter den gegebenen Voraussetzungen) für H, gilt: 


2 
% 
falls: | 
bzw. 
| 
4,215 7.4(3)15 —eé (44) 
dv f Z 
falls: 


f) Abschließende Betrachtung der Ergebnisse 


Praktisch wird die hier gemachte Annahme, daß die Symmetrieachsen der 
Einschlüsse nur längs der drei magnetischen Vorzugsrichtungen des kubisch an- 
genommenen Mutterkristalls liegen, nicht immer erfüllt sein. Es müßte also ge- 
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gebenenfalls noch eine statistische Verteilung der Symmetrieachsen von Ein- 
schlüssen berücksichtigt werden. Man sieht aber bereits aus den angestellten 
Betrachtungen, daß eine wesentliche Verminderung der Koerzitivkräfte unter 
diesen Umständen erst für d bzw. L > 6 zu erwarten ist. Das gleiche gilt auch 
bei der Verschiebung von 90°-Wänden. In dieser Arbeit sollte zunächst nur 
gezeigt werden, daß Abschätzungen für H, unter Berücksichtigung von Teil- 
chen mit d<6, wie sie in den letzten Jahren erneut von Dijkstra und 
Wert) durchgeführt wurden, nur für exakt kugelige Partikel richtige Ergebnisse 
liefern. Um ein ungefähres Maß für den Fehler früherer Abschätzungen zu be- 
3 N, —*) 2. 

2 4 ö 


für verschiedene Achsenverhältnisse = der Fremdkörper mit rotationselliptischer 


kommen, wurde fiir Eisen in Abb. 3 das Verhiltnis von ( 


Form aufgetragen. Dabei wurde N, nach der fiir Rotationsellipsoide bekannten 


Gleichung: 
| 
= 32[1-15 an (45) 
mit 
d yı-e 
ermittelt. 
| (N,-N;) =(3N,-477) | 
BP | 
| | 
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Abb. 3. Verhältnis von ——*;— "7; zu z=% in Abhängigkeit von T I m Eisen 


13) L. J. Dijkstra u. C. Wert, Physic. Rev. 79, 979 (1950). 
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Man erkennt, daß schon für 4 ew 0,96 bzw. 3 x 1,04 die Streufeldenergie den 
gleichen Einfluß auf H, hat wie die Einsparung an Wandenergie. 

Eine kurze Bemerkung sei noch zu dem Sonderfall gemacht, daß längere Teil- 
chen —‘z. B. Nadeln — gleichmäßig ausgerichtet auftreten. Es würde sich dann 
nämlich wegen der Streufeldenergie eine magnetische Vorzugsrichtung längs der 
Nadelachsen bemerkbar machen, die charakteristische Effekte, wie größere 
Wahrscheinlichkeit von 180°-Wänden und auch feinere Lamellierung zur Folge 
hätte. Nach Beobachtungen von Nesbitt™) kommt die Anisotropie der 
Al—Ni—Co-Magnete nach dem Tempern im Magnetfeld durch die Ausrichtung 
von sehr feindispersen, nadelférmigen stark Co-haltigen Ausscheidungen zu- 
stande, was nach den hier angestellten Betrachtungen energetisch verständ- 
lich ist. Allerdings muß man beachten, daß im Falle des Al—Ni—Co höchst- 
wahrscheinlich die Ummagnetisierung durch Drehprozesse besorgt wird. Ein 
direkter Vergleich der H,-Werte dieses Materials mit den hier gewonnenen Er- 
gebnissen ist also nicht möglich. 


14) R. D. Heidenreich u. E. A. Nesbitt, J. appl. Phys. 28, 352 (1952). 


Jena, Forschungsinstitut für magnetische Werkstoffe. 


(Bei der Redaktion eingegangen am 16. Juni 1952.) 
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Schrittweiser Ubergang von der Wellenoptik zur Strahlenoptik 
in inhomogenen anisotropen absorbierenden Medien. I 


Gleichungen fiir Wellennormale, Brechungsindex 
und Polarisation 


Von Kurt Suchy 
(Mit 2 Abbildungen) 


Inhaltsiibersicht 


Aus den Gleichungen der Wellenoptik werden Gleichungen fiir Wellennormale, 
Brechungsindex und Polarisation unter Bedingungen angegeben, die einem teil- 
weisen Ubergang von der Wellen- zur Strahlenoptik entsprechen. Diese Bedin- 
gungen gelten allgemein fiir inhomogene anisotrope absorbierende Medien. Sie 
beschränken die Krümmung der Inhomogenität des Mediums bzw. des Verlaufs 
der Wellennormalen sowie die Polarisation. 


(10) Einleitung 

In der theoretischen Physik wird der Ubergang von der Wellenoptik zur 
Strahlenoptik mathematisch meist dadurch gewonnen, daß man die Wellen- 
frequenz f nach Unendlich, d. h. die Wellenlänge A nach Null gehen läßt. Dies ist 
zwar mathematisch einwandfrei, jedoch physikalisch nicht exakt realisierbar, da 
Wellenfrequenzen immer endlich bleiben. Physikalisch kann der Übergang also 
nur genähert vor sich gehen, was sich in den mathematischen Formeln durch 
Näherungen bzw. Vernachlässigungen ausdrücken läßt. Dies führt zu denselben 
Formeln wie für {— 00, nur daß diese Formeln jetzt auch mathematisch appro- 
ximativen Charakter haben. Für ihre Gültigkeit ist die Kenntnis der vorausge- 
setzten Vernachlässigungen notwendig. Für manche physikalischen Gegebenheiten 
sind nicht alle diese Vernachlässigungen möglich, und man kann für sie deshalb 
nicht alle Schritte des Übergangs Wellenoptik—Strahlenoptik durchführen, 
sondern muß nach einem bestimmten Schritt aufhören, ehe der vollständige Über- 
gang zur Strahlenoptik erreicht ist. Daher soll der schrittweise Übergang hier be- 
sonders betont werden. 


(11) Aufstellung des Wellengleichungssystems aus den Maxwellschen 


Gleichungen 
Ausgangspunkt der Wellentheorie sind die beiden Maxwellschen Gleichungen 
158 _ y 
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Setzt man 8 = uw: und kann die Zeitabhängigkeit durch D ~ Bm eiet 
eliminieren, so erhält man aus beiden Vektorgleichungen bei örtlich konstantem 


w? 
—VxVx)+Gu-(D+ (11,02) 
Mit Hilfe der Materialgleichungen 
ij=o-G, (11,03) 
worin x und a Dyaden (Tensoren 2. Stufe) seien, und den Bezeichnungen 
2 
k= 5" = = é=el+x+—o}) (11,04) 
0 


ergibt sich aus (11,02) das System linearer homogener partieller Diff.-Gleichungen 
2. Ordnung fiir die E- Koordinaten 


(11,05) 
Mo £o 
Ein analoges System für die 8-Koordinaten 
Vx (Vx 8) + = x(V xB) + BE ze =0 (11,06) 


gilt nur für skalare u, x und o und soll deshalb hier nicht weiter behandelt werden. All- 
gemein kann der B-Vektor aus dem ©-Vektor durch eine der obigen Maxwellschen 
Gl. (11,01) berechnet werden. 


Zur Schreibvereinfachung beim Ana zu den Maßzahlen wird in (11,05) 


für die Dielektrizitätsd yade z = einfach . °- geschrieben. Den Ausgangspunkt 
für die Wellenoptik bildet dann "das „W: 


VX + -E=0. (11,07) 


(12) Einschränkung der Koordinatenkriimmungen im Wellengleichungssystem 


Zunächst soll im Wellengleichungssystem (11,07) der Ausdruck —V x (V x €) 
in AG—VV -€ umgewandelt werden. Dies geht bekanntlich bei Verwendung 
cartesischer Koordinaten ohne weiteres. Mit allgemeinen krummlinigen orthogo- 
nalen Koordinaten miissen jedoch fiir diese Umwandlung bestimmte Vernach- 
lässigungen méglich se‘n. 

Die u, seien die Koordinaten, die e, die Einheitsvektoren (e,-e, = Ö,,), 
und das Linienelement sei 


ey 9, du,. (12,01) 
V-Operator und € Vektor werden hiermit 
3 & @ 3 
V= und E = zeb- (12,02) 


Nach elementaren Differentiationsregeln erhält man mit (12,02) für das Wellen- 


N e,e, ist die Einheitedyade. Die Matrix ihrer Maßzahlen ist 
v 
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gleichungssystem (11,07) 
r3 e eu a 
die Gestalt 


Die Umwandlung von —V x (V x ©) in AE —VV-€ beruht auf der Ver- 
tauschbarkeit der durch die beiden V-Operatoren vorgenommenen partiellen 
Differentationen der &-Koordinaten. Diese Vertauschung ist in (12,03) nur im 
8. (vorletzten) Term möglich. Deshalb müssen die ersten 7 Terme in (12,03) gegen 
die beiden letzten vernachlässigt werden. 


Um diese Vernachlässigung zu deuten, sei eine Hilfsbetrachtung bezüglich der (geo- 


metrisch zunächst unübersichtlichen) Variation der = vorausgeschickt. 
n 


A, gu Au, 


9,Au, . 
(UU, Uy) 
Abb. 1. Allgemeine krummlinige orthogonale Koordinaten 


Man wählt Au, so, daß |g, Au, | = | e,| = 1 ist. Sodann entnimmt man der Figur 
14,9, 4u,,| A, el 
| A, e | 


A 
14, 9, | | | 
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Division durch | g, du, | = 1 ergibt 


| A, In | A, e, 
| 9, Au; | Au, 
2: Kr oe, | 
9, Mu; ou, | 
Durch die Beziehung 
1 1 Qu 
‘va Zu; ce 
J Iu u » 


ist die örtliche Variation der = auf die geometrisch leicht deutbare der e, 
A 
übergeführt. 


Es sollen nun die 3 ersten Terme von (12,03) gegen den 9. (letzten), die nächsten 


4 gegen den 8. (vorletzten) vernachlässigt werden. Da der 9. Term vom Index A 
unabhängig ist, bedeutet die Vernachlässigung des 1. Terms gegen den 9. unter 
Berücksichtigung von (12,04) dasselbe wie die des 2. (für A= u). Da ferner im 
8. Term die Indizes A und u vertauschbar sind, so bedeutet die Vernachlässigung 
des 4. Terms gegen den 8. dasselbe wie die des 7. und die des 5. unter Berück- 
sichtigung von (12,04) dasselbe wie die des 6. Es bleiben also die Vernachlässigungen 
des 2., 3., 6. und 7. Terms zu interpretieren: 


~ 


1 @E, | 
1 @, 9y Ou, 
<= (13,08) 


Faßt man A, V ® „| als Maß für die Wellenlänge im betrachteten Medium auf, 
so bedeutet die Vernachlässigung (12,05): 


Die örtlicke 1. bzw. 2. Variction der Einheitsvektoren e, auf dem Weg einer 
Wellenlänge muß von 2. bzw. 1. Ordnung sehr klein gegen I sein. 

Gröber ausgedrückt: Die östliche Variation der e, auf dem Weg einer Wellen- 
länge darf nur „sehr langsam‘‘ erfolgen, 

Die Vernachlässigung (12.06) bedeutet: 


Die örtliche 1. Variction der Einheitsvektoren muß sehr klein sein gegen 
die relative örtliche 2. Variction der &-Koordincten. 


Ist in (12,05) nur die Länge der Wellen maßgebend, so ist es in (12,06) die 
Veränderung der Elongation, wofür Länge und Amplitude der Wellen sich aus- 
wirken. 

Die Vernachlässigungen (12,05), (12,06) schränken die Krümmung der Koor- 
dinatenlinien ein und sollen deshalb „Krümmungsbedingungen‘‘ genannt werden, 


and 
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_ Mit den Kriimmungsbedingungen (12,05, 12,06) wird das Wellengleichungs- 
system (12,03) zu 
r In ou, ou, 
anders geschrieben 


)) 


& 
9, = const. 
= - (12,07) 
0 
Ix 
Halt man beim Differenzieren die orthogonalen Einheitsvektoren konstant, so gilt 
{-VxVx®}-, = {4E—VV (12,08) 
Ine Ine 


so daß das Wellengleichungssystem (12,07) nunmehr in der Form geschrieben 
werden kann: 


{| 47-77 + =o. (12,09) 
Ix 


Fiir cartesische Koordinaten sind die Kriimmungsbedingungen (12,05) (12,06) 
exakt erfiillt, so daB das Wellengleichungssystem (12,09) hierfiir exakt gilt. 

Hat man jedoch ein krummliniges Koordinatensystem, so bedeuten die Kriim- 
mungsbedingungen (12,05) (12,06) den 1. Schritt von der Wellen- zur Strahlen- 
optik. 

Eine Unterdrückung des Terms VV -€ im Wellengleichungssystem (12,09) würde 


den vollständigen Übergang zur reinen Strahlenoptik bedeuten, wie in einer folgenden 
Arbeit gezeigt werden wird. 


(13) Das Eikonal 

Eine wichtige Hilfsfunktion für die Betrachtung des Überganges Wellen- 
optik—Strahlenoptik ist das Eikonal, und zwar in der Art, wie es von Sommer- 
feld und Runge?) — auf eine Bemerkung von Debye zurückgehend — 
verwendet wurde. Die Anwendung erstreckt sich dort allerdings nur auf eine 
skalare Wellenfunktion, also z.B. eine Koordinate des E-Vektors. Es soll hier 
jedoch die Methode auf die vektorielle Wellenfunktion des gesamten E-Vektors 
erweitert werden. Diese Erweiterung berührt die Polarisationserscheinungen und 
führt nach gewissen Vernachlässigungen, die für die Polarisation bestimmend 
wirken, auf eine skalare Wellenfunktion — eine &-Koordinate — zurück. 

Für den Ortsanteil des Vektors € = führt man die „Ertkonal- 
transformation‘: 


3 
(13,01) 


‘mit komplexen L, durch. Was wird nun aus dem Wellengleichungssystem (12,09) ? 


Aus 
1 a 
— — E, = — LE, 


2) Sommerfeld u. Runge, Ann. Physik 85, 277 (1911). 
3) Zwei gleiche Akzente verbinden einen Differentialoperator und die damit zu dif- 
ferenzierende Funktion. 
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ergeben sich unter Einführung der Dyade 


3 
A= = DL, ey ey (13,02) 
die Beziehungen 
= ike WA)e, und Ee) = tke V -A)e, 
Ix Ix 
Die Auswirkung der zweiten V-Faktoren ergibt sich analog, wobei das Konstant- 
halten der e, und g, zu beachten ist. Setzt man 


4 -A}\., (13,03) 
Gy, 


als ,,Hikonaldyade‘, so wird aus dem Wellengleichungssystem (12,09) das System 
linearer Gleichungen für die &-Koordinaten 


|v—=]-6 (13,04) 


Aus dem Wellengleichungssystem (12,09), einem System linearer homogener 
partieller Diff.-Gleichungen 2. Ordnung für die 3 abhängig Variablen Z,, ist also 
durch die Eikonaltransformation (13,01) ein lineares Gleichungssystem für die E, 
geworden. Die Koeffizienten dieses Systems sind Ausdrücke aus 1. und 2. par- 
tiellen Ableitungen der 3 abhängig Variablen Z,. Zur Bestimmung dieser zu- 
nächst auch unbekannten Koeffizienten steht nur die Bedingung des Verschwin- 
dens der Koeffizientendeterminante zur Verfügung: 


0. (13,05) 


Dies ist eine nichtlineare partielle Diff.-Gl. 2. Ordnung für die 3 abhängig Variablen 
Z,. Mach ihrer Lösung kann man die Maßzahlen der Eikonaldyade M (13,03) 
durch die Maßzahlen der Dielektrizitätsdyade & ausdrücken. Ist dies gelungen, 
so können aus dem ,, Polarisationsgleichungssystem“ (13,04), dessen Koeffizienten 
jetzt nur noch Ausdrücke aus den Maßzahlen von £ sind, die Verhältnisse £,:E,:E, 
bestimmt werden, die die Polarisation charakterisieren. 

Die „Eikonalgleichung‘‘ (13,05) bildet den Ausgangspunkt für alle weiteren 
Rechnungen. Zunächst ist nach dem eben Gesagten ihre Lösung für die Be- 
stimmung der Polarisation aus dem Polarisationsgleichungssystem (13,04) not- 
wendig. Weiterhin ergeben sich aus der Eikonalgleichung (13,05) die Brechungs- 
erscheinungen. 


(14) Einschränkung der Polarisation 
Die Eikonalgleichung (13,05) enthält 3 abhängig Variable L,. Dies erschwert 


ihre Lösung bedeutend, da es keine weitere Gleichung zwischen den drei L, gibt. 


Es soll daher versucht werden, die Eikonalgleichung (13,05) durch Vernachlässi- 
gungen auf eine Diff.-Gleichung für nur eine abhängig Variable zu reduzieren. 
Hierdurch werden dem Zusammenhang der L,, der die Polarisation bestimmt, 
Bedingungen auferlegt. : 

Man wählt ein Z, aus, das nicht auf Grund irgendwelcher Bedingungen überall 
konstant ist. Es sei dies Z,. Setzt man nun 
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so schreibt sich die Eikonaltransformation (13,01) 
En = (€ + hoes + @5) (14,02) 


02 
Hieraus ersieht man, daß L,, und L,, Lage und Form der Polarisationsellipse und 
L deren Größe bestimmen. Stellt man nun die TEN) 
OLys | | Lvs | | 
nt- so enthält die Eikonaldyade M (13,03) und PR Eikonalgleichung (13,05) nur 
noch die eine abhängige Variable L. 
33) Die Polarisationsbedingungen (14,03) bedeuten: 
1. Die 1. bzw. 2. örtliche Variation von Lige und Form der Polarisations- 
a ellipse ist sehr klein gegen die 1. bzw. 2. örtliche Variation der Größe der He 


Polarisationsellipse. 
Gröber ausgedrückt: Die 1. bzw. 2. örtliche Variation der Polarisation ist 
14) sehr klein gegen die 1. bzw. 2. örtliche Variation der Amplitude. 


2. Für die Größen L,, und L,5, die die Polarisation bestimmen, erhält man keine 


er Integrationskonstanten, d.h. keine freien Anfangsbedingungen. Alle Formeln 
so unter den Polarisationsbedingungen (14,03) gelten also nur fiir bestimmte Anfangs- 
E, bedingungen der Polarisation, die durch Lösungen des Polarisationsgleichungs- 
ar systems (14,05) im Wellenanfangspunkt festgelegt sind. . 
ur Man kann sich die Polarisation einer im Anfangspunkt beliebig vorgegebenen 
n- Welle aus 2 Komponenten zusammengesetzt denken, deren eine Komponente — 


als eine Lösung von (14,05) — diese bestimmten Anfangsbedingungen erfüllt. Sie 
soll „Hauptpolarisation‘‘ genannt werden. Mit Hilfe dieser Bezeichnung kann man 


5) Absatz 2 wie folgt formulieren: 
en 2a. Alle Formeln unter den Polarisationsbedingungen (14,03) gelten nur 
3) für die Haurtpolarisation *) . 
n, Die Erfüllung der Polarisationsbedingungen (14,03) bedeutet bei Benutzung 
" eines krummlinigen Koordinatensystems den 2. Schritt‘ von der Wellen- zur 
’ Strahlenoptik, ansonsten den 1. Schritt. 
Unter den Polarisationsbedingungen (14,03) wird VA zu VLI,V-A zu VL 
he und die Eikonaldyade (13,03) zu 
t- My =|VL + (14,04) 
S- 
und das Polarisationsgleichungssystem (13,04) wird zu 

i- a Der weitere Rechnungsgang besteht in der Lösung der Eikonalgleichung 
L. 

+ [VLI—1VL)., — =| = (14,06) 


Ix 


die nur noch eine abhängige Variable Z enthält. 


) 4) Die Ausführungen der Absätze 2. und 2a. verdanke ich Dr. Rawer 


. 
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Die hier mit ,,Eikonalgleichung* bezeichnete Gl. (14,06) ist eine Verallgemeinerung 
der bekannten Eikonalgleichung 


| i ky 
für skalares =. Der Übergang von (14,06) zu (14,07) wird in’ einer folgenden Arbeit 
dargelegt werden. 


Es ist zweckmäßig, in der Eikonalgleichung (14,06) an Stelle des Eikonals I 
als abhäng:ge Veränderliche den „Eikonalgradienten“ 


zu verwenden. Hiermit schreibt sich die Eikonaldyade Mp (14,04) 
In + [nI —In]e. (14,09) 
t ko 
und die Eikonalgleichung (14,06) 
Mp |= [n + In / = (14.10) 


Die physikalische Bedeutung des Eikonalgradienten n ergibt sich aus fol- 
gendem: Nach der Kikonaltransformation (14,02) kann man die Flächen L(r) = 
const. als Wellenfronten auffassen. Demnach ist der Einheitsvektor 

A 

|VL| In 

die Normale auf den Wellenfronten, die ,,Wellennormale“. Den Betrag 


n=|n|=Yn-n (14,12) 


bezeichnet man als „Brechungsindex“. 


(14,11) 


(15) Lösungswege für die Eikonalgleichung durch verschiedene 
Koordinatenwahl. Sommerfeld-Rungesches Brechungsgesetz 


Durch die Eikonalg'eichung (14,10) kann der Eikonalgradient n nicht ein- 
deutig in allen dıei Koordinaten n,, n,, n, bestimmt werden. Man hat vielmehr 
nur die beiden folgenden Möglichkeiten: 


a) Eine Koordinate n, des Eikonalgradienten n wird in Abhängigkeit von’ 


den beiden anderen (n, und n,) ausgedrückt. 


b) Der Betrag n (= Brechungsindex) wird abhängig von der Richtung l 
(= Wellennormale) gefunden. 


Diese beiden Méglichkeiten lassen sich durch passende Wahl des Koordina- 
tensystems ausführen: 

a) In einem „raumfesten‘‘ System richtet man sich nach den äußeren phy- 
sikalischen Bedingungen, d. h. nach der Inhomogenität des Mediums und den 
Anfangsbedingungen für die Welle. 

Liegt z. B. ein ,,geschichtetes‘‘ Medium vor, so wird man eine Schar Koor- 
dinatenflächen u, = const. in die Flächen £ = const. legen und somit & örtlich 
nur von einer Koordinate u, abhängen lasen. Für ebene und Zylinder-Wellen in 
einem geschichteten Medium wird man ferner eine weitere Koordinatenflächen- 
schar z = u, = const. parallel zur Einfallsebene legen, so daß überall !, = 0 wird. 
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b) In einem „wellenfesten“ System (l,, v,, h,) richtet man sich nach dem 
Verlauf der Wellenfortpflanzung und legt eine Koordinatenrichtung I, in die 


ov, 

Wenn jedoch die Anfangsbedingungen der Welle berücksichtigt werden sollen, 
muß in einem vom Wellenverlauf, also von |, unabhängigen raumfesten 
System weitergerechnet werden. 

Im Fall a) beschränken die Krümmungsbedingungen (12,05) (12,06) die Krüm- 
mung der Inhomogenität des Mediums, im Fall b) die Krümmung des Verlaufs 
der Wellenfortpflanzung. 

Zur eindeutigen Bestimmung des Eikonalgradien‘en n muß noch das von 
Sommerfeld und Runge?) aus VYxVL= 0 kergeieitete Brechungsgesetz 


Vxn=Vxnl=0 (15,01) 


herangezogen werden. In Koordinaten: 


Richtung der Wellennorma!en |, so daß überall : O= im wird. 


My 

My 

0 (15,03) 
_ 

eu, 0 . (15,04) 

Die Anfangsbedingungen 

J» = Jvo n= Mh (15,05) 


legen die Integrationskonstanten des Brechungsgesetzes (15,01) (15,02) (15,03) 
(15,04) fest. 

Für den Fall a) liefert das Brechungsgesetz (15,01) (15,02) (15,03) (15,04) 
die fehlende Beziehung zwischen den 3 Koordinaten des Eikonalgradienten n, 
für den Fall b) die Beziehung zwischen Bre.hungsindex n und Wellennormale I. 


(16) Koordinatenanpassung an äußere physikalische Bedingungen 


Für ein geschichte‘es Medium lassen sich unter besonderen Voraussetzungen 
über seine Struktur die Kcordinaten n, und n, des Eikonalgradienten n schon 
aus dem Brechungsgesetz V x n = 0 (15,01) (15,02) (15,03) (15,04) herleiten. 

Wählt man das Koordinatensystem (e,, u,, gy) so, daß 


€=€(u,), (16,01) 
ferner 
Jv = Jv (u,) ’ (16,02) 
so ergibt sich zunächst aus der Eikonalg!eichung (14,10) 
n=n(u) (16,03) 
und dann aus dem Brechungsgesetz (15,03) (15,04) 
J2 Mz = Cy = Jao Ma, d.h. = = Gap la (16,04) 
= Cy = Jao Myo, d.h. Jy nly = Cy = No la. (16,05) 
Hiermit wird 
= Os 16,06 


| 
beit, 
08) 
99) 
LO) 
ol- 
1) 
n- 
| 
“ 
a- | 
y- 
r- 
in 
n- 
d. 


122 Annalen der Physik. 6. Folge. Band 11. 1952 


Legt man die Koordinatenflächen u,= const. im Wellenanfangspunkt parallel 
zur Einfallsebene, so bedeutet dies 


womit sich aus dem Brechungsgesetz (16,05) ergibt 
C,=0=1, (16,08) 
und aus (16,06) der Eikonalgradient 
C 
(u) = en,(u) + eg —— (16,09) 
(uy) 
Mit cartesischen Koordinaten 
(16,10) 
=1 J = 1 92 = 1 
und Zylinderkoordinaten 
Uy = ¢ Ug (16,11) 
lassen sich die Bedingungen (16,02) erfüllen, jedoch nicht mit Kugelkoordinaten 
9,=1 T= 9 = rsind = u, Sin Ug. 


Abb. 2. Kugel-, Zylinder- und cartesische Koordinaten 


Für cartesische Koordinaten bedeutet (16,04) das Snelliussche Brechungsgesetz 


. = Ng loo, (16,13) 
fiir Zylinderkoordinaten das Bourguersche Brechungsgeset: 
qnl, = qo No log - (16,14) 


Mit den Eikonalgradienten n (16,09) unter den Bedingungen (16,01) (16,02) 
(16,07) erhält die Matrix der Eikonaldyade M, (14,09) die Form 


2 
1 0 
92 92 
C 
(Mp) ny ike th ny 0 (16,15) 
0 0 tating 1 
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Die Eikonalgleichung mM ri = 0 (14,10) wird eine gewöhnliche Differen- 
tialgleichung 1. Ordnung für die Koordinate r, des Eikonalgradienten n: 
| Mp— =| = + 4 + + (16,16) 
+ The, My My + My My + + (Gk)? Xo = 
Die Unbekannte Mu läßt sich mit dem Brechungsgesetz g. nl, = C (16,04) in 


(16,17) 
4 


umformen. Das Verhältnis — der beiden Wellennormalenkoordinaten legt wegen 


I, = 0 (16,08) die Wellennormalenrichtung fest. Also kann die Eikonalg!eichung 
in der Gestalt (16,16) als Bestimmungsg'eichung für die Wellennormalenrichtung 


i gelten und soll deshalb ,,Wellennormalengleichung* heißen. 
Die Koeffizienten «,, von (16,16) ergeben sich mit der Bezeichnung 


| Ex 
(16,18) 
Ede 
wie folgt: 
~ Cc? 
C + En 


Die Wellennormalengleichung (16,16) hat im allgemeinen 4 Lösungen n,. Sie 
entsprechen den einfallenden und reflektierten Wellen der beiden durch Doppel- 
brechung entstehenden Wellenzüge. 

Die letzten Schritte zur Strahlenoptik bestehen in der näherungsweisen Lösung 
der Wellennormalengleichung (16,16) durch die WBK-Methode. Dies soll in einer 
folgenden Arbeit durchgeführt werden. 


(17) Koordinatenanpassung an die Wellennormale 


Für diesen in (15) unter b) erwähnten Fall ist die Richtung von n maß- 
gebend, nach (14,11) also die Wellennormale |. 


Wegen der zentralen Bedeutung von [ bei der Koordinatenwahl in diesem Absatz 
sei zunächst etwas näher auf die physikalische Bedeutung von | eingegangen. 
Da L im allgemeinen eine komplexe Größe ist, sind durch (14, 11) (14,12) auch die 
Größen [ und n im allgemeinen komplex. Ein komplexer Vektor ® in einem Koordinaten- 
system mit reellen Einheitsvektoren e, hat komplexe Koordinaten V,. Es ist dann 


Re® = Le, Re V,, Im 8 = Le, Im y,. 


) 
)) 
) 

| 

| 
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In Anlehnung an Försterling, S. 108)5), ergeben sich für | und m folgende Re- 
lationen: 


Wegen ist 
Rel-ImI=0. (17,01) 
Wegen 
L(t) =fdr-VL =f dt-nt=f de-(RenRel—Imn Im!) + ¢f 


und 


eb —koImZ ei 


sind fiir reelles dr die reellen Vektoren 
RenRel—ImnIm! bzw. A= (17,02) 
den Phasen- bzw. Amplitudentrajektorien parallel. Wegen (17,01) ist 


=€ 


P-A=RenImn. (17,03) 
Ferner sei bemerkt: Nach (14,11) (14,12) wird 
Im! = 0, falls Imn=0. - (17,04) 
Daher verschwinden nach (17,02) die Amplitudentrajektorien: 
A=0, falls Inn = 0. (17,05) 


Wenn also Im n = 0 ist, so bleibt die Amplitude konstant und es findet keine Absorption 
statt. 

Die Anpassung des Koordinatensystems an die Wellennormale | geschieht 
folgendermaßen: Man wählt (das gegebenenfalls komplexe) { als 1. Einheitsvektor |,, 
die beiden anderen Einheitsvektoren [,, [, beliebig senkrecht dazu, eventuell I, so 
daß es reell bleiben kann. In diesem System von ,,Wellennormalenkoordinaten* 
(I, v», hy) wird dann — siehe (12,02) — 


I, aL (v,) 


h, (v1 vg, dv, (v1, V2, Us), (17,06) 
und 
[72 3-2 Li, 1 


Die Matrix (Mp) der Eikonaldyade Mp (14,09) wird mit (17,06) (17,07) 


0 0 0 
1 1m 1 
(Mp)=|0 0 = (n+ (17,08) 


Die Form (17,08) der Eikonaldyade Mp für Wellennormalenkoordinaten ergibt 
sich aus der Fom (16,15) durch 


(17,09) 


8) pases Lehrbuch der Optik, 1. Aufl., Leipzig 1928. 
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Die Koeffizienten «,, der Eikonalgleichung in Wellennormalenkoordinaten 
erhält man mit (17,09) entsprechend (16,19) zu 


C= = 0} = 2a, 
=0}=0 n{C=0}=0 

= 0) = (17,10) 

Om { = 0} = — =4 = 0} = a {C = 0} = ay. 


. Hiermit läßt sich die Eikonalgleichung für Wellennormalenkoordinaten wie folgt 


schreiben: 


i 18 
Man ersieht hieraus sofort, daß es bequem ist, zuerst die algebraische ,,Doppel- 
brechungsgleichung‘‘ 


|w |= + +a = 0 (17,12) 
zu lésen und dann die nichtlineare baa Diff.-Gl. 
1 on 
die „Brechungsindexgleichung‘‘. Ein Spezialfall von (17,13) erstmalig bei 
Försterling®). 
Die beiden Lösungen Ri 
Ag 


“der Doppelbrechungsgleichung (17,12) entsprechen den beiden Hauptpolari- 


sationen (14) der durch Doppelbrechung erzeugten Wellenzüge, wobei jedoch 
nicht im ganzen Wellenverlauf w, dem einen, w_ dem anderen Wellenzug zu 
entsprechen braucht, worauf Poeverlein’) im Spezialfall der Ionosphäre hin- 
wies. Man hat sich vielmehr die Doppelbrechungsgleichung (17,12) als Defini- 
tionsgleichung einer zweiwertigen algebraischen ‚Doppelbrechungsfunktion“ 
w= w(r) vorzustellen, wobei die Zuordnung der w, zu den beiden Wellen- 
zügen auf dem Verzweigungsschnitt der zweiblättrigen Riemannschen r-,,Flache“ 
wechselt. (In einem geschichteten Medium mit &=E(u,) wird w=w(u,) und 
man hat eine zweiblättrige w,-Flache.) 

Nach Hines®) kann man die Diskriminante der Doppelbrechungsgleichung 
(17,12) in der Form schreiben: 


ote as — ey 


*) Försterling, Ann. Phys. 11, 1 (1931). 
7) Poeverlein, Z. angew. Phys. 1, 517 (1949). 
8) Hines, J. Geoph. Res. 56, 63 (1951). 
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Hiermit und mit (17,10) wird die Doppelbrechungsfunktion (17,14) zu 


In einem nichtabsorbierenden Medium ist © hermitisch, d. h.2,,=®%,. Hiermit 


wird gemäß (16,18) = = er und daher D = 0 (reell). Außerdem wird |?] reell und nach 
(17,10) werden hiernach alle a, reell. Also ergeben sich für ein nichtabsorbierendes Me- 
dium zwei reelle Elemente w, der Doppelbrechungsfunktion w 


In einem isotropen Medium ist é ein Skalar. Die hierfür geltenden Formeln gewinnt man 
aus denen für dyadisches # durch den Übergang ¢->@J. Hiermit wird iu, = €6,, und 


= = e706, . ö,,, daher D= 0 und w= - . Das Verschwinden der Diskriminante D 


und das Einwertigwerden der Doppelbrechungsfunktion w für skalares & bedeutet das 
Verschwinden der Doppelbrechung in isotropen Medien. 

Nach.der Brechungsindexgleichung (17,13) gehören zu jedem Doppelbrechungs- 
Funktionselement w. zwei Lösungen n, die zu einer Koordinatenrichtung I, ge- 
hören und daher 2 in entgegengesetztem Richtungssinne laufenden Wellen mit 
I=-[I, entsprechen. Einfallender und reflektierter Welle entsprechen sie nur 
in denjenigen Fällen, in denen die Normale der einfallenden und die der reflek- 


tierten Welle genau entgegengesetzt verlaufen. Dies geschieht nur in Spezial- . 


fällen und im allgemeinen sind die beiden Lösungen n der Brechungsindexgleichung 
(17,13) Ausdruck eines Reziprozitätsgesetzes, wonach in einem System von Wellen- 
normalenkoordinaten (f,, v,, A,) in beiden Richtungen + I, ein Wellenverlauf 
möglich ist. 

Die Brechungsindexgleichung (17,13) kann man unter Berücksichtigung von 
(17,09) formal als Spezialfall der Wellennormalengleichung (16,16) auffassen mit 


= = App = = Ay = = 0, = X= 1, = — 


Ihre Lösung durch die WBK-Näherung ist daher ein Spezialfall der WBK~ 
Näherungslösung von (16,16), wie in einer folgenden Arbeit gezeigt werden wird. 

Die Methode, mit Wellennormalenkoordinaten zu rechnen, hat den Nachteil 
der Unanschaulichkeit, wenn die Wellennormale komplex ist, da man dann in 
einem Koordinatensystem mit (mindestens zwei) komplexen Einheitsvektoren 


zu rechnen hat. 
i nichtabsorbierenden Medien mit reellem und positivem w, ferner mit Unter- 


drückung des Terms + ie = in der Brechungsindexgleichung (17,13), also im Fall 
reiner Strahlenoptik, erhält man zwar einen reellen Brechungsindex n und damit laut 
(17,04) auch eine reelle Wellennormale [, jedoch nicht in absorbierenden Medien. Darum 
ist auch das in (18) erwähnte graphische Verfahren von Poeverlein ®) exakt nur möglich 
im Fall der reinen Strahlenoptik in nichtabsorbierenden Medien. 


Ein weiterer Nachteil ist im allgemeinen die Unkenntnis von h, in der Brechungs- 
indexgleichung (17,13). h, hängt von der geometrischen Gestalt des Koordinaten- 
systems (l,, v,, ,) ab, für dessen Kenntnis man aber erst den Verlauf der Wellen- 


8) Poeverlein, S.-B. Bayer. Akad. Wiss. 1948, 175. 
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normalen [ = I, in einem raumfesten System bestimmen muß. Dies kann man aber 
nur mit der Methode a) aus (15). Nur in 2 Fällen hat man diesen Umweg nicht 
nötig: 

1. Man weiß von vornherein, daß die Wellennormale immer konstant bleibt. 
(Z. B., wenn eine Welle in ein geschichtetes Medium mit € = E (u,) genau senkrecht 
zur Schichtung einfällt, also mit [= e,.) Das dem Wellennormalenverlauf ange- 
paßte Koordinatensystem (I,, v,, h,) ist dann dasselbe wie das dem Medium ange- 
paßte System (e,, %,, g,) und daher wird h, = g,. Für geschichtete Medien 
unter den Bedingungen g, = g, (u,) (16,02) wird A, = h, (v) und nun hängt n nur 
noch von v, = &, ab. Hierdurch wird die Brechungsindexgleichung (17,13) zu einer 
gewöhnlichen riccatischen Diff.-Gleichung. 


2. Man kann rein strahlenoptisch rechnen, d.h. das Glied es 2 in in der 


tke © 
Brechungsindexgleichung (17,13) unterdriicken. 


Zur Aufstellung der Doppelbrechungsgleichung (17,12) sei noch eine Bemerkung 
angefiigt. 
Gemäß (17,10) benötigt man ‘fiir die Koeffizienten a, die Kenntnis der Maßzahlen — 
der Dielektrizitätsdyade 


(11,04) 


Nun erhält man aus der Dispersionstheorie zunächst die Dyaden x"! bzw. o! und müßte 
zur Bestimmung der Maßzahlen von £ erst mühsam die Reziproken von x! bzw. r! 
bilden. Man kann dies wie folgt az 
Setzt man 

n=x*+ (17,18) 

so wird aus dem Polarisationsgleichungssystem (13,04) 
[ae =] -€=0, (17,19) 
€ ~ 

durch linksseitige skalare Multiplikation mit [M — IT! 


Ip -[Z 1-4 .€=0. (17,20) 
Aus (17,08) und (17,13) entnimmt man fiir Wellennormalenkoordinaten 
M, = w {I —1, 4), (17,21) 
womit man nach kurzer Rechnung erhält: 
Hiermit wird das Polarisationsgleichungssystem (17,20) zu 
1 
[- Lh — [2] | (17,28) 
0 


Dieses Gleichungssystem kann nun an Stelle des Polarisationsgleichungssystems 
€ = 
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benutzt werden, wenn die MaBzahlen von n7! leichter zu erhalten sind als die von €. An 
Stelle der Doppelbrechungsgleichung 


| e | 
| — =| -0 (17,12) 
{ fo | | fo | 
tritt dann 


w— 1)? 


4 <w 1)-1 [I L h] | = 0, (17,25) 


was ausmultipliziert wieder die Doppelbrechungsgleichung (17,12) mit denselben Koeffi- 
zienten a, ergibt, die diesmal nicht aus den Maßzahlen von = = &,/ +.n, sondern aus 
denen von n7! berechnt werden. 


(18) Bestimmung der Wellennormalen durch das Brechungsgesetz von 
Sommerfeld-Runge 


Hat man die Lösungen n = n(v,) der Brechungsindexgleichung (17,13) für 
alle Richtungen der Wellennormalen |, an einem Ort r bestimmt und will nun mit 
Hilfe der Wellenanfangsbedingungen die gesuchte Wellennormalenrichtung | und 
das hierzugehörige n am Ort r finden, so muß man zunächst in ein raumfestes 
Koordinatensystem (e,, u,, g,) transformieren, in dem die Wellenanfangsbedin- 
gungen (15,05) und der Verlauf von | festgelegt werden können. 


Es sei 


u, 
I, V2 = "I Us], (18,01) 
es Us 


wobei 7'-! die zu 7 kontragrediente Matrix ist. Die Elemente der Transforma- 
tionsmatrix 7' sind die Koordinaten der |, im System (e,, u,, g,). Also lassen sich 
nach den Transformationsformeln (18,01) die v, durch die I, (w,) und die uw, aus- 
driicken. Es ist also 


n (vy) = n (L,(u,y), Uy). (18,02) 


Zur Bestimmung von | mit Hilfe von n ({,, u,) benützt man das Sommer- 
feld-Rungesche Brechung:gesetz 


Vxnt=0 (15,01) 


Für geschichtete Medien unter den Bedingungen (16,01) (16,02) (16,07) ergibt 
sich aus (16,04) (16,08) 


Ll 
(18,03) 


Kann man |Im n|< |Ren| annehmen, so erhält man nach (14,11) (14,12) auch an- 
nähernd rein reelle | und g,. (Wie in (17) ausgeführt, werden im nichtabsorbierenden 
Medium mit positivem w für reine Strahlenoptik n und I reell und somit auch die g, reell.) 


In diesem Fall kann man aus (18,02) (18,03) graphisch sehr übersichtlich für jeden Ort 
u, denBrechungsindex n und die Wellennormale { bestimmen [Poeverlein °) !0)]. 


10) Poeverlein, Z. angew. Phys. 2, 152 (1950). 
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(19) Die Polarisation 
Zur Untersuchung der Polarisation muß natürlich das Koordinatensystem der 
Wellennormale angepaßt werden. Aus dem Polarisationsgleichungssystem 
M = e = 
(17,24) 


für die ©-Koordinaten erhält man unter Benutzung von 2. und 3. Zeile 
-ı2 


€ € € € € 
E,:E,:E, = =) : (= w— =) 


4 a): +): 


(Die 1. Zeile ergibt Ausdrücke mit w*, was die Formeln unnötig komplizieren 
würde.) Die Polarisation wird durch das Verhältnis 


charakterisert. Mit a w. aus (17,16) erhält man 


Da im Polarisationsgleichungssystem (17,24) die Polarisationsbedingungen (14,03) 
vorausgesetzt sind, handelt es sich in (19,03) um die Hauptpolarisationen der beiden 
durch Doppelbrechung entstandenen Wellenzüge. Aus (19,03) entnimmt man sofort 


>) ‘ 3) 
2). = (19,04) 
In nichtabsorbierenden Medien wird für hermitisches = 
_ 


d.h. die beiden Hauptpolarisationsellipsen sind senkrecht zueinander gekreuzt. Für 
symmetrisches reelles 2 wird (Försterling S. 180)°) 


2) (2) || E. E 

d.h. die beiden Hauptpolarisationsellipsen sind senkrecht zueinander gekreuzt und wer- 
den im gleichen Sinne durchlaufen. 


Rechnet man bequemer mit den Maßzahlen von nt anstatt von = = eI +n, so 
bestimmt man zunächst aus (17,23) die Koordinaten des Vektors 


i ou 
und erhält hieraus die Koordinaten von © wie folgt: Nach (17,19) ist 
= [M —1]- (19,08) 
Mit 
M, = w (I (17,21) 
wird aus (19,08) 
G6 (19,09) 


Aun. Physik. 6. Folge, Bd. 11 9 
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und somit ° 
= —G, 
& B, = (19,10) 
& Es = <w — 13710, 
Die Größe <w — 137! bestimmt man aus der Doppelbrechuitgsgleichung (17,25). 
Zum Schluß sei noch erwähnt, daß man für die Polarisation die Brechungs- 
indexgleichung 


gar nicht zu lésen braucht, sondern nur die Doppelbrechungsfunktion w zu er- 


mitteln hat. Der Übergang zur reinen Strahlenoptik (Unterdrückung des Terms 


: 2 =) hat also mathematisch auf die Polarisation keinen Einfluß, was 
0 1 1 

physikalisch dadurch plausibel erscheint, daß die Polarisation als reiner Wellen- 
begriff in der reinen Strahlenoptik sinnlos wird und nicht etwa einem Grenz- 


wert zustrebt. 


Die Anregung zu vorliegender Arbeit verdanke ich Herrn Prof. Seeliger. 
Ihm und Herrn Dr. Rawer danke ich für wertvolle Hinweise und Diskussionen. 


Freiburg/Br., Tivolistr. 36. 


(Bei der Redaktion eingegangen am 25. Juni 1952.) 
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Die Eigenladungs-Verbreiterung 
eines fokussierten Elektronenbündels 


Von Hans Griimm 


(Mit 5 Abbildungen) 


Inhaltsiibersicht 


Die bisherigen Arbeiten iiber die Verbreiterung, welche ein Elektronenbiindel 
durch die Wirkung der eigenen Raumladungerfährt, sind zumeist von geometrisch- 
optischen Gesichtspunkten ausgegangen. In den Vereinigungsgebilden von 
Elektronenbiindeln treten jedoch erhebliche Abweichungen der Intensitätsver- 
teilung von geometrisch-optischen Verhältnissen auf, die der Größenordnung nach 
die Eigenladungs-Effekte weit übersteigen. Es wird daher hier der Versuch unter- 
nommen, von einer durch die Wellenmechanik bestimmten Bündelstruktur 
auszugehen. Das Ziel besteht darin, einfache Formeln für die oberen Grenzen der 
Eigenladungseffekte (Bündelverbreiterung, Verschiebung der Bildebene, Ver- 
schiebung der Fraunhofer-Ebene usw.) aufzustellen. 


1. Einleitung 


Bei der Durchrechnung elektronenoptischer Systeme sieht man meist davon 
ab, daß die Elektronen (bzw. Ionen) als geladene Teilchen aufeinander wirken 
und im allgemeinen zu einer Verbreiterung der Elektronenbündel führen. Diese 
Verbreiterung wurde in mehreren Arbeiten!)?)®)*) unter vereinfachenden An- 
nahmen untersucht. Es ergab sich dabei, daß man angesichts der außerordentlich 


geringen Strahlstromstärken in den meisten elektronenoptischen Geräten (beim - 


Mikroskop z. B. Mikroampere) in der Tat oft von der Eigenladung des Bündels 
absehen kann. In den zitierten Arbeiten wurden ein paralleles Bündel, ein gegen 
einen Punkt konvergierendes Bündel und schließlich Bündel mit ,,geometrisch- 
optischer Struktur‘ untersucht. Der Raum wurde in allen Fällen als feldfrei, die 
Stromdichte in jedem Strahlquerschnitt als konstant angenommen. Diese -An- 
nahmen haben den Vorteil, daß die Lösungen auf tabuliert vorliegende Funktionen 
zurückgeführt werden können. Arbeitet man konsequent geometrisch-optisch, 
dann treten aber in den Strahlvereinigungs-Gebilden (Brennpunkte, Kaustik usw.) 


1) E. E. Watson, Philos. Mag. 3, 894 (1927). 

2) M. Knoll u. E. Ruska, Ann. Physik (5) 12, 604 (1932). 

3) B. v. Borries u. J. Dosse, Arch. El. Techn. 32, 221 (1938). 

4) G. Wendt, Ann. Physik 5—6, 256 (1948). Über eine experimentelle Ermittlung 
von Raumladungswirkungen im Gummimembran-Modell vgl. B. W. Bobykin, B.M. 
Kelman und D: L. Kaminski, J. techn. fis., Moskau, XXII, 736 (1952). 
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ohne Berücksichtigung der Raumladung unendlich hohe Stromdichten auf und 
stellen die Berechtigung dieser Betrachtungen in Frage‘). 

Es liegt daher nahe, die Eigenladungs-Effekte vom Standpunkt der Wellen- 
mechanik aus zu untersuchen, die in allen Bündelquerschnitten zu physikalisch 
sinnvollen Stromdichteverteilungen führt. Außerdem wollen wir uns von der 
einschränkenden Bedingung des feldfreien Raumes befreien und Verhältnisse 
betrachten, die den in der Elektronenmikroskopie vorliegenden möglichst ähnlich 
sind. Die Möglichkeit zu dieser Untersuchung wurde in letzter Zeit dadurch ge- 
geben, daß W. Glaser®) und seine Mitarbeiter die wellenmechanischen Strom- 
dichteverteilungen in verschiedenen typischen fokussierten Elektronenbündeln 
berechnet haben. 


2. Verteilung der Stromdichte (Raumladung) in einem typischen fokussierten 
Elektronenbündel 


Um die durch die Eigenladung bewirkten Effekte in einem fokussierten Elek- 
tronenbündel kennenzulernen, betrachten wir eine einfache, aber für viele elek- 
tronenoptische Geräte (vor allem für das Mikroskop) typische Anordnung. Ein 
zylindrisches Elektronenbündel, beschleunigt durch die Spannung U, trete in der 
Objektebene z = z axial in ein rotationssymmetrisches Abbildungssystem ein 
und werde durch dieses modifiziert. Wir beschränken uns im folgenden auf den 
paraxialen Bereich, beachten aber, daß sich unsere Überlegungen sinngemäß 
auch auf die Seidelsche Dioptrik übertragen lassen (allerdings werden sich die 
Eigenladungseffekte dabei nur um vernachlässigbare Größen ändern). Auch 
andere ,,Beleuchtungsverhiltnisse“ in der Objektebene lassen sich leicht berück- 
sichtigen, ohne wesentlich Neues zu liefern. Schließlich sei noch bemerkt, daß 
wir hier nur die elektrostatische Abstoßung zwischen den Elektronen im Bündel 
berücksichtigen und nicht relativistisch rechnen’). 

Vom Standpunkt der geometrischen Elektronenoptik aus sind die Paraxial- 
bahnen ohne Berücksichtigung der Eigenladung des Bündels bestimmt durch die 
“-Differentialgleichung 

"+ + + 5; Bi) r = 0. (1) 
Dabei bedeutet ®(z) den Verlauf des Linsenpotentials längs der Feldachse, während 
B,(z) die Axialkomponente der magnetischen Feldstärke, ebenfalls längs der 
Feldachse bedeutet. Mit zwei partikulären Integralen von (1), s(z) und ¢(z), welche 


5) In der Arbeit von G. Wendt ist dementsprechend die konsequent geometrisch- 
optische Betrachtungsweise verlassen worden, um zu physikalisch sinnvollen Formeln 
zu gelangen. G. Wendt legt einen geometrisch-optischen Strahlengang zugrunde, nimmt 
aber auch in den kritischen Bündelquerschnitten (wo Einhüllende auftreten) eine endliche, 
über den ganzen Strahlquerschnitt konstante Verteilung der Stromdichte an. Wir werden 
sehen, daß dieses Vorgehen vom Standpunkt der Wellenmechanik aus als Näherungs- 
verfahren völlig gerechtfertigt ist. 

6) W. Glaser, Grundlagen der Elektronenoptik, Springer-Wien (im Druck). Der 
Verf. dankt Herrn Prof. Glaser an dieser Stelle für den Einblick in das noch unver- 
öffentlichte Buch. 

?) Die relativistische Korrektur kann ohne weiteres berücksichtigt werden und führt 
lediglich zu komplizierteren Ausdrücken. In den in der Praxis auftretenden Fällen ist 
sie jedoch vernachlässigbar klein. Im übrigen führt sie zu kleineren Effekten, als die 
klassische Rechnung. Da wir uns vor allem für die oberen Schranken der Raumladungs- 
effekte interessieren, sind wir um so mehr berechtigt, klassisch zu rechnen. 
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den Bedingungen 
8(%) = 1; 4%) = 0 
8'(%) = 0; t'(%) = 1 
geniigen, lassen sich die Paraxialbahnen in der Gestalt 
r(z) = 9° 8(2) + (3) 
schreiben. Das in z = z parallel zur Achse einfallende Bündel vom Radius R 
ist also durch die Bahn 


(2) 


r(z) = R- s(e) (4) 
gegeben. Das fokussierte Biindel stellt ein doppelkegelartiges Gebilde dar. Die 
Abb. 1 zeigt den Meridianschnitt eines solchen- Bündels, erzeugt durch ein rein 
magnetisches Abbildungssystem vom Glaserschen Glockentypus (Linsenstärke 
k* — 1) bei einer Beschleunigungsspannung U = 100 kV und einer Vergrößerung 
V ——10. Der Eintrittsradius R des Bündels wurde mit 10-5 d angenommen (die 
Feldhalbwertsbreite d dient als Maßeinheit). In die Abb. 1 wurde auch die (in allen 
achsensenkrechten Biindelquerschnitten konstante) Stromdichte J, (z) eingetragen. 
Wie man sieht, versagt die geometrisch-optische Lösung physikalisch an der 
durch s(zp) = 0 definierten Stelle (Ort der „„Fraunhofer-Ebene‘“). 


Abb. 1. Meridianschnitt durch ein fokussiertes Elektronenbündel nach der klassischen 
Mechanik. Die Stromdichte ist in allen Querschnitten konstant. Der Deutlichkeit hal- 
ber wurde die Ordinate 5 - 10?-fach gestreckt 


Bei der wellenmechanischen Betrachtung wird die Stromdichte an keiner 
Stelle des Bündels unendlich. W. Glaser und P. Schiske®) haben eine den 
Zwecken der Elektronenoptik angepaßte Lösung der Schrödinger-Gleichung 
angegeben, die es gestattet, aus der in der Dingebene z = z vorgegebenen Wellen- 
funktion p(z) die Stromdichteverteilungen in den einzelnen Einstellebenen z = 
const. zu berechnen. Spezialisiert man auf das paraxiale Gebiet (paraxiale Schrö- 


‚dinger-Gleichung) und läßt man die Bündelladung außer acht, so ergibt sich in 


unserem Fall nach den erwähnten Autoren die Stromdichteverteilung: 


= 28 (Ui (a, br) + Ui (a, br)]. (5) 
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Dabei bedeutet J, (z) = J, die (konstante) Stromdichte bei z = z. Die Funktionen 
a(z) und b(z) sind gegeben durch 
2 
mit p= \V2emU. (6) 
Die Funktionen s(z) und ¢(z) sind identisch mit den durch (2) festgelegten parti- 
kulären Lösungen des „zugeordneten geometrisch-optischen Problems“ (1). 
Unter U, und U, sind die Lommelschen ,,U-Funktionen“ zu verstehen: 


U,= 


Mit J wurden hier die gewöhnlichen Besselfunktionen bezeichnet. In Abb. 2 
ist die wellenmechanische Stromdichteverteilung zum Vergleich mit der geo- 
metrisch-optischen (Abb. 1) dargestellt (die Auswertung erfolgte durch E. Giitter *) 
Außer den Dichteverteilungen in den verschiedenen charakteristischen Einstell- 
ebenen sind die Bahnen R - s(z) eingezeichnet, welche die geometrisch-optische 
Begrenzung des Bündels darstellen. Man bemerkt, daß die Stromdichte in der 
Fraunhofer-Ebene z = 2p endlich bleibt. Es ist daher durchaus möglich, wie 
das G. Wendt?) getan hat, auch in den Vereinigungsgebilden von Elektronen- 
bündeln zur näherungsweisen Berechnung eine über das ganze Bündel konstante, 
endliche Dichteverteilung anzumehmen. 


| VR 


7, 


Abb. 2. Meridianschnitt durch dasselbe Bündel wie in Abb. 1. Die Stromdichtever- 
teilungen sind hier auf Grund der Wellenmechanik bestimmt 


Eine eingehende Diskussion von (5) ergibt in der Objektebene z = z für das 
(geometrisch-optische) „Lichtgebiet“ 


J, (20 r) oii Jo (r s Rk), (8) 


®) Vgl. Anm. 5). Die geometrisch-optische Bahn R - s(z) entspricht etwa der „inneren 
Bahn“ bei Wendt. 
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während im ‚Schatten‘, keine Elektronen auftreten. In der Gaußschen Bild- 
ebene z = 2,, die durch ¢(z,) = 0 definiert ist, herrscht im Lichtgebiet die Strom- 
dichte 

J,(2,r) = J/V? (rs RV), (9) 
während sie im Schatten ebenfalls verschwindet. In der Fraunhofer-Ebene 
2 — zp schließlich zeigt sich das Airysche Beugungsscheibchen: 


=Io( by RY 


brr 
Die hier geschilderte Intensitätsverteilung im Bündel soll uns bei den folgenden 
Betrachtungen als konkretes Beispiel dienen, ohne daß wir uns an die speziellen 
Verhältnisse gebunden fühlen. 


(10) 


3. Verschiebung der Bildebene und Änderung der Vergrößerung 


Der grundlegende Gedanke unserer Untersuchungen besteht in folgendem. 
In jeder achsensenkrechten Ebene hinter der Objektebene herrscht eine durch die 
paraxiale Schrödinger-Gleichung bestimmte Stromdichteverteilung J,(z, r). 
Das Bündel kann als stationäre Raumladungswolke aufgefaßt werden, wobei die 
örtliche Raumladungsdichte mit der Stromdichte nach 


m 


zusammenhängt. Diese Raumladungswolke kann gewissermaßen als Raum- 
ladungs-Zerstreuungslinse aufgefaßt werden, die auf den Abbildungsgang ein- 
wirkt. Diese „Eigenladungslinse“ wollen wir zu den „äußeren“ Feldern D(z) 
und B,(z) dazuschlagen und das wellenmechanische Problem unter diesen ver- 
änderten Bedingungen erneut durchrechnen. Es ergibt sich auf diese Weise ein 
Iterationsverfahren, das man aber wegen der von vornherein als sehr klein voraus- 
gesetzten Wirkung der Eigenladung mit dem ersten Schritt abbrechen kann 
(eine Verbesserung der so gewonnenen Näherung, insbesondere bei hohen Strahl- 
strömen, ist prinzipiell möglich). 

Wie W. Glaser®) gezeigt hat, ist das wellenmechanische Problem (im par- 
axialen Bereich) im Prinzip gelöst, wenn die zugeordnete geometrisch-optische 
Gl. (1) gelöst ist. Wir wollen für das Folgende annehmen, daß wir über die beiden 
partikulären Integrale s(z) und ¢(z) ohne Berücksichtigung der Raumladung ver- 
fügen. Unsere Aufgabe reduziert sich dann darauf, die Gl. (1) für den Fall zu lösen, 
daß außer dem elektrisch-magnetischen Abbildungsfeld noch eine (z. B. durch (5) 
gegebene) wellenmechanisch ermittelte Raumladungswolke vorhanden ist. 

Bei Vorhandensein einer von r unabhängigen Raumladungsverteilung o(z) 
lautet die paraxiale Bahngleichung bekanntlich: 


(2,7) = (11) 


[73 (77 e 2 1 
t+ + ee |r = 0. (11) 


Dabei bedeutet ¢ die absolute Induktionskonstante (praktisches Maßsystem). 
In unserem Fall haben wir aber auch eine r-Abhängigkeit der Raumladung zu 
berücksichtigen, die ja innerhalb des paraxialen Gebietes stark schwankt. Um 
die Verallgemeinerung von (11) zu finden, betrachten wir eine achsensenkrechte 
zylindrische Scheibe der Breite Az und vom Durchmesser r. Der Fluß durch die 
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Oberfläche dieser Scheibe muß gleich sein der e,-fachen Gesamtladung im Innern 
der Scheibe. Die Radialkomponente der Feldstärke ergibt sich auf diese Weise zu 


z 


(r, 2), (12) 


wobei F(r,z) die Gesamtladung im Innern der Scheibe, pro Längeneinheit der 
z-Achse bedeutet: 


F (r,z) = 22 fot u)u du. (13) 


Mit Hilfe des auf die paraxiale Bewegung spezialisierten Energiesatzes eliminieren 
wir nun aus der Bewegungsgleichung 


7=——E, | (14) 


die Zeit®) und erhalten, wenn wir die Beziehung EZ, = —®’ (z) beachten, an Stelle 
von (1) die ee 
Diese Gleichung bestimmt die Bahnen 7 (z) unter Berücksichtigung der Biindel- 
ladung. Die Gl. (15) hat auf der Achse r = 0 keine Singularität, wie es auf den 
ersten Anblick scheinen möchte. Wie man durch Einsetzen einer konstanten 
Raumladung statt o(z,%) erkennt, ist das Integral dann proportional zu r?, was 
sich gegen das r? im Nenner hebt 
Voraussetzungsgemäß soll die Raumladungswirkung gering sein. Wir können 
daher das Raumladungsglied als Störglied auf die rechte Seite schaffen und die 
Integration bis zur ungestörten Bahn r ARE: 


S(r,2)= -F(r, z). (16) 


47 & r 
Die Lösung der homogenen Gl. (1) sei in der Gestalt (3) vorgegeben, dann ergibt 

sich als angenäherte Lösung von (15): 

(17) 


r(z)=s n— : dz 1 
r 


Dabei haben wir bei der Berechnung der ERS ER von der be- 
kannten Beziehung 


(st' — 8’ t) = const. (18) 
(Helmholtz-Lagrangescher Satz) Gebrauch gemacht. 


Ein erster, durch den Einfluß der Bündelladung bewirkter Effekt ist die Ver- 
schiebung der Gaußschen Bildebene. Sie wird bestimmt durch die Ver- 


®) Wir haben dabei wie üblich die durch die Bündelladung bewirkte Veränderung der 
Axialgeschwindigkeit nicht in Betracht gezogen. Gegenüber den in der Praxis auftretenden 
hohen Beschleunigungsspannungen kann dieser Effekt vernachlässigt werden. 

10) Wir setzen dagegen voraus, daß das Linsenpotential zwischen Ding- und Bild- 
ebene nirgends verschwindet und damit, daß die Axialkomponente der Geschwindigkeit 
in diesem Bereich nirgends Null wird (Spiegelwirkung). 
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schiebung der ersten auf z = z folgenden Nullstelle der modifizierten Lösung 
“ t(z). Nach Abb. 3 ergibt sie sich zu 


Mit r = ¢ und (2) erhalten wir weiter aus (17) und (19): 


v:yo, 


Vo 


(> 0). (20) 


Es ergibt sich also eine RESCUE der Bilddistanz. 
In der neuen Gaußschen Bildebene wird sich eine Veränderung des Bündel- 
durchmessers gegenüber dem in der alten ergeben. Der Abb. 3 können wir die 
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Abb. 3. Die Veränderung der charakteristischen Lösungen ¢(z) und R-s(z) (voll aus- 
gezogen) durch die Wirkung der Bündelladung (gestrichelt) 


Beziehung RV = Rs, + RsiAz, entnehmen. Indem wir s, aus (17) einführen, 
erhalten wir für die relative Veränderung der Vergrößerung die Formel 


1 F(Rs,z) t 


Wenn wir vom zweiten Glied in dieser FUN absehen, das sich in der Regel 
nur wenig auswirkt, setzt sich die Änderung der Vergrößerung aus einem posi- 
tiven und einem negativen Summanden zusammen, da das Vorzeichen von s(z) 
an der Stelle z = zy wechselt. Durch die Wirkung der Raumladung muß es also 
keineswegs in allen Fällen zu einer Verbreiterung des Bündels kommen. Das 
ist auch anschaulich klar: der Raumladungseinfluß vor der Fraunhofer-Ebene 
bewirkt eine Hebung der Randbahn R - s(z), also, wenn man nur ihn allein in Be- 
tracht zieht, eine Verengung des Bündels hinter der Fraunhofer-Ebene. 
Die Bündelladung im Bildraum hat die gegenteilige Wirkung. Das Endergebnis 
ist ein Kompromiß zwischen beiden Effekten. Wir werden im Abschnitt 5 auf 
diese Erscheinung zurückkommen. Beim Zutreffen der Bedingung 
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wird sich in der neuen Gaußschen Bildebene keine Verbreiterung des Bildflecks_ 


gegenüber dem alten ergeben. 


4. Verschiebung der Fraunhofer-Ebene und Veränderung des Brennflecks 
Auch in der dureh 
s(zp) = 0 (23) 
definierten Fraunhofer -Ebene z = zy treten durch die Wirkung der Eigenladung 


Veränderungen auf. Zunächst ergibt sich eine Verschiebung der Fraun- 


hofer-Ebene (vgl. Abb. 3): 


Asp ty D, 


Azy = —Asp/R sp = _ (24) 
0 
Mit (17) erhalten wir 
ty V2, F(R 8, 2) 
Azr = Tne, | Vo dz (>0). (25) 


20 


Durch die Wirkung der Eigenladung vergrößert sich also die Entfernung Objekt- 
ebene — Fraunhofer-Ebene. Außerdem bewirkt die gegenseitige AbstoBung der 
Elektronen im Bündeleine Veränderung des Brennflecks. Nach (10) tritt 
in der Fraunhofer-Ebene die Airysche Intensitätsverteilung auf. Die Größe 
des Brennflecks kann durch das zentrale Beugungsscheibchen, das bis zum ersten 
Minimum dieser Verteilung reicht (und in dem etwa 84% des gesamten Strah- 
stroms konzentriert sind, charakterisiert werden. Der Radius r, des zentralen 
Beugungsscheibchens ohne Eigenladungseinfluß ist durch 

bp ry = _ 3,83 (26) 

ty 

gegeben. Eine analoge Beziehung gilt aber auch für das durch die Eigenladung 
modifizierte Scheibchen in der neuen Fraunhofer-Ebene: 


by 7, = _ 3,88. (27) 
tr 
Aus (26) und (27) erhalten wir 

(28) 


Mit (17) ergibt sich die relative Veränderung des zentralen Beugungsscheibschens 
in der neuen Fraunhofer-Ebene: 


Ar, = 
ro 47V ®, yo t 
2 
Hier läßt sich von vornherein keine allgemeine Aussage über das Vorzeichen dieser 
Veränderung machen. Das zentrale Maximum wird dabei nach (10) eine Ver- 
änderung erfahren, die durch 


on 


dz + 9) 


bestimmt ist. 
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Die Beziehungen (19—30) für die Eigenladungseffekte gelten allgemein und 
hängen nicht von den speziellen Feld- und Größenverhältnissen ab (wir haben von 
der Verteilung (5), die nur ein Illustrationsbeispiel ist, bisher noch keinen Gebrauch 
gemacht). Voraussetzung ist jedoch, daß wir uns auf den paraxialen Bereich 
beschränken und die Raumladungseffekte als kleine Störungen aufassen 
dürfen. 


5. Näherungsweise Auswertung der gewonnenen Formeln. Obere Schranken 
für die Raumladungseffekte 


Um die eben hergeleiteten Formeln auszuwerten, müssen wir vor allem die 
durch (13) definierte Funktion F(r,z) bestimmen, indem wir über die Raum- 
ladungsverteilung (in unserem Beispiel über 5) integrieren. Schon die Verteilung 
(5) ist durch eine komplizierte Reihe dargestellt. Es ist von vornherein kaum zu 
erwarten, daß die Integration zu übersichtlicheren Ausdrücken führen wird. Die 
Integration läßt sich zwar durchführen (s. Anhang), doch ist das Ergebnis alles 
andere als einfach. Für bestimmte spezielle Felder (insbesondere für die kurze und 
schwache Magnetlinse und für den feldfreien Fall) kann man die strenge Rechnung 
mit einigem Aufwand weitertreiben. Es wird sich also empfehlen, nach plausiblen 
Näherungen für F(r,z) zu suchen. Dies kann in mannigfacher Weise geschehen. 
Man wird sich dabei von der physikalischen Bedeutung von F(r, z) leiten lassen: 
F(r, z) bedeutet, wie schon gesagt, die gesamte, auf die Längeneinheit der z-Achse 
bezogene Raumladung, die sich in einer flachen, achsensenkrechten zylindrischen 
Scheibe vom Radius r an der Stelle z befindet. 

Wenn wir z.B. die Funktion F(t,z) näherungsweise bestimmen wollen, be- 
denken wir, daß {(z) (wie aus Abb. 1 hervorgeht) in z = 0 die Steigung Eins hat 
und schnell zu einem Maximum ansteigt, das die Größenordnung von d (in der 
elektronenoptischen Praxis etwa 1mm) hat. Die Lösung ¢(z) wird also sehr 
schnell — etwa bei a = z + R— aus der geometrisch-optischen Bündelbegrenzung 
austreten, sich wegen des kleinen Biindelradius (R = 10-*d bis 10-*d) hoch 
dariiber erheben und im Bildraum, etwa an der Stelle 


b=2—V Rit (31) 


wieder in das Bündel zurückkehren. Zwischen z = a und z = 6 umschließt daher 
die Scheibe mit dem Radius ¢(z) praktisch die gesamte Raumladung. Für diesen 
Bereich haben wir also zu setzen: 


Fit, z) = R?20, D/P (a <z <b). (32) 
In den Randabschnitten setzen wir dagegen F(t, z) proportional zu ¢°(z): 
Fit, z) = 272 (0<z<a) 
(33) 


70, (Py\3/2 
Diese Annäherung wird um so besser sein, je kleiner die anfängliche Bündelöffnung R 
gegenüber der Halbwertsbreite d ist. Wenn wir die Randabschnitte vernach- 
lässigen und annehmen, daß ¢(z) überall die gesamte Raumladung umfaßt, er- 
halten wir einen zu großen Wert für F(t, z), den wir später zur Abschätzung brau- 
chen werden: 


b<z< 4). 


F(t,2) < . (34) 
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Auf ähnliche Weise können wir uns brauchbare Näherungsformeln für F(Rs; z) 
verschaffen. Zweckmäßigerweise wird man von einem Ansatz der Gestalt 


F(Rs,z) = (Rs)? (35) 


ausgehen Für z = 2, z, und zp hat man die durch (8—10) zusammen mit (11) 
gegebenen Raumladungsdichten fiir o(z) einzusetzen. Die Zwischenwerte kann 
man linear oder parabolisch interpolieren. Zur Abschätzung kann 


F(Rs,z) < (Rs)’noy (36) 


F (Rs, z) > (Rs)? V D/P 


und 


dienen, wobei 0,7 und o, die maximale, bzw. minimale Raumladungsdichte © 


in dem betrachteten Feldabschnitt darstellen. In der Regel wird es sich dabei um 
einen der Werte an den Stellen z = z,, %, zp handeln, der dann nach (8—10) be- 
stimmt werden kann. Die eingangs zitierten, vom geometrisch-optischen Stand- 
punkt ausgehenden Arbeiten nehmen durchwegs an, daß sämtliche Elektronen 
innerhalb der von den geometrisch-optischen Randstrahlen begrenzten Röhre 
verbleiben, verwenden also den Ansatz F — R?r 0, und haben daher die Bedeutung 
von Abschätzungen nach oben hin. 

Nun haben wir die als bekannt vorausgesetzten Lösungen s(z) und ¢(z) des 
zugeordneten geometrisch-optischen Problems, zusammen mit einer der Nä- 
herungen (32—36) in die Formeln (19—30) einzuführen und die Integrale nach 
einem passenden Verfahren auszuwerten. Im folgenden sollen einige handliche 
Formeln zur Abschätzung der Eigenladungseffekte angegeben werden. 

Betrachten wir zunächst die Verschiebung der Gaußschen Bildebene. 
Über den Verlauf des elektrischen Linsenpotentials in (20) können wir nichts all- 
gemeines aussagen und schaffen daher den entsprechenden Ausdruck mit Hilfe 
des Mittelwertsatzes aus dem Integral heraus. Mit der guten Näherung (32—33) 
erhalten wir dann, wenn wir die gesamte Strahllänge (z,—2)) mit Z bezeichnen: 


®,, = OS (38) 


mit 


Bei einer reinen Magnetlinse haben wir fiir®,,®, und®,, die Strahlspannung U 
einzusetzen. 

Mit der Abschätzung (34) gelangen wir zu einer oberen Schranke fiir die 
Bildverschiebung, die ganz allgemein gilt. Wir wollen unter ®, nunmehr den 
kleinsten Wert verstehen, den das elektrische Potential innerhalb des Abbildungs- 
feldes annimmt und nach (11) die Gesamt-Strahlstromstärke 


i = (0) (39) 


einfiihren. Es ergibt sich auf diese Weise: 


Az, | 
7 < 152-10 V%,. (40) 


Wir wollen kurz ein konkretes Zahlenbeispiel für (40) angeben, das sich mög- 
lichst eng an die in der Übermikroskopie auftretenden Verhältnisse anschließt. 


V 
ni 
ni 
| el 


Ög- 
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Das Abbildungssystem sei eine reine Magnetlinse vom Glaserschen Glockenfeld- 
typus mit k®—1. Ferner wählen wir U = 100 kV und in der Dingebene einen ° 
Bündelradius R = 10"? d (d ist die Halbwertsbreite des Feldes)'!). Bei einer 
Vergrößerung V =—-10 ergibt sich dann nach der genaueren Formel (37) eine 
Bildverschiebung von Az,/Z — 4,79 - 10-27. Bei der hohen Vergrößerung V = —100 
erhalten wir Az,/Z = 4,787: Die nach (40) ermittelten oberen Grenzen weichen 
nicht stark davon ab und betragen 4,81 - 10-27%, bzw. 4,8 7. 

Untersuchen wir nun die Bündelverbreiterung in der Bildebene. Zu- 
nächst sei bemerkt, daß der Summand sj Az, /V in (21) in der Regel kaum ins 
Gewicht fallen wird. (In dem 
eben gegebenen Zahlenbeispiel be- 
sitzt er die Größenordnung 10~*7.) tz 

Um die physikalische Be- 
deutung der Kompensations-- 
bedingung (22) zu erkennen, t 
betrachten wir speziell eine kurze 
und schwache Magnetlinse. 

Die Linse befinde sich bei z — z,,, z; 
die Objektebene bei z, — 0. Mit — 

f= %p—2y bezeichnen wir die Abb. 4. Strahlengang in der kurzen und schwachen 
Brennweite!?). Die Bahnen sind Magnetlinse 

im Dingraum als Gerade zu be- 

trachten und erfahren nur in der unmittelbaren Umgebung des Linsenortes z,; eine 
Richtungsänderung. Für den Dingraum haben wir also zu setzen (Abb. 4): 


(41) 
für den Bildraum dagegen: 


e= (zp 2); t= (21 =)». (42) 
In der Nähe der Ding- und der Bildebene befindet sich innerhalb der durch R-s(z) 
begrenzten Röhre fast die gesamte Raumladung R?ro,. Eine Abweichung davon 
tritt in der Umgebung der Fraunhofer-Ebene auf. Wir gelangen schnell zu 
einem annähernden Ergebnis, wenn wir annehmen, daß die Beiträge zyr bis zp 
und zp bis zy + f zum Integral (in Abb. 4 schraffiert) in (21) einander kompen- 
sieren. Es bleiben also nur die Beiträge von Ding- und Bildraum. Mit (41) und (42) 
gelangen wir so zu der Formel 


A = 1,52. 10°, [In — ]- (43) 


Beim Verschwinden des Ausdrucks in der eckigen Klammer wird in der Bildebene 
keine Verbreiterung des Bündels auftreten. Das ist für V, ~ — 4,8 der Fall. Dieser 
Wert ist nicht sehr exakt, da die Neigung von ¢(z) in den einander kompensierenden 
Integrationsstücken unberücksichtigt geblieben ist. Es zeigt sich aber, daß das 


11) Dieses Beispiel weicht insofern von dem in Abb. 2 dargestellten ab, als dort R = 
10-5 d gewählt wurde. 

12) Da die Linse ,,kurz‘ ist, gilt hier die Newtonsche Linsengleichung und die 
Brennebene fällt mit der Fraunhofer-Ebene zusammen. Vgl. W. Glaser u. O. Berg- 
mann, ZAMP. I., 363 (1950). ' 
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Integral in (21) gegeniiber Anderungen der Raumladungsverteilung nicht allzu 
' empfindlich ist und daß auch andere Abschätzungen zu kritischen Werten zwischen 
V, =—4 und —6 führen. Wir können also für die kurze und schwache Magnetlinse 
folgern: bei Vergrößerungen, die wesentlich über eine fünffache hinausgehen, wird 
sich eine Verbreiterung des Bündels in der Bildebene einstellen. Bei kleineren 
Vergrößerungsmaßstäben, oder gar bei Verkleinerungen wird der dingseitige Raum- 
ladungsteil den Ausschlag geben und in der Gaußschen Bildebene kommt es zu 
einer Verengung des Bündels. Benutzt man also die Linse zur Verkleinerung 
(Konzentrierspule, Kondensor), so unterstützt die Raumladung die Absichten des 
Konstrukteurs, der ja an möglichst kleinen Bildflecken interessiert ist. 


Nach dieser Plausibilitätsbetrachtung gehen wir zur Aufstellung von Ab- 
schätzungsformeln über, die für beliebige elektrische und magnetische Ab- 
bildungsfelder gelten. Es liegt nahe, dazu die Lösungen s(z) und ¢(z) in geeigneter 
Weise durch Tangenten, bzw. Sehnen zu ersetzen. Bei reinem Magnetfeld können 
wir uns dabei auf die aus (1) folgende Tatsache stützen, daß diese Lösungen stets 
gegen die Feldachse hin gekrümmt sind. Beim elektrischen Abbildungsfeld muß 
dies jedoch nicht der Fall sein. Wir gehen daher zunächst mit der bekannten 
Transformation ®) 


r(z) = R(z) - (44) 
von (1) zu der Differentialgleichung (1) zur Gleichung 
” 3 
R”=—QR mit Q(@)= [55 5) + und Q20 (45) 


über. Von den partikulären Lösungen S(z) = ®~/* s(z) und T(z) = Dt t(z) 
können wir nun wieder aussagen, daß sie 


a“ stets gegen die Achse hin gekrümmt sind. 

“2 Wir setzen mittlere bis hohe Vergröße- 

4 a’ | Fiz es rungen voraus, nehmen also an, daß die 
<T, N 8 


Verbreiterungstendenz iiberwiegt. Nun 
Kz führen wir S(z) und 7(z) in (21) ein und 
ersetzen diese Funktionen nach Abb. 5 


md im Bildraum durch ihre Tangente (Sz 
und 7,) und im Dingraum durch ihre 
P Sehnen (Sp und 7'p). Wir haben dann 


die Garantie, im ersten Fall zu große, 
im zweiten zu kleine, im ganzen also zu 
große Werte zu erhalten, wenn wir die 
Abschätzungen (36) benützen. Indem 
wir mit (44) und (18) wieder zu s(z) und ¢(z) zurückkehren, ergibt sich die 
allgemeine Abschätzung: 


Abb. 5. Zur Abschätzung der Raum- 
ladungseffekte 


: (2 — Om?ptr 
AV M_ 1/4 46 
24 .| ® (PD, Py) | (46) 


13) Vgl. auch die Abschätzung von oberen Grenzen des Farbfehlers, W. Glaser, 
Z. Phys. 116, 56 (1940). 


Fi 


Fi 


| Fi 

0% 
Fi 

mi 
rel 

V 

Se 

Di 
be 
du 

Fi 

Sp N 

In 

; 

2, Zr N tre 

Ai 

“ 

sti 
ha 
lac 

da 

vo 

St 


ser, 


H. Grimm: Die Eigenladungs-Verbreiterung eines fokussierten Elektronenbiindels 143 


Für das rein magnetische Abbildungsfeld, und wenn wir 0y = 9) und @,, = 0, = 
o/|V| setzen dürfen, vereinfacht sich (46) zu 


> i _ 2rtr 
AV <38.10 - (47) 


Für das auf $. 140 angeführte Zahlenbeispiel, das den Verhältnissen in der Uber- 
mikroskopie entspricht, ergibt sich bei einer Vergrößerung von V =—10 eine 
relative Bündelverbreiterung AV < 1,23 - 10*i und bei der hohen Vergrößerung 
V =—100 der Wert AV < 1,91. 108%. 

Für die Verschiebung der Fraunhofer-Ebene erhalten wir eine obere 
Schranke, wenn wir in (25) F(Rs,z) durch die gesamte Raumladung ersetzen: 


Azp 
= < 1,52 - 10% Ro, V®r- (48) 
Für ein rein magnetisches Abbildungefeld ergibt sich daraus: 
1,52 10: (49) 


Die relative Verschiebung der Fraunhofer-Ebene ist bei V—=—10 nach (49) 
beschränkt durch Azp/zp < 1,56 - 10°?i und bei hoher Vergrößerung (V = —100) 
durch Azp/zp < 1,13 - 1037. 

Nun haben wir noch die Veränderung des Brennflecks in der neuen 
Fraunhofer-Ebene zu bestimmen. 

Dazu schätzen wir F(t,z) mit (32—33) ab. Es ergibt sich: 


1 Rz 2 \ 
+ 27 (opin 24+ |. 50) 
D,, VP, E ( r) ( 
Fiir ein rein magnetisches Abbildungsfeld ergibt sich speziell: 
Rzr 
Ary < 1,52 - 108 +7 +# (In — )|- (51) 


In unserem Zahlenbeispiel liefert der Ausdruck tr Azp/tp praktisch keinen Bei- 
trag und für V = —10 erhalten wir Ar, < 7,26 - 103. Für V = —100 ergibt sich 
Ar, < 5,41- 

Nimmt man an, daß durch die Blende vom Radius R = 10-*d ein Strahl- 
strom von 10-7 A tritt!) und setzt man die Halbwertsbreite d= Imm, so er- 
hält man in unserem Zahlenbeispiel ‘die folgenden Abschätzungen für die Raum- 
ladungseffekte (V = 100): . 


o <4,8-10-7; AV < 0,19 


y 


< 1,13. 104; Ary< 5,41. 10-10, 


14) B v. Borries, Hab. T. H. Berlin 1945, Optik 8, 321, 389 (1948) gibt an, 
daß bei U = 90kV ein engster Radius des vom Elektronenstrahler erzeugten Bündels 
von r=0,02mm und ein Strahlstrom von 50 yA auftritt. Das entspricht einer 
Stromdichte von etwa 4. 10-2 A/mm?., 
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Anhang 
Bestimmung der Funktion Fer, z) 
Die Funktion F(r, z) ist definiert durch (13) und lautet mit (5) und (7) 
[To 1 n )y2n+2 
F(r,2)= 220) la (2) ea (2) 


n=0 
(52) 
(n-») +2 (2) dx, 
wobei wir br — x gesetzt haben. Die darin auftretenden Integrale sind vom Typus 
S a de (53) 


und lassen sich mit der von Lommel®) angegebenen Formel 


I) 1 1 | | 


im Prinzip leicht auswerten. Wir verzichten darauf, den sich ergebenden umständ- 
lichen Ausdruck anzuschreiben. 


») E. Lommel, Math. Ann. 14, 530 (1897). 


Wien 6, Joanelligasse 7/9. 


(Bei der Redaktion eingegangen am 26. Juni 1952.) 
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Zur Umkehrbarkeit der Lichtwege in geschichteten Medien 
Von Hubert Pohlack 


Inhaltsiibersicht 


Das optische Durchlässigkeitsvermögen einer beliebigen Anordnung geschich- 
teter Medien mit planparallelen Begrenzungsflächen ist unabhängig vom Richtungs- 
sinn der senkrecht auftreffenden Lichtstrahlung, wenn die Außenmedien reelle 
Bre:hzahlen oder den gleichen Absorptionsindex besitzen. Unterscheiden sich die 
beiden Absorptionsindizes, so gilt der Umkehrungssatz nicht. Es wird gezeigt, daB 
dieser letzten Aussage jedoch keine optisch-meßteshnische Bedeutung zukommt. 
Insbesondere wird nachgewiesen, daß eine einseitige „Entspiegelung‘‘ eines abge- 
schlossenen, in einem einheitlichen Medium befindlichen optischen Schicht-Systems 
nicht zum unterschiedlichen Durchlässigkeitsvermögen führen kann. 


In zwei neueren Veréffentlichungen')*) hat v. Fragstein das Problem des 
Energieübergangs elektromagnetischer Strahlung an der Grenze zwischen ab- 
sorbierenden Medien behandelt und in diesem Zusammenhang die Reziprozität 
der Strahlungsdurchlässigkeit einer solchen Grenzfläche untersucht. Die am 
Schluß der einen Arbeit aufgeworfene Frage, wie es sich mit der Umkehrbarkeit 
der durchgelassen:n Strahlungsenergie im Falle komplizierterer Schichtanordnun- 
gen verhält, gibt mir Veranlassung, die Ergebnisse einiger Untersuchungen 
mitzuteilen, die ich kürzlich auf eine ähnliche Fragestellung hin durchgeführt habe. 


Problemstellung 


Eine beliebige Anordnung planparalleler optischer Schichten wird von einer 
ebenen Lichtwelle durchsetzt. Gesucht sind die allgemeinen Beziehungen zwischen 
den elektromagnetischen Feldgrößen unmittelbar vor und hinter den 4: Beren 
S:hichtgrenzen, wenn das System von räumlich ausgedehnten optischen Medien 
eingeschlossen ist. Das Intensitätsdurchlässigkeitsvermögen ist zu definieren 
und der Einfluß der Strahlungsumkehr auf diesen Term zu prüfen. Schließlich 
ist zu untersuchen, inwieweit die formalen Voraussetzungen für das Resultat 


.mit gewissen idealisierenden Annahmen über Strahlungsquelle und -empfänger 


vereinbar sind. 


Die Wellenausbreitung in Interferenzschichtsystemen 
Wir gehen von der üblichen Betrachtungsweise aus und setzen eine ebene, un- 
endlich ausgedehnte, linear polarisierte Lichtwelle an, die sich in Richtung der 


1) v, Fragstein, Ann. Physik (6) 7, 63 (1950). 
2) v.. Fragstein, Phys. Bl. 6, 153 (1950). 


Ann. Physik. 6. Folge, Bd. 11 
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z-Achse eines rechtwinkligen Koordinatensystems ausbreite. €- und $-Vektor 
mögen in der x- bzw. y-Richtung schwingen®). Dann kann der Schwingungs- 
vorgang in jeder Fläche z = const des ausgedehnten Mediums 0 (mit der im all- 
gemeinen komplexen Brechzahl ty = ny (1— j %)) wie folgt beschrieben werden: 


Ein = eo 

=0 

(2) = 0 2; 

9! = (A auf Vakuum bezogen) 
2) — () 
0 

= n, ei (ot-#,2) 

=0 


mit ¢, und e, als Einheitsvektoren in der 2- bzw. y-Richtung und Z als Amplitude 
des elektrischen Vektors in der Fläche z = 0, oder, wenn man zeitfreie Größen be- 
nutzt und die Vektorkennzeichnung sowie die durchweg unveränderliche Indi- 
zierung für die Polarisationsrichtungen wegläßt: 


> 


H, =n, Ey e’%:. 


Bringt man jetzt im Medium 0 senkrecht zur Strahlungsrichtung ein System 
dünner, „interferenzfähiger‘“, planparalleler, isotroper Schichten 1, 2,..., m an 
(mit den im allgemeinen komplexen Brechzahlen n,,...n,) und schließt diese 
Anordnung ab durch das unendlich ausgedehnte Medium m + 1 (n„..ı), so über- 
lagern sich in den Schichten 1,..., m die in positiver und negativer z-Richtung 
laufenden Teilwellen (Reflexionen), und man erhält auf Grund der Stetigkeits- 
bedingungen die bekannten Terme für die komplexen Amplituden in den einzelnen 
Grenzflächen: 


Ky), Ky +, = + 
Ho = n, (EY — = u, B,) 
e itd + E, etitid, — 4 k, 
y, (ET B,) (1). 
Enm ° kne j tn din t E,, et) fudm = 
Ham i= Um (Ey, e itm dm — E,,, e* tmdm) 


*) Annahmen, die gegenüber natürlichem Licht für isotrope Medien keine Beschrankung 
der Allgemeinheit bedeuten. 
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Hieraus gewinnt man fiir die Beziehungen zwischen benachbarten Grenz- 
flächen die folgenden linearen Gleichungspaare: 


Eon = cos hd + sinhd 
Hop = jn,sin cos Hy) 


= sin ft, + cos Ho); (2) 
Em -ijn = COS Em Emim+ i+) sin +1 


H„n-ım = 7 Mn sin Em Am Emm +ı 4 cos £,, dy, * 
oder allgemein fiir die Grenzflächen der i-ten Schicht: 


E;-ıji = KO E;jizı + H;jirı 


3 
Ay = KY + RD 
mit 
RKO = KYO = cos f,d, 
RO = j sin 
1; 
= jn, sin 
Ks erweist sich als zweckmäßig, die Gleichungspaare (3) in Matrizenform zu 
schreiben ®): 
| 0 | oF aL. 


Dann vereinfacht sich nämlich, wenn man an den Beziehungen unmittelbar 
zwischen den Feldgrößen der ersten und letzten Grenzfläche interessiert ist, die 
Elimination der Zwischengleichungen zur formalen Matrizenmultiplikation: 


Ey), 0 E 10 
| mim 
so daß mit den ne der Produktmatrix 
Mor | 
zu setzen ist: 
Kon = Mn rit Rig Ham +1 
(0) 


Hy), Ma Moe H im 


4) Vgl. u.a.: H. Weinstein, J. opt. Soc. Amer. 87, 576 (1947); H. Pohlack, Ann. 
Physik (65 5,311 (1950); F: Abelés, J. Physique Radium 11, 307 (1950). 
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Bei Benutzung der ersten und letzten Terme des Gleichungssystems (1) kann (5) 
dann auch geschrieben werden: 


Ey + Eo = (Ku + 
m (Eg — Eo) = (War + Thm +1 Ra) (6) 


oder nach Elimination von E,: 


Ei = + “2? Ei (7) 


Der Umkehrungssatz fiir Interferenzschichtsysteme 
Die aus Gl. (7) gebildeten Quotienten 


E; Ry + + — 2 a+ 
(8) 
Me 
= 
und 
+1 _ +1° R41 
H+ My Et 
5 Mm +1 8 
- (83) 
Vit, + Rn + = = + Var 9 


stellen die Verhältnisse zwischen den RATE und einfallenden Feld- 
größen dar und können als komplexe Durchlässigkeitskoeffizienten der Feld- 
stärken bezeichnet werden. 

Das Intensitätsdurchlässigkeitsvermögen wird dann definiertals das Ver- 
hältnis der in das Abschlußmedium eintretenden Strahlungsenergie zu derjenigen 
Energie, welche durch eine gedachte Fläche am Ort der vorderen Schicht- 
begrenzung strömt, wenn das Schichtsystem beseitigt wird und das vor- 
gelagerte Medium den gesamten Raum einnimmt). 

Die „mittleren Energieströmungen‘ ergeben sich durch Integration des Vektor- 
produkts der zeitabhängigen Feldstärken über die Schwingungsperiode®). Man 
kann aber auch, da bisher bereits mit den zeitfreien Größen operiert wurde, 
auf diese Mittelwertbildung verzichten und bekommt die gleichen Ergebnisse, 
wenn man (in Analogie zur Wechselstromtechnik) die ‚effektiven‘‘ Feldstärken 


Bett V2 | Emax| 


5) Diese etwas umständliche Ausdrucksweise ist der üblichen Formulierung (,,ein- 
stromende“ Energie) vorzuziehen, da sich bei absorbierendem Außenmedium die Energie- 
strömungsverhältnisse unmittelbar vor der Grenzfläche verkomplizieren (vgl. v. Frag- 
stein, a.a.O.). Wir werden jedoch im folgenden aus Gründen der Einfachheit die 
kürzere Formulierung verwenden, ohne Unklarheiten befürchten zu müssen. 

6) Vgl. u. a.: Born und Ladenburg, Physik. Z. 12, 198 (1911), sowie die zitierten 
Arbeiten v. Fragsteins. 
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und 


1 N 
Hess = y2 | Hmax! 
mit dem cosinus der Phasendifferenz multipliziert’). 
Wir betrachten zunächst den Fall nichtabsorbierender Außenmedien: Wegen 
H = nE wird die Phasendifferenz Null, so daß dann mit (8) (8a) für das Inten- 
ee zu schreiben ist: 


| 
Dojm+1 = | 
(9) 


Es wird jetzt nach den ,,Ubertragungseigenschaften“ des Schichtsystems ge- 
fragt, wenn man die Strahlungsrichtung umkehrt. Behält man das bisher benutzte 
Koordinatensystem bei, so lautet die Matrizengleichung entsprechend (4): 


Es): 0| ell-| Emm+ı 0| 
Das negative Vorzeichen bei H* ergibt sich dabei aus der Umkehrung der Fort- 
pflanzungsrichtung, da ja der $-Vektor immer durch positive 90°-Drehung aus 
dem &-Vektor hervorgeht. Im übrigen ändert sich an den Beziehungen, wie sie 
von den aus den Stetigkeitsbedingungen hergeleiteten Übertragungsgleichungen 


beschrieben werden, nichts, insbesondere ist die Produktmatrix ||@|| in (4a) 
identisch mit dem Matrizenprodukt in (4): 


| = KO) KOM). 
Der Matrizenkalkiil erlaubt nun die folgenden Umformungen der Gl. (4a): 
I— m/m+1 0 [8] — Ru 


(mit | 8 | als Determinante der Produktmatrix | 8 |), wofür auch geschrieben werden 
Roe 


kann: 
| 
Ko Ku | 


Mit den entsprechenden Überlegungen wie oben folgt dann aber für die ,,In- 
tensitätsdurchlässigkeit‘‘ in Richtung (m + 1) —0: 


(4a) 


Ein 0 
am 0 


| Enm+ı ol = 
0 (ei |X| 


0 
0 . 


D | Eo 
m+1/o = 
\2 


und Dy, 41/0 unterschieden sich also nur durch den Faktor | |*. 


?) Die „effektiven“ Feldstarken entstehen natürlich auch erst durch zeitliche Mittel- 
wertbildung. 
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Es soll jetzt mittels Induktionsschluß gezeigt werden, daß die Determinante 
|| der Produktmatrix immer gleich 1 ist: 


Für die Einzelschichtmatrix war gefunden worden: 
| =| 
|jn,sinf,d, cost,d, 


. sin f,d 


Um den allgemeinen Fall (beliebige Kombinationen absorbierender und nicht- 


absorbierender homogener Schichten) zu erfassen, haben wir die Wellenzahl 


komplex angesetzt: 
= k, x) = (IJ x). 
Wegen 
cos f,d; = cos [k; (l1—jx,)d,] = cos k,d, -Cojk,x,d, + jsink,d,- Sink, x,d, 


und 


kann die Schichtmatrix dann auch geschrieben werden: 


|KO || = cos p; + jsing;: Sin gy; iz (sin Coj y; —j cos Sinp) 


jn, (sing, -Coj cos p; -Coj + j sing; - Sin g,) 
(y= kd, 9, — kh, x, 
Damit erhält man für die Determinante 
= (cosp, -Coj + j sin y; Sin g,)? + (sin y; — j cos y; Sin g,)* 
= (sin? + cos* p,) (Coj? Sin? = 1. 
Mit der Voraussetzung 
||| || || | 1) GREY | | l<k<m 
wird aber auch 
= | | | +D) | 


gleich eins, was zu beweisen war, und es folgt aus (9) und (10): 


Dojm+1 = 


Wir formulieren mit diesem Ergebnis den folgenden ,,Umkehrungssatz“: Fällt 
eine ebene elektromagnetische Welle senkrecht auf ein planparalleles Interferenz- 
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schichtsystem, das von nichtabsorbierenden Außenmedien eingeschlossen wird, 
so ist das Verhältnis der durch die Schichten hindurchgelassenen zur auftreffenden 
Strahlungsintensität unabhängig vom Richtungssinn der Strahlung®). 


Es wird jetzt die Umkehrbarkeit für den Fall untersucht, daß mindestens eins 
der beiden Außenmedien absorbiert. Definiert man wieder die einfallende und 


durchgehende Intensität als proportional zu &, x 9, bzw. ©, x 9, oder (wie oben) 
als Produkt der effektiven Feldstärken und dem cosinus des Phasenwinkels, so ist 
unter Benutzung der Gln. (8) (8a) fiir die Durchlässigkeiten zu schreiben: 


a) in Richtung 0-—>m + 1 


| 008 +1 Hm +1) 


Et Ht ‘cos (Ei, 

4 

VE + Vito ty +1 + + 

. 008 (En +1» Hin +1) (11) 
eos(E}, Hi) 
1+ 
V: 

Var, Kn + Viton, +1 Re + Ny 


b) in Richtung m + 1 —0 


Es Hs cos (Es, Be) 
En+ı Hn+ı 008 (Ey, +7, Hm +1) 
Lax? 9 
V (12) 
1+% 
+ Mn +1 Vi tm. Ry, ] Stop | 
1 
daH=nH=nyV1+ x? E-eiartgx, also cos (E, H) =-——— -. Dem- 
(B, H) = 


nach ist die ‚Durchlässigkeit‘ nicht invariant gegen Seitenvertauschung, wenn 
eins der beiden Außenmedien absorbiert oder sich die beiden Absorptionsindizes 


‘(%, %m+1) voneinander unterscheiden®). Für den Spezialfall verschwindender 


Dicken des Schichtsystems folgt die Nichtumkehrbarkeit für einfache Grenz- 
jlichen zwischen zwei Außenmedien mit unterschiedlichen Absorptionsindizes 
(vel. v. Fragstein, a.a.0.). 


®) Herleitung und Formulierung stehen in weitgehender Analogie zum Umkehrungs- 
satz für das Betriebsübertragungsmaß elektrischer Vierpole (vgl. auch Anm. 9). 

®) Wenn der analoge Umkehrungssatz in der Vierpoltheorie auch im allgemeinen Fall 
komplexer Außenwiderstände Gültigkeit besitzt, so liegt dies daran, daß hier Schein- 
leistungen verglichen werden und somit der ,,Leistungsfaktor“ (also der cosinus der „ 
Phasendifferenz) unberücksichtigt bleibt, während die entsprechende Größe bei den uns 
interessierenden „Wirkintensitäten‘“‘ gerade für die Nichtumkehrbarkeit verantwortlich 
ist. (Vgl. Anm. 8.) 


a 
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Beliebige Schichtungen (‚„‚dieke‘“ Schichten) 


Bisher wurden interferenzfähige Schichten behandelt. Die Zusammensetzung 
der Schichtsysteme kann nun dadurch weiter verallgemeinert werden, daß 
man auch „dicke“ (nicht interferenzfähige) Schichten zuläßt, ohne daß sich die 
Voraussetzungen für die Umkehrbarkeit ändern. Das soll an einem Beispiel klar- 
gemacht werden: 


Wir denken uns eine Interferenzschichtanordnung an einer beliebigen Stelle 
aufgetrennt und durch eine dicke planparallele Schicht erweitert. Dann können die 
bisher für Interferenzschichtsysteme durchgeführten Überlegungen auf die beiden 
Einzelsysteme übertragen werden, indem man die eingefügte, nicht interferenzfähige 
Schicht als Abschlußmedium bzw. als vorgelagertes Medium der Teilsysteme auf- 
faßt. Das ist sicher zulässig, wenn die Zickzackreflexionen im Inneren der dicken 
Schicht berücksichtigt werden, die ja nicht „phasengerecht‘‘ zu summieren sind, 
so daß die Intensitäten der einzelnen Teilstrahlen betrachtet werden können. 


— — 
Seien S, und S, die beiden getrennten Einzelsysteme, Ds, und Ds. bzw. Ds, und 


Da ihre Durchlässigkeiten!®) für beide Strahlungsrichtungen, Rs, und Rg, 
entsprechend das Reflexionsvermögen (an den inneren Grenzflächen) und bedeute 
x den auf einmaliges Durchlaufen der Zwischenschicht bezogenen Dämpfungs- 
faktor (für nichtabsorbierende Zwischenschichten wird x = 1), so bekommt man 
für die Strahlungsdurchlässigkeit D};, . ı des erweiterten Systems: 


Zus 
"Ti, Rg 


Und wenn man beachtet, daß man mit (11) und (12) erhält: 


Ds, = ite Ds, und = 


(der Index e bezieht sich auf die eingeschaltete Schicht), so ist fir die 
Strahlungsrichtung m + 1—0 zu schreiben: 


Ds, Ds. m +1 


Der Umkehrungssatz besitzt also fiir die durch die dicke Schicht erweiterte An- 
ordnung unter den genau gleichen Voraussetzungen Giiltigkeit wie fiir ein einfach 
zusammenhängendes Interferenzschichtsystem. Man verifiziert die entsprechende 
Aussage leicht für beliebige andere Erweiterungen. 


e 


10) die ja gemäß den Feststellungen im vorigen Abschnitt i.a. vom Richtungssinn 
der Strahlung abhängen. 
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Kritik der Ergebnisse 


Die Notwendigkeit einer Unterscheidung zwischen absorbierenden und nicht- 
absorbierenden Außenmedien in bezug auf die Gültigkeit des Umkehrungssatzes 
erscheint vom Standpunkt der praktischen Optik aus sowie in Hinblick auf das 
Helmholtzsche Reziprozitätstheorem nicht überzeugend, und es liegt die 
Vermutung nahe, daß die bei der formalen Herleitung des Resultats gemachten 
Voraussetzungen in teilweisem Widerspruch zur realisierbaren Versuchsanordnung 
oder wenigstens zum idealisierten Gedankenexperiment stehen. 


Wenn wir, um die Reziprozität der Durchlässigkeit einer Schichtanordnung zu 
prüfen, bisher so vorgingen, daß wir uns die Strahlung einmal vom Medium ny 
aus auf die Vorderfläche, das andere Mal vom Medium n,, ,, aus auf die Rückfläche 
des Systems auftreffend dachten, so haben wir bei dieser Betrachtungsweise über 
Ort und Struktur von Strahlungsquelle und -empfänger keine Abmachungen 
getroffen, allenfalls stillschweigend vorausgesetzt, daß die Schwingungserzeugung 
und der eigentliche Meßvorgang von der zu prüfenden Anordnung ‚weit genug‘“ 
entfernt stattfinden, um unbeeinflußt von Rückwirkungen der reflektierten Ener- 
gie zu bleiben"). Nunstellt aber die Lichtquelle selbst wieder irgendein (unter 
Umständen sehr kompliziertes) System optischer Medien dar; das gleiche gilt für 
den Strahlungsempfänger. Idealisieren wir Quelle @ und Empfänger E durch je ein 
beliebiges Medium der Brechzahl n, bzw. n, mit planparalleler Begrenzung nach 
dem zu untersuchenden Schichtsystem hin, und nehmen wir weiter an, auf die 
Begrenzungsfläche von Q treffe eine ‚irgendwo im Inneren‘ erzeugte ebene Welle 
auf, so wird, wenn Q und E die Plätze vertauschen, in Wirklichkeit nicht das ur- 
sprüngliche System S, sondern der gesamte Bereich zwischen @Q und E auf Umkehr- 
barkeit untersucht. Brächte man aber Q und £ in optischen Kontakt mit S, um 
die Zwischenbereiche auszuschalten, so hätte man in Gl. (9) bzw. (11) n, und n, 
als Brechzahlen der Außenmedien einzusetzen, und in Gl. (10) bzw. (12) wären 
diese beiden Brechzahlen miteinander zu vertauschen. Dies würde aber zu einer 
von der bisherigen sich wesentlich unterscheidenden Reziprozitätsbedingung führen: 
Jetzt müßten nämlich entweder die Brechzahlen n, und n, oder die Elemente 
R,, und $,, der Produktmatrix des Schichtsystems einander gleich sein, wenn 
Umkehrbarkeit vorliegen soll'*). Andererseits wurde jedoch in der ursprünglichen 
Aufgabe vorausgesetzt, daß das zu prüfende Schichtsystem seine relative Lage zu 
den beiden Außenmedien beibehalten soll, wenn die Strahlungsrichtung umgedreht 
wird. Diese Forderung führt aber nach den obigen Überlegungen immer dann zu 
einem Widerspruch, wenn ein optischer Meßvorgang eingeschaltet wird. Lediglich 
für die formale sukzessive Durchrechnung einer Schichtanordnung (vgl. das in 
den v. Fragsteinschen Arbeiten diskutierte Beispiel einer lichtdurchlässigen 
Aluminiumschicht) ist die ursprüngliche Betrachtungsweise sinnvoll. 


1) Eine ähnliche Forderung formuliert W. O. Schumann (Elektrische Wellen, 
Carl Hansen Verlag, München 1948) bei der Herleitung des mit unserem Problem ver- 
wandten Sommerfeldschen Reziprozitätstheorems derdrahtlosen Telegraphie: ,,Stérende 
Ungleichmäßigkeiten des Dielektrikums seien immerhin so weit von der Antenne ent- 
fernt, daß deren Einfluß wegen des relativ starken Antennenfeldes sich nicht bemerkbar 
macht“, 

12) Die Bedingung n, = n, bzw. &,, = 85, führt zu interessanten Aussagen über die 
„Anpassung“ von Strahlungsquelle und -empfänger bzw. über die Symmetrieeigen- 
schaften des zu prüfenden Schichtsystems. Ich beabsichtige, die entsprechenden Er- 
gebnisse in einer weiteren Arbeit mitzuteilen. 
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In praxi liegen die Verhältnisse glücklicherweise meist einfach: Ist die zu 
prüfende Anordnung ‚abgeschlossen‘‘ und eingebettet in einem einheitlichen Medium 
( = Nm+1) so ändert sich bei Vertauschung von Strahlungsquelle Q und -emp- 
fänger E an der Energieübertragung zwischen Q bzw. E und den Einbettungs- 
medien nichts, diese können selbst als Sender bzw. Empfänger und als ,,AuBen- 
medien‘‘ aufgefaßt werden. Damit sind aber sämtliche Voraussetzungen für die 
Gültigkeit des Umkehrungssatzes erfüllt, unabhängig von der Struktur des Schicht- 
systems und der Art des Außenmediums. In diesem Sinne ist auch die Frage, ob 
durch einseitige ‚„‚Entspiegelung‘‘ ein unterschiedliches Durchlässigkeitsvermögen 
erzielt werden kann, eindeutig zu verneinen "®). 


13) Die bei gewissen optischen Anordnungen (beschichtete Brillengläser, Trennungs- 
platten für ungestörte Beobachtung usw.) erzielte „unterschiedliche Durchsicht“ steht 
nur scheinbar im Widerspruch zu dieser Aussage. Sie ist dadurch bedingt, daß die von 
der Anordnung getrennten Räume unterschiedlich hell beleuchtet und die nach beiden 
Seiten reflektierten Intensitäten verschieden sind. 


Jena, Carl-Zeiss-Werke. 


(Bei der Redaktion eingegangen am 4. August 1952.) 
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Zur Struktur des ferromagnetischen Zustandes in «-Fe © 
und einigen Fe-Legierungen 


Von Gerhard Heber 


Inhaltsiibersicht 


Es wird auf einige gegenwärtig wohl schwer zu beseitigende Schwächen der 
Zenerschen Theorie der Sättigungsmagnetisierung des «-Fe und einiger Fe- 
Legierungen hingewiesen. Es sollten deshalb weitere experimentelle Argumente 
für oder gegen diese Theorie gesucht werden, wozu Möglichkeiten gezeigt und 
Formeln bereitgestellt werden. Die möglichen Arten der dem Ferromagnetismus 
in den erwähnten Substanzen zugrunde liegenden Austauschkräfte werden kritisch 
diskutiert. 


$1. Kürzlich hat Zener!) darauf hingewiesen, daß man in überraschend ein- 
facher Weise die beobachteten Sättigungsmagnetisierungen von x-Fe sowie einiger 
Fe-Legierungen (z. T. mit Überstruktur) verstehen kann. Es handelt sich im 
wesentlichen um folgende Gesichtspunkte: 

I. Eine gewisse ‚„‚Disproportionierung‘‘ der Fe-Atome bezüglich der 3d-Elek- 
tronen findet?) im Metall statt. Sie ist dadurch bestimmt, daß die innere Aus- 
tauschenergie aller 3d-Elektronen am gleichen Atom maximal negativ für den 
ganzen Kristall wird. Dies führt z. B. im x-Fe dazu, daß 50% der Fe-Atome im 
Metall 2 3d-Elektronen mehr, 50% aber 2 3d-Elektronen weniger besitzen als im 
metallischen Zustand ohne ,,Disproportionierung“ (..... 3d’4s). Es bilden sich 
also Fe?-- und Fe*+-Ionen aus, in welchen die 3d-Elektronen folgende Spins an- 
nehmen: Fe? — ——); Fe*+(+++++). Es ist anzunehmen, 
daß die Fe-Ionen zur größtmöglichen Reduktion der Coulomb-Energie eine Über- 
struktur bilden derart, daß im Idealfalle alle Fe*--Ionen Gitterpunkte des einen, 
alle Fe®+-Ionen die Gitterpunkte des anderen einfach-kubischen Gitters besetzen 
und so insgesamt das kubisch-raumzentrierte Gitter des x-Fe bilden. — Ähnliche 
Strukturen gibt Zener für Fe,N, Fe,Si und Fe,Al an. 

II. Die Austauschkräfte zwischen den magnetischen Ionen sollen im Falle 
a-Fe und Fe,N Anlaß zu einer antiparallelen Einstellung der resultierenden Spins 
der beiden verschiedenen magnetischen Fe-Ionen-Sorten geben, so daß eine Art 
Ferrimagnetismus entsteht, wie bei Ferriten und anderen Substanzen bekannt. 
Im Falle Fe,Si und Fe,Al sollen die Austauschkräfte eine Parallelstellung der 
magnetischen Fe-Ionen bewirken. Unter diesen Annahmen werden die beobachte- 
ten Sättigungsmagnetisierungen N'(0) recht gut verständlich. 

Es scheint angebracht, im Anschluß an die Zenerschen Überlegungen einige 
Bemerkungen und Betrachtungen in die Diskussion zu werfen. 


1) C. Zener, Physic. Rev. 85, 324 (1962). 
2) Nach Zener. 
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Zunächst muß betont werden, daß die erwähnten Annahmen leider quantitativ 
nicht sehr sicher untermauert sind, was an ihrer Plausibilität natürlich nichts 
ändert. So trägt der Nachweis der energetischen Bevorzugung der , dispropor- 
tionierten Überstrukturen nur Abschätzungscharakter; Vorzeichen und Größen- 
ordnung der Austauschwechselwirkung benachbarter Ionen werden notwendiger- 
"weise (wegen des Fehlens.von Berechnungen der Austauschintegrale) nur plausibel 
gemacht. Schließlich scheint es fraglich, ob man die 3d-Elektronen des Ferro- 
magnetikums einigermaßen adäquat mit dem ionogenen Modell erfassen kann. 

Wie man sieht, sind die Schwächen der Zenerschen Überlegungen solcher Art, 
daß sie wahrscheinlich in naher Zukunft nicht mit Sicherheit überwunden werden 
können, so daß die genauere quantitative Begründung der Zenerschen Theorie 
auf beträchtliche Schwierigkeiten stoßen dürfte. Andererseits ist die experimentelle 
Basis dieser Überlegungen bisher einzig M(0) und deshalb sehr schmal. Verfasser 
möchte daher vorschlagen, zu versuchen, die experimentelle Basis zu erweitern. 
Man sollte vielleicht zuerst versuchen, Folgerungen über M(T) und (7'), spontane 
Magnetisierung M und Suszeptibilität y in Abhängigkeit von der Temperatur 7', 
aus Zeners Annahmen zu ziehen und diese mit der Erfahrung zu vergleichen. 
Die hierzu nötigen Formeln werden im nächsten Abschnitt bereitgestellt. 

Eine solche Erweiterung der experimentellen Basis würde auch in folgender 
Hinsicht von Bedeutung sein: Bekanntlich haben Rechnungen von Manning?) 
über die Strukturen der Energiebänder in x-Fe zu einem befriedigenden Verständnis 
der Sättigungsmagnetisierung in diesem Metalle geführt. Zeners Vorschlag 
ist also nicht die einzige Möglichkeit, das M(0) in x-Fe zu erklären. Um die beiden 
Vorschläge gegeneinander abwägen zu können, sollte man andere Eigenschaften 
des Metalles zu Rate ziehen; vielleicht sind M(7) und x(7') hierzu geeignet. 

Im Schlußabschnitt wird darauf hingewiesen, daß nach Meinung des Verfassers 
die kürzlich ebenfalls von Zener eingeführte) indirekte Austauschwechselwirkung 
(an welcher bekanntlich die 4s-Elektronen des Ferromagnetikums wesentlich be- 
teiligt sind) nicht bei allen diskutierten Substanzen zum Verständnis der postu- 
lierten Spinanordnung notwendig ist. 

$2. Während die von Zener!) vorgeschlagene Struktur der ferromagnetischen 
geordneten Legierungen Fe,Al, Fe,Si kein anomales Verhalten von M(7) und 
x(T) erwarten läßt, ist dies bei x-Fe und Fe,N sehr wohl der Fall. Ferrimagnetische 
Stoffe kann man ja, wie N&el°) u. a. gezeigt haben, i. a. daran erkennen, daß das 
Curie-Weißsche Gesetz für 7(7') nicht erfüllt ist und daß auch M(T) nicht 
den z. B. von Ni her bekannten und von der Weißschen Theorie des Ferromagne- 
tismus gut dargestellten Verlauf hat. Nur in gewissen Spezialfällen sind diese 
Anomalien nicht vorhanden; wir werden im folgenden zeigen, daß diese Spezial- 
fälle bei x-Fe und Fe,N nicht zufällig verwirklicht sein können. Wir werden uns 
dazu der Näherung des Molekularfeldes bedienen. Dabei ist es allerdings sehr 
fraglich, welchen Einfluß die von Zener zusätzlich postulierten Fluktuationen der 
3d-Elektronen zwischen den beiden Teilgittern im x-Fe haben; diese Frage muß 
hier offen bleiben. 

Die Molekularfeldmethode ist bereits von Néel5) u. a. auf Substanzen ange- 
wandt worden, deren magnetische Ionen auf zwei getrennten Teilgittern ange- 
ordnet sind. Dabei ist die Theorie allerdings explizit nur entwickelt für solche 


3) M. F. Manning, Physic. Rev. 68, 190 (1943). 


4) C. Zener, Physic. Rev. 81, 440; 82, 403; 83, 298 (1951). 
5) L. Néel, Ann. de phys. [12], 3, 137 (1948). 
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Kristalle, welche aus gleichartigen magnetischen Ionen bestehen; die Ausdeh- 
nung auf Stoffe, welche Ionen verschiedenen magnetischen Momentes auf den 
beiden Teilgittern besitzen, ist von Néel wohl schon durchgefiihrt und auf einige 
Substanzen angewandt worden; er hat aber darauf verzichtet, die entsprechenden 
Formeln anzugeben. Es scheint deshalb berechtigt, diese Verallgemeinerung hier 
in Kürze anzudeuten. 

Die Ionen der Sorte a mögen das atomare magnetische Moment u, = 28, “pz, 
die der Sorte b jedoch u, = 2S, u besitzen. wg sei das Bohrsche Magneton, 
S, und S, die Spinmomente von a bzw. 6 in Einheiten %/2. Ein Bahnmoment 
sollen die Ionen nicht besitzen. M,(T), M,(7) seien die zur Temperatur 7 ge- 
hörigen spontanen Magnetisierungen aller N,a- bzw. N„b-Ionen pro Volu- 
meneinheit; wir nennen u, N, = Mao, Hy Ny = Myo. Das auf eines der a-Ionen 
wirkende Molekularfeld 9, setzen wir ganz wie Néel5) an in der Form: 


= 2 (eM, + M,) | 
und entsprechend : D, = (eM, + BM,). (1) 


n, €, x, ß sind hier Parameter und müßten durch eine tiefergehende Quanten- 
theorie bestimmt werden; ¢ soll beschränkt werden auf ¢ = + 1. 
M, und M, sind dann in bekannter Weise zu bestimmen aus dem System: 


Ma(T) _ B . = Bs, : (2) 
bo 


wo Bg, die Brillouin-Funktion des Index S;, $,, $, die durch (1) definierten 
Molekularfelder und § ein äußeres Magnetfeld sind. 

In bekannter Weise zeigt sich, daß nach (2) unterhalb einer bestimmten Tem- 
peratur 7’,, der Curie-Temperatur, auch für § = 0 eine spontane Magnetisierung 
M,, Mt, + 0 der beiden Teilgitter existiert; in elementarer Weise lassen sich die 
Gln. (2) allerdings nur oberhalb 7’, lösen, wo es bei kleinen Feldern § genügt, die 
Brillouin-Funktion durch das erste Glied der Reihenentwicklung 


= (3) 


anzunähern. In dieser Näherung erhält man die paramagnetische Suszeptibilität 
x(T) in der Form: 


mit: 
me 2E +B) + A + Ap 
+ « B); 
=— 
+ + 2A pe); 
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Wir haben als Abkürzungen eingeführt: 
N,, = Zahl der a-Ionen bzw. -Atome im Einheitsvolumen ; 
N, = Zahl der b-Ionen bzw. -Atome im Einheitsvolumen ; 
N = Zahl aller 'magnetischen Ionen bzw. Atome im Einheitsvolumen. 


4 N, 4 N 
A= Sa(Sa+ y So + 1); 


Man sieht, daß in dieser Bezeichnung cA = C,; c-u = C, die Curie-Kon- 
stanten der Ionen a bzw. 6 sind. 

Man erkennt an (4), daß o + 0 fiir die Abweichungen vom Curie-Weißschen 
Gesetz verantwortlich ist. Durch Einsetzen der für x-Fe bzw. Fe,N gültigen Werte 
von A, u und ¢ = — | sieht man, daß o > 0, solange «, ß > 0 bleiben. D.h. falls 
die Zenersche Struktur in der Natur verwirklicht ist, müßte man dies bei Fe,N an 
x(T) erkennen, während man bei x-Fe wegen des unberücksichtigten Einflusses 
der Ladungsfluktuationen diesen Schluß nicht mit Sicherheit ziehen kann. Mes- 
sungen von 7(7') liegen bei Fe,N meines Wissens leider nicht vor; am Forschungs- 
institut für magnetische Werkstoffe in Jena wird an diesem Problem gearbeitet. 
Es muß allerdings betont werden, daß solche Experimente recht schwierig sind 
vor allem aus 2 Gründen: Fe,N in völlig reiner Form herzustellen, ist nicht ein- 
fach; der geringste Fe-Überschuß aber würde die Messungen stark verfälschen, 
wenn nicht unmöglich machen. Ferner ist Fe,N nur unterhalb 670° C beständig, 
sein Curie-Punkt liegt bei 480° C, also kann man x(7') nur in einem relativ schma- 
len Temperatur-Bereiche messen. Einfacher zu erhalten sind Informationen über 
M(T);die vorhandenen Daten®) zeigen eine deutliche Abweichung vom Weißschen 
Gesetze in ähnlicher Richtung wie bei anderen Ferrimagnetika. Beim x-Fe zeigen 
M(T) und 7(7) keine merkliche Abweichung vom Weißschen bzw. Curie- 
Weißschen Gesetz. Ob dies mit Zeners Struktur verträglich ist oder nicht, 
muß noch dahingestellt bleiben. — (4) gestattet noch die Bestimmung zweier 
Punkte, welche von einigem Interesse sind: 1/y wird 0 am-Curie-Punkt "= T,, 
also: 


Ferner ist die fiir sehr große 7' verwirklichte Curie-Weißsche Gerade: 7 ats: 


Sie spielt die Rolle einer Asymptoten an 1/y(T); ihre Verlängerung schneidet die 
T-Achse im Punkte ©,, dem asymptotischen Curie-Punkte, welcher bestimmt ist 
durch: 


9,= + mß+2ine). (6) 


Leider ist, wie schon erwähnt, das Verhalten von M(7') für 7’ < 7, mit elemen- 
taren Mitteln aus (2) nicht zu gewinnen. Man kann aber leicht zeigen, was für uns 
gegenwärtig am wichtigsten ist, daß (2) in unseren Fällen nicht zufällig Lösungen 
haben kann, welche mit dem alten Weißschen Gesetz übereinstimmen. Wäre 
nämlich das Gegenteil der Fall, so müßten M, und Mt, die gleiche Temperatur- 
abhängigkeit zeigen: M (T) = ¢ -M,(T); es folgten dann aus (2) die Gleichungen: 


[My | = Be Mel (e+ = Bu +B)). 


SC. Guillaud, C. R. Acad. Sci. Paris 223, 1110 (1946). 


| 
we 
We 
sc 
A 
m 
di 
wi 
ha 
Fe 
pa 
a, 
di 
be 
F 
YA 
sä 
E 
K 
WwW 
n 
N 
A 
H 
st 
k 
d 
x 
k 
W 
d 
st 
b 
it 


G. Heber: Zur Struktur des ferromagnetischen Zustandes in x-Fe 159 
welche wegen S, + S, nicht gleichzeitig erfüllt sein können. — Die M(T')-Gesetze, 
welche man statt des Weißschen zu erwarten hat, sind sehr verschieden für ver- 
schiedene &, ß, A, u; wir wollen auf eine allgemeine Diskussion dieser Gesetze im 
Augenblick verzichten. Es soll ja hier nur darauf hingewiesen werden, daß es 
möglich sein muß, aus y und M weiteres experimentelles Material für oder gegen 
die Zenerschen Thesen zu gewinnen. 

Die Berechtigung einer von Zener stillschweigend gemachten Annahme können 
wir aber vom Standpunkt der Molekularfeld-Theorie noch leicht nachprüfen. Es 
handelt sich um die Annahme, daß im energetisch Zen Zustand von «-Fe und 
Fe,N die beiden Teilgitter bei völliger Sättigung (M.| = Mao; |M,| = Myo) anti- 
parallel magnetisiert seien. Wir wissen ja aus Néels Arbeiten), daß bei gewissen 
«, B, € auch teilweise Sättigung energetisch am günstigsten sein kann. Zur Prüfung 
dieser Frage hat man die Energie des Molekularfeldes 


W=— 9 + M, H,) (7) 
bezüglich |M,| und |M,| zu variieren und das tiefste Extremum aufzusuchen. 
Führt man dies aus, so sieht man, daß für «, 8B >0; e=—1 (welche Werte für 


Zeners Strukturen!) von ~«— Fe und Fe,N notwendig zu sein scheinen) tat- 
sächlich energetisch am günstigsten der Zustand |M,| = Moo, |M,| = Moo ist. 
Es sei jedoch bemerkt, daß diese Aussage der Molekularfeld-Theorie unter Um- 
ständen in einer Quantentheorie des Problemes nicht völlig erhalten bleiben wird; 
Überlegungen von Kramers und Mitarbeiter”) zeigen nämlich, daß die lineare 
Kette im Falle der isotropen Heisenberg-Austauschkopplung (im Gegensatz 
zur anisotropen Ising-Kopplung etwa) keine Sättigung der Teilgitter aufweist, 
wenn man die Quantentheorie anwendet. — Es wäre wünschenswert, diese Rech- 
nungen auf 3-dimensionale Gitter auszudehnen*). 

$3. Schließlich sollen noch einige Bemerkungen gemacht werden, welche die 
Natur der Austauschkräfte zwischen benachbarten magnetischen Ionen bzw. 
Atomen betreffen. Wie schon erwähnt, nimmt Zener‘) an, daß die direkte, seit 
Heisenberg ausschließlich untersuchte Wechselwirkung zwischen 3d-Elektronen 
stets negativ, antiferromagnetisch sei und die positive, ferromagnetische Spin- 
kopplung indirekt durch die Leitungselektronen bewirkt wird. Es ist klar, daß 
diese Annahme unabhängig ist von derjenigen über die Disproportionierung im 
x-Fe und seinen Legierungen. Es scheint mir von Bedeutung, dies zu betonen, da 
kürzlich von Kaplan’) Rechnungen über 3d-Austauschintegrale referiert wurden, 
welche positives Austauschintegral und damit ferromagnetische direkte Kopplung 
der Atomspins im Ferromagnetikum ergaben. Zudem hat Verfasser?) kürzlich 
auf der Grundlage der Spinwellentheorie nachgewiesen, daß die indirekte Zener- 
sche Kopplung zwar prinzipiell möglich ist, aber bei Ni und wahrscheinlich auch 
bei Fe und Co keine nennenswerte Stärke erreichen dürfte.. 

Tatsächlich kann man ja die zur Stabilisierung der Zenerschen Spinstrukturen 
in x-Fe und Fe,N nötigen Spinkopplungen ohne Schwierigkeiten aus den von 


*) H.A. Kramers, Physica 18, 101 (1952); P. W. Kasteleijn, Physica 18, 104 (1952). 

8) H. Kaplan, Physic. Rev. 85, 1038 (1952). 

®) G. Heber, Wiss. Z. der Univ. Jena, Jahrg. 1, Heft 1, S. 79; Heft 3, S. 95 (1952); 
Ann. Physik im Druck. 

*) Anmerkung bei der Korrektur: Inzwischen sind Verf. Arbeiten zur Kenntnis 
gekommen, welche die gewünschte Ausdehnung auf 3 Dimensionen erhalten; vgl.: 
P. W. Anderson, Physic. Rev. 86, 694 (1952). 
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Zener meines Erachtens zu Unrecht kritisierten!) N&elschen Betrachtungen") 
verstehen: Im «-Fe ist die direkte Wechselwirkung zwischen nächsten Nachbarn 
negativ, zwischen übernächsten Nachbarn aber positiv; entsprechendes müßte 
dann für Fe,N gelten. 

Eine ähnliche Bemerkung scheint mir notwendig zu sein bezüglich der Uber- 
legungen von Carr?) über die Sättigungsmagnetisierung gewisser binärer Le- 
gierungen. Auch dort werden Annahmen nur über die Koeffizienten des Molekular- 
feldes gemacht (dort mit W44, Wa»; W 42 bezeichnet); ob diese dem direkten oder 
indirekten Austausch entspringen, ist damit nicht im geringsten festgelegt! 

Anders liegen die Verhältnisse im Fe,Si und Fe,Al. Dort sind die magnetischen 
Ionen durch unmagnetische (Si bzw. Al) Ionen getrennt. Es käme also, wenn man 
auf den Zener-Austausch verzichten wollte, nur Super-Austausch in Frage, 

Man kann sich leicht überlegen, daß nach den von Anderson aufgestellten) 
und an einer ganzen Reihe von Experimenten bestätigten Regeln über das Vor- 
zeichen des Superaustausches in diesen Substanzen ein negativer Superaustausch 
zu erwarten wäre. Die beteiligten Ionen haben nämlich hier genau die gleichen 
Elektronenkonfigurationen wie im FeO z. B.; es müßte also wie dort Antiferro- 
magnetismus resultieren. Hier scheint also tatsächlich die Zenersche indirekte 
Wechselwirkung maßgebend am Zustandekommen des Ferromagnetismus be- 
teiligt zu sein. 


10) Verf. hält die erwähnte Kritik für nicht berechtigt, weil ja z. B. auch Zener selbst 
verschiedentlich fordern muß, daß die direkten Austauschintegrale genügend schwach 
negativ seien, damit der indirekte Austauscheffekt sie überwiegen kann. Derartige An- 
nahmen sind unumgänglich, solange keine quantitativen Rechnungen über Austausch- 
integrale verfügbar sind. 

11) L, Néel, Ann. de phys. [11] 8, 237 (1937). 

12) W. J. Carr jr., Physic. Rev. 85, 590 (1952). 

18) P. Anderson, Physic. Rev. 79, 350 (1950); vgl. auch: J. H. van Vieck, J. 
Physique Radium 12, 262 (1951). 


Jena, Theoretisch-Physikalisches Institut der Universitat. 


(Bei der Redaktion eingegangen am 21. Juli 1952.) 
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Demodulation von Ultraschallwellen in Flüssigkeiten 
Von Kurt Altenburg und Siegfried Kästner‘) 
(Mit 4 Abbildungen) 


Inhaltsübersicht 


Es wird die Demodulation einer Ultraschallwelle in einer Flüssigkeit unter- 
sucht. Die Theorie zeigt, daß hierbei der Schallstrahlungsdruck eine wesentliche 
Rolle spielt. Es zeigt sich: 

1..Die NF-Schallamplitude wächst bei konstanter Energiedichte des Trägers 
linear mit dem Modulationsgrad. 

2. Die NF-Schallamplitude wächst bei konstantem Modulationsgrad linear 
mit der Energiedichte des Trägers. 

3. Die 1. Oberwelle der Niederfrequenz wächst bei konstanter Energiedichte 
des Trägers quadratisch mit dem Modulationsgrad. 

Diese Ergebnisse wurden experimentell quantitativ bestätigt. Das Ver- 
fahren ist geeignet die Energiedichte in Schallfeldern zu bestimmen, wobei sich 
gegenüber den üblichen Methoden, die mit dem Radiometer arbeiten, gewisse 
Vorteile ergeben. 


Einleitung 


Breitet sich eine modulierte Ultraschallwelle in einer Flüssigkeit aus, so tritt 
eine Demodulation ein, und man kann die Niederfrequenz mit einem Mikrophon, 
das für Ultraschall unempfindlich ist und nur auf Niederfrequenz anspricht, 
empfangen. Rein subjektive Untersuchungen haben R. W. Boyle und Mit- 
arbeiter?) durchgeführt. Nach Abschluß unserer Arbeit wurden Messungen von 
R. Esche?) bekannt, der diesen Effekt zur Bestimmung der Ultraschallabsorption 
in tierischen Geweben und Kunststoffen — in einer unserem Aufbau ähnlichen 
Apparatur — benutzt, ohne jedoch auf den Mechanismus der Demodulation ein- 
zugehen. Im folgenden soll diese Erscheinung näher untersucht werden. 


Theoretischer Teil 


Die Demodulation von Ultraschallwellen in einer Flüssigkeit hat zwei Ur- 
sachen. Einmal zeigen Analogiebetrachtungen bei entsprechenden elektrischen 
Problemen, daß durch die nichtlinearen zwischenmolekularen Kräfte eine Demo- 


1) Experimentelle Teilergebnisse der Diplomarbeit von S. Kästner, Institut für 
allgemeine Elektrotechnik der T. H. Dresden (1951). 

2) R. W. Boyle, D.K. Framan and G. 8. Field, Canad. J. Res. 6, 102 (1932); 
zit. L. Bergmann, Der Ultraschall, 5. Aufl. (Stuttgart 1949) S. 137. 

3) R. Esche, Akustische Beihefte 1952, 71, Heft 2. 
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dulation eintritt*). Eine Abschätzung zeigt jedoch, daß dieser Effekt von unter- 
geordneter Bedeutung ist!). Er wird vollkommen verdeckt durch die Demo- 
dulation, die durch den Schallstrahlungsdruck bedingt ist. Wir betrachten zu- 
nächst den Fall, daß sich eine unmodulierte Ultraschallwelle ausbreitet. Bekannt- 
lich wirkt auf ein Hindernis im Schallstrahl, dessen Abmessungen groß gegenüber 
der Wellenlänge sind und das den Ultraschall reflektiert, neben dem Ultraschall- 
wechseldruck ein Gleichdruck. Dieser „Langevinsche Strahlungsdruck“ ist 
quantitativ gleich der mittleren Schallenergiedichte e vor dem Hindernis, unab- 
hängig von der Zustandsgleichung des Mediums’). Die Energiedichte e ist pro- 
portional dem Quadrat der Amplitude der Teilchenauslenkung 


e=}ow?a? (1) 
o = Dichte, w = Kreisfrequenz, a = Ampl. d. Teilchenauslenkung. Wird die 
Ultraschallwelle niederfrequent moduliert, dann schwankt die Energiedichte im 


‚Rhythmus der Niederfrequenz. Damit erhalten wir nach (1) an einer Membran 
im Schallfeld folgenden Strahlungsdruck*) P 


P =e (1+ mcos w, t)? (2) 
m = Modulationsgrad € = Energiedichte fiir m = 0 


== Energiedichte des Trägers 
Nach dem Quadrieren erhält man 


2 m cos w,t + cos 20, t) 
(3) 
Am Mikrophon wirkt neben einem erhöhten Strahlungsgleichdruck (P,) 
ein NF-Druck der Grundschwingung (P,), sowie ein NF-Druck der ersten Ober- 
schwingung (P3). 
Experimenteller Teil 


Die Meßapparatur ist in Abb. 1 dargestellt. Ein zweistufiger Hochfrequenz- 
generator erzeugt die niederfrequent modulierte Hochfrequenzwechselspannung, 
die dem Ultraschallgeber, einem Piezoquarz, zugeführt wird. Die Energiedichte 

wurde indirekt ermittelt 


= durch Messung der Hoch- 
frequenzwechselspannung 

Flüssigkeit mit einem Röhrenvoltmeter. 

Steuer. Modular: Uitraschallgeber | Aöten. Der Modulationsgrad wurde 

Sender Stufe vom mit einem speziellen Röhren- 

voltmeter für Spitzenspan- 

Abb, 1. Blockschaltbild nungen nach K.Altenburg’) 


4) K. Altenburg, Frequenz 5, 285 (1951). 

5) G. Richter, Z. Physik 115. 97 (1940). r 

6) Der Schallstrahlungsdruck wird im wesentlichen durch eine Änderung der Dichte 
im Schallfeld verursacht. Dazu ist ein Zu- oder Abfließen von Flüssigkeit erforderlich. 
In Richtung quer zum Schallstrahl findet dieser Ausgleich mit Schallgeschwindigkeit c 
statt. Hat der Schallstrahl einen Durchmesser d, so ist dazu eine Zeit t = d/c erforderlich. 
Der Schallstrahlungsdruck kann dem Rhythmus der Niederfrequenz folgen, solange die 
Periodendauer .7' groß gegenüber der Zeit ¢ ist. Diese Bedingung ist unter den allge- 
meinen Versuchsverhältnissen erfüllt. 

7) K. Altenburg, Frequenz 4, 315 (1950). 
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gemessen. Hierbei läßt sich aus den beiden Spannungswerten mit und ohne 
Modulation der Modulationsgrad berechnen. Bei geringem Modulationsgrad wurde 
zur Erhöhung der Meßgenauigkeit eine Kompensationsschaltung angewendet. 
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Abb. 2. Abhängigkeit der NF-Schallamplitude vom Modulationsgrad 


In 8 cm Abstand vom Sehallgeber befindet sich der niederfrequente Schall- 
empfänger, ein Seignettesalz-Biegeschwinger. Er erzeugt eine geringe nieder- 
frequente Wechselspannung, die 
nach Verstärkung über einen 
Resonanzverstärker®) miteinem 97 
Kopfhörer abgehört bzw. mit uf 
einem Röhrenvoltmeter gemes- 
sen werden kann. Sämtliche 
Messungen wurden mit de- 35} 
stilliertem Wasser und, um den 
Temperaturgang des Mikro- m 
phoniibertragungsmaBes aus- 
zuschalten, bei einer Tempe- Mm 7 
ratur 9 —= 18° durchgeführt. 

Um die Bildung stehender 


m=const (02) 


Wellen zu verhindern, wurde ” 

das Mikrophon gegen den 

Schallstrahl geneigt und die BR 
reflektierte Welle durch Wb 
Schwammgummi absorbiert. Es 

wurden folgende Messungen Abb. 3. Abhängigkeit der NF-Schallamplitude von 
durchgeführt: der Energiedichte des Trägers 


*) K. Feher u. G. Kurtze, Frequenz 4, 72 (1950). 
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1. NF-Schallamplitude für konstante Energiedichte des Trägers e, als Funk- 
tion vom Modulationsgrad m (Abb. 2). 


2. NF-Schallamplitude für konstanten Modulationsgrad m als Funktion der 
Energiedichte des Trägers e, (Abb. 3). 


3. 1. Oberschwingung für konstante Energiedichte des Trägers e, als Funktion 
vom Modulationsgrad m (Abb. 4). 


In Übereinstimmung mit der Theorie zeigen sich folgende Ergebnisse: 


1. Die NF-Schallamplitude wächst bei konstanter Energiedichte des Trägers &, 
linear mit dem Modulationsgrad m. 


2. Die NF-Schallamplitude wächst bei konstantem Modulationsgrad 72 linear 
mit der Energiedichte des Trägers e,. 


3. Die 1. Oberwelle der Niederfrequenz wächst bei konstanter Energiedichte 


des Trägers e, quadratisch mit dem Modulationsgrad m. 


4. Das Verhältnis Oberschwingung zu Grundschwingung wird gleich =. 


018 - 
06 - Umv Const. (4,8646) 
0% + 
- 


002 + 
Modulationsgrad 


02 04 06 08 
Abb. 4. Abhängigkeit der NF-Schallamplitude der 1. Oberschwingung vom Modulationsgrad 


Das Verfahren stellt eine neue Methode dar, um die Energiedichte einer Ultra- 
schallwelle zu messen. Es hat gegenüber dem üblichen Verfahren, bei dem mit 
dem Radiometer der Strahlungsgleichdruck gemessen wird, zwei Vorteile: 


1. Die Messungen des Strahlungsgleichdrucks werden durch auftretende Gleich- 
strömungen in der Flüssigkeit gefälscht. Diese Gleichströmungen können bei 
diesem Verfahren keine Fehler verursachen, da der Strahlungswechseldruck ge- 
messen wird. 


2. Das Meßverfahren ist empfindlicher. L. Bergmann?) gibt für ein be- 
sonders empfindliches Radiometer an, daß noch Energiedichten von 0,8 vb ge- 
messen werden können. Mit der hier beschriebenen Apparatur, die nicht für hohe 


%) L. Bergmann, Der Ultraschall, 5, Aufl. (Stuttgart 1949), S. 133, 
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Empfindlichkeit angelegt ist, war es möglich, noch Energiedichten von 0,25 ub 
zu messen. 

Diesen Vorteilen steht der Nachteil gegenüber, daß der HF-Generator modu- 
liert werden muß. Falls jedoch der Anodengleichspannung eine gewisse Brumm- 
spannung überlagert ist, kann die Brummfrequenz mit dem Mikrophon abgehört 
werden, so daß sich eine besondere Modulation erübrigt. Zum Schluß sei noch 
darauf hingewiesen, daß der hier geschilderte Demodulationseffekt auf Grund 
seines geringen Wirkungsgrades für eine technische Anwendung kaum in Frage 
kommt. 


Z.Z. Berlin-Adlershof, Deutsche Akademie der Wissenschaften zu Berlin, 
Laboratorium für Kunststoffe. 


(Bei der Redaktion eingegangen am 2. Juli 1952.) 
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Rauschmessungen an dünnen Bleisulfid- Phot schichten 
Von F. Eckart 
(Mit 3 Abbildungen) 


Inhaltsiibersicht 


Es wird iiber Rauschmessungen an chemisch hergestellten und einer Sauer- 
stoff-Behandlung unterworfenen Bleisulfid-Schichten berichtet. Es zeigt sich, daß 
der Absolutwert der Rauscheinströmung proportional dem Zellenstrom und unab- 
hängig davon ist, ob die Elektronen thermisch oder lichtelektrisch abgelöst werden. 
Der Dunkelstrom ist ebenso wie die mittlere Rauscheinströmung der Feldstärke 
proportional. Die mittleren Lebensdauern der Ladungsträger werden unter Zu- 
grundelegung einer aus Halleffektsmessungen ermittelten mittleren Beweglichkeit 
abgeschätzt. 


I. Einleitung 


Die elektrischen und lichtelektrischen Eigenschaften von Festkörpern können 
durch ,,Stérstellen‘‘ wesentlich beeinflußt werden, wobei unter Störstellen ganz 
allgemein Fehlordnungszustände verstanden werden sollen, an denen relativ 
nur wenig Atome, Ionen oder Elektronen des Gitters sowie zusätzliche Elektronen 
oder Elektronenlücken beteiligt sind. Der Einfluß chemischer, d.h. gittereigener 
oder gitterfremder Störstellen ist zweifellos am ausgeprägtesten und auch bereits 
weitgehend untersucht. Auch beim Bleisulfid sind die elektrischen und lichtelek- 
trischen Eigenschaften wesentlich von Art und Konzentration der Störstellen ab- 
hängig. Ultrarotempfindliche Photoschichten sind für die Anwendung dann be- 
deutungsvoll, wenn es gelingt, Schichten herzustellen, deren Nachweisempfindlich- 
keit hinreichend groß und deren Konstanz für Meß- und Registrierzwecke aus- 
reichend ist. Dabei wird als Nachweisempfindlichkeit jene Bestrahlungsstärke 
definiert, bei der die Zahl der lichtelektrisch abgelösten Elektronen oder Defekt- 
elektronen gleich der mittleren Schwankung der Trägerdichte ist. Die Schwan- 
kungserscheinungen werden allgemein als Rauschen bezeichnet. Die vorliegenden 
Untersuchungen sollten einen Beitrag zur Kenntnis der Rauscheigenschaften von 
Bleisulfid-Photoschichten liefern. 


II. Meßanordnung und Ergebnisse 


Die Messung der Rauschströme erfolgte mit einem 3stufigen Batteriever- 
stärker mit einer Halbwertsbreite von etwa 6,8 kHz und einer mittleren Verstär- 
kung von 1,500fach. Der Eigenstörpegel des Verstärkers ist durch das Nyquist- 
Rauschen des ersten Gitterableitwiderstandes bestimmt und liegt für 10 M2 
bei 2-10-12 [A]. Als Rauschstrom wird stets die Wurzel aus dem mittleren Rausch- 
stromquadrat verstanden, Eine optische Einrichtung, in die Zelle und Vorver- 
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stärker unmittelber eingebaut sind, gestattet auch Rauschmessungen bei Belich- 
tung durchzuführen. Die Ausgangsspannungen werden mit einem in Effektiv- 
werten geeichten Röhrenvoltmeter gemessen. Für verschiedene lichtelektrisch 
empfindliche mit Sauerstoff behandelte Bleisulfid-Schichten sind die Rauschein- 
strömungen ermittelt und in Abb. 1 im doppelt logarithmischen Maßstab in Ab- 
hängigkeit vom Gleichstrom und in Abb. 2 mit und ohne Belichtung aufgetragen. 
Der mittlere Rauschstrom ist demnach proportional dem Zellenstrom, wobei es 


1 5 
N lA] 
fal Dunkelstrom # 
Chemisch verspiegelte 
- 
mit Sauersto 0 
; sensible 
Schicht Net 
oben 
jedoch mitwentg § im Dunkeln 
| Souerstoff 2} bei Belichtung 
10 
Elektronenstrom —> 


Abb. 1. Rauschstrom in Abhängigkeit vom Abb. 2. Rauschstréme im PbS mit 
Dunkelstrom fiir PbS-Schichten und ohne Belichtung 


offenbar gleichgiiltig ist, ob der Zellenstrom von thermischen Elektronen oder 
von lichtelektrisch abgelösten Leitungsträgern herrührt. Diese Verhalten zeigen 
allerdings nur sauerstoff-sensibilierte Schichten. Nichtsensibiliserte Schichten 
verhalten sich offenbar anders. Hierüber wird in Kürze noch berichtet. 


Ill. Diskussion der Ergebnisse 


Die mit Rauschen bezeichneten Schwankungserscheinungen kommen durch 
zeitliche statistische Dichteschwankungen von Ladungsträgern zustande. Gi- 
solf!) hat für Halbleiter dieses mittlere Schwankungsquadrat durch Fourier- 
zerlegung der durch Ablösung und Wiederanlagerung der Ladungsträger hervor- 
gerufenen Rechteckimpulse und durch Summierung der dadurch bedingten 
Rauschströme das effektive Rauschstromquadrat ermittelt, und zwar gilt für 
niedrige Frequenzen (ft < 1): 


= eb “Af = = t Af = - Af {Al; 
wobei e = Elektronenladung [A sec] = mittlere Lebensdauer in [sec] 
b = Beweglichkeit [Y% = | f = Frequenz [sec-!] 
& = Feldstärke [Y/em] V = Volumen der Schicht [em*] 
Af = Bandbreite [sec~*] n = Elektronendichte in [em”?] 
1 = Länge des HL [em] I = Gleichstrom in [A] 


q = Querschnitt [cm?]. 


1) Gisolf, Physica 15, 826 (1949). 
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Die obige Beziehung gilt zwar nur für homogene Schichten, bei denen die Felda 
stärke nicht wesentlich vom Ort abhängt, doch sind geringe Unterschiede vom 
geringem Einfluß, da in obiger Beziehung lediglich das Verhältnis von E/yR 
eingeht. Wie Abb. 3 zeigt, ist die 
Rauscheinströmung der wirksame 
Feldstärke gleichfalls proportional. Aug 
den Rauschmessungen folgt unmittel 
bar das Verhältnis der Lebensdaue 
zur Dichte der Ladungsträger (Tabelle lj 
Im allgemeinen ist r eine Funktiog 
von n. Die Beweglichkeit kann da 
gegen im untersuchten Bereich alg 
nahezu konstant angenommen werden 23 
Mit einem aus Halleffekts-Messungeg 
von Sauerstoff behandelten Bleisulfid 
Schichten ermittelten mittleren Wert 


der Beweglichkeit von 10 Vjem ergeben 


sich bei Gültigkeit des Ohmschen Gee 
setzes als mittlere Lebensdauer 
10  Ladungsträger Werte zwischen 1,3-104 


Abb. 3. Rausch winged Abh k und 3,7- 10-5 [sec], die mit Literatur 
von der Feldstärke für PbS-Schichten angaben übereinstimmen ®)®). 


Schicht Zellen- Rauschstrom 


| Lebensdauer r Träger- ‘ mittlere 
Nr. strom[A], - [A] I 
| | 


Trägerdichte n dichte n | Lebensdauer #3 
[em? sec] [em?] [sec] 


7,4 . 10-12 1,5 - 10-20 2.5 - 1015 
.8- 10712 2.1018 
10% 13.102 5.10% 7.1028 
10% 12 . 10-12 10 . 10-20 14 . 1015 


2) H. Diedrich, Diplomarbeit Berlin 1950. 
8) H. Pick, Z. Physik 126, 12 (1949). 
4) H. Miser, Z. Physik 129, 504 (1951). 


Berlin-Buch, Institut für Festkörperforschung der Deutschen Akademie det 
Wissenschaften. 


(Bei der Redaktion eingegangen am 1. Juli 1952.) 
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Leuchtstoff-Schichten und Folien mit Polyvinyl-Alkohol 
als Schutzkolloid und Bindemittel 


Von F. Eckart 
(Mit 4 Abbildungen) 


Inhaltsiibersicht 


Versuche zeigen, daß es möglich ist, an Stelle des bisher als Bindemittel‘ ver- 
wendeten Wasserglases für im Vakuum ausheizbare hochauflösende sedimen- 
tierte Leuchtschirme wasserlöslichen und nachträglich härtbaren Polyvinyl- 
Alkohol zu verwenden. Der Polyvinyl-Alkohol wirkt dabei gleichzeitig als Schutz- 
kolloid. PVA-Leuchtschirme sind gegen mechanische Beschädigungen wesentlich 
unempfindlicher und können bei höherem PVA-Gehalt, analog dem PVA i ea a 
stoff-Folien, abwaschbar gemacht werden. 


I. Einleitung 


Zur unmittelbaren physiologischen Beurteilung von Elektronenbildern dienen 
Leuchtschirme und Leuchtfolien. Als Erkennbarkeitsschwelle bezeichnet man 
jene Anregungsenergie, die notwendig ist, um ein gegebenes Bild nicht nur als 
Lichteindruck wahrzunehmen (Wahrnehmbarkeitsschwelle), sondern bei der das 
Testbild auch als solches erkannt werden kann. Eine niedrige Erkennbarkeits- 
schwelle setzt daher hohes Auflösungsvermögen und hohe Lichtausbeute (Leucht- 
dichte je eingestrahlte Energieeinheit) voraus. Diese Forderungen sind gleich- 
zeitig optimal nicht erfüllbar. Es wurden daher je nach Verwendungszweck ver- 
schiedene Herstellungsverfahren entwickelt, die man folgendermaßen kenn- 
zeichnet: 


1. Perl-Verfahren 


Der Leuchtstoff wird auf eine Bindemittel-Folie aufgebracht und ergibt wegen 
der fehlenden Absorption der erregenden Strahlung im Bindemittel Leucht- 
schirme höchster Leuchtdichte aber geringerer Auslösung, da der Leuchtstoff 
nicht kolloidal sedimentiert wird. 


2. Sedimentations-Verfahren 


Der Leuchtstoff wird zusammen mit Bindemittel und Schutzkolloid kolloidal 
suspendiert und auf die Unterlage sedimentiert. Da der Leuchtstoff im Binde- 
mittel eingebettet ist, so ist die Leuchtdichte geringer als bei Leuchtschirmen 
nach dem Perl-Verfahren. Die Auflösung ist wegen der kolloidalen Suspension 
des Leuchtstoffes größer als bei Leuchtschirmen nach dem Perl-Verfahren. 


Ann, Physik. 6, Folge, Bd. 11 11b 
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3. Aufdampf-Verfahren 
Der Leuchtstoff wird entweder aus der Dampfphase niedergeschlagen oder 
dünnste Aufdampf-Schichten nachträglich leuchtfähig gemacht. Bei feinkristalli- 
nem Niederschlag der Substanz erhält man Schirme höchster Auflösung. Dieses 
Herstellungs-Verfahren ist bisher auf ZnO-Schirme beschränkt, deren Leucht- 
dichte die der sedimentierten Schirme jedoch nicht erreichen. 


4. Einkristall-Leuchtschirme 


Bisher ist lediglich in einer Veröffentlichung!) auf das hohe Auflösungsver- 
mögen derartiger Kristalle hingewiesen worden. 


5. Uranglas-Leuchtschirme 
Uranglas der Firma Schott & Gen., Nr. GG 17, ist bei Elektronenstrahl- 
Anregung leuchtfähig. Die Auflösung liegt gleichfalls in der Größenordnung der 
der Aufdampf-Schichten*), die Leuchtdichte in der Größenordnung der Auf- 
dampf-Leuchtschirme. 


Das Auflösungsvermögen von Leuchtschirmen ist im wesentlichen durch die 
Streuung der einfallenden Elektronenstrahlen ‘und durch die Lichtstreuung in 
der Schicht begrenzt*). In der Technik wird heute für hohe Anforderungen fast 
ausschließlich nach dem Sedimentations-Verfahren gearbeitet. Im folgenden soll 
über ein Herstellungs-Verfahren für Leuchtstoff-Schirme und Folien berichtet 
werden, das gegenüber den bisherigen Verfahren offenbar Vorteile aufweist, die 
sich vor allem bei niedrigen Anregungs-Spannungen auswirken. 


II. Herstellung von Leuchtschirmen und Folien mit Polyvinyl-Alkohol als 
Bindemittel und Schutzkolloid 


Nach dem in der Technik gebräuchlichen Sedimentations-Verfahren wird die 
Leuchtsubstanz in einem Wasserglas-Essigsäure-Gemisch kolloidal suspendiert 
und sedimentiert. Durch Zusatz von Essigsäure entsteht freie Kieselsäure, die 
kolloidal in Lösung bleibt und als Schutzkolloid für die Leuchtstoff-Teilchen 
wirksam ist. Wasserglas dient als Bindemittel. Dieses Verfahren liefert gegenüber 
dem Perl-Verfahren zwar gleichmäßigere Schichten, vor allem auch deshalb, 
weil man größte und kleinste Teilchen vor dem eigentlichen ‚Gießen‘ der Schirme 
heraussedimentiert, aber Leuchtschichten geringerer Ausbeute, da die Substanz- 
teilchen allseitig in das energie-absorbierende Bindemittel eingebettet sind. Da 
die Absorption vom spez. Gewicht des verwendeten Bindemittels abhängig ist, 
so konnte eine Steigerung der Lichtausbeute erwartet werden, wenn man zu ge- 
eigneten organischen Bindemitteln übergeht. 

Versuche haben nun gezeigt, daß sich wasserlöslicher Polyvinyl-Alkohol (PVA) 
als Schutzkolloid und Bindemittel für Leuchtschichten und Folien ebenfalls 
gut eignet. 


Der polymere Vinyl-Alkohol (-08,-69-) wird bis zu hohem Polymeri- 


n 


sationsgrad technisch hergestellt. Zur Herstellung des Polyvinyl-Alkohols geht 


1) M. v. Ardenne, Z. techn. Physik 119, 397 (1942). 
2) Kopp u. Möllenstedt, Optik 2, 283-300 (1947), 
*) O. Hachenberg, Telefunken-Bericht (unveröffentlicht). 
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man von organischen Polyvinyl-Estern halogenierter Säuren aus, die mit Ammoniak 
in alkalischer Lösung verseift werden. 

PVA ist geruch- und geschmacklos, seine scheinbare Dichte schwankt von 
0,054 bis 0,360 g/em*. Die Dichte einer 5proz. wäßrigen Lösung beträgt beispiels- 
weise 1,0079 g/cm*. Bereits bei einem Festkörpergehalt von 30% geliert die 
Lösung. Die Dichte fester PVA-Folien beträgt etwa 1,2 g/cm*. Die Dichte von 
festem Natrium-Silikat beträgt demgegenüber etwa 2,4 g/cm?, ist also um etwa 
den Faktor 2 größer als die einer PVA-Folie. 

PVA ist in den meisten anorganischen und organischen Lösungsmitteln un- 
löslich. Beim Verdunsten von wäßrigen PVA-Lösungen erhält man äußerst zähe 
Filme. PVA verhält sich als Kolloid und erleidet beim Eintrocknen keine irre- 
versiblen Zustandsänderungen. Im Gegensatz zu Gelatine quillt PVA in Wasser 
nur sehr wenig und löst sich bereits in der Kälte auf. 

Wasserunlöslicher PVA wird erhalten, wenn man das Harz auf höhere Tem- 
peraturen erhitzt oder wenn man die Filme mit einer wäßrigen Cupro-Ammonium- 
Hydroxyd-Lösung behandelt. 

Zur Herstellung von PVA-Leuchtschirmen wurde die Leuchtsubstanz in einer 
etwa 3proz. PVA-Lösung suspendiert und auf Glasplatten sedimentiert. Zum 
Vergleich wurden jeweils analog hergestellte Wasserglas-Essigsiure-Leucht- 
schirme verwendet. 


III. Eigenschaften von PVA-Leuchtschirmen bei Elektronen-Anregung 


PVA-Leuchtschirme und -Leuchtfolien sind im Vakuum bis 400° C beständig. 
Beim Ausheizen an Luft tritt bereits bei 300° C eine schwache Gelbfärbung auf, 
die indessen die Eigenschaften nicht wesentlich ‚beeinflußt. Immerhin ist es 
zweckmäßig, beim Ausheizen an Luft eine Temperatur von 300° C möglichst 
nicht zu überschreiten. 

Zur Prüfung der Schirme bei Elektronen-Anregung wurde das in Abb. 1 skiz- 
zierte Versuchsrohr benutzt. Der Leuchtschirm konnte bei geeigneter Defo- 
kussierung wenigstens bei Strahlspannungen bis etwa 3 kV gleichmäßig aus- 
.geleuchtet werden. Bei den Messungen, über die wir im folgenden berichten, 
wurde der Elektronenstrahl auf den Leuchtschirm punktförmig abgebildet und 
mit magnetischen Horizontal- und Vertikal-Ablenkspulen rasterförmig über den 
Leuchtschirm bewegt. 

Der Versuchsschirm bestand jeweils aus zwei Hälften, von denen die eine 
Hälfte den unter sonst gleichen Bedingungen (gleiche Belegungsstärke und Leucht- 
substanz) hergestellten Wasserglas sedimentierten Vergleichs-Schirm tragt.(Abb. 1). 

In Abb. 2 sind die relativen Leuchtdichte-Unterschiede (Leuchtdichte des 
Wasserglas-Schirms willkürlich = 1 gesetzt) von PVA- und Wasserglas-Leucht- 
schirmen in Abhängigkeit von der Strahlspannung für zwei verschiedene Leucht- 
stoffe aufgetragen. Bei Strahlspannungen von 500—1000 V beträgt der Leucht- 
dichte-Unterschied von PAV- und Wasserglas-Leuchtschirmen mehr als 100%. 
Wie zu erwarten, nehmen mit steigender Strahlspannung die Leuchtdichte-Unter- 
schiede immer mehr ab. Bei einer Strahlspannung von etwa 7 kV sind die Leucht- 
dichten beider Schirme etwa gleich. Bei noch höheren Strahlspannungen wird der 
Wasserglas-Schirm dem PVA-Schirm überlegen. 

In Abb, 3 ist die Leuchtdichte in Abhängigkeit von der Strahlspannung für 
die beiden Leuchtschirme aufgetragen. Demnach steigt bei Strahlspannungen bis 
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zu etwa 8 kV die Leuchtdichte etwa quadratisch mit der Strahlspannung an (bei 
den in der Abb. 3 wiedergegebenen Kurven ist die Schwellenspannung U, nicht 
berücksichtigt, weshalb die Neigung der Geraden etwas größer als 2 ist). Bei 
noch höheren Spannungen strebt die Kurve immer mehr einem „Sättigungswert“ 
zu, der allerdings für Wasserglas-Leuchtschirme offenbar entsprechend der großen 
Absorption höher liegt als für PVA-Schirme. 


Ablenk - Magn. Konzentrier- 
spule 


Up = 6.3V Leuchtschirm 
Up =-100 v 
U = 1 AV 
Ug = 10kKV (15 kV) 
Ug | In |tsec| Bemerkung 
hne Beleucht. 
b| 70 | 30 | 60 
C| 10 | 40 | 60 
@| 10 | 40 | 60 
e| 10 | 40 | 120 
Fi 10 | 40 | 200 
g| 10 | 40 | 400 
h| 15 | 40 | 400 
f Fokussie 
20 | eo 


Abb. 1. Versuchsrohr zur Untersuchung der Leuchtschirme bei Elektronen-Anregung 
und Photographie eines Wasserglas-PV A-Schirmes 


Vergleichsversuche mit Al-belegten Schirmen (Abb. 4) zeigten, daß die mit- ° 


geteilten Ergebnisse nicht durch Aufladungen verfälscht sind. In der Abb. 4 
sind für die Leuchtsubstanz H,-weißblau (ZnS-ZnSe-Ag)] die Leuchtdichten in Ab- 
hängigkeit von der Strahlspannung für PVA- und Wasserglas-Schirme, und zwar 
jeweils mit und ohne Al-Folie aufgetragen. Die Al-Folie wird mit zunehmender 
Strahlspannung durchlässiger, d. h., die Werte nähern sich mit steigender 
Strahlspannung mehr und mehr den Werten unfolierter Schichten. 
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Die Leuchtdichte B ist in Abhängigkeit von der Stromdichte j und der Strahl- 
spannung U darstellbar durch die Gleichung: 
B=A-j(U-U,); (1) 
wobei A = Konstante 
U, = Schwellenspannung, bei der das Leuchten des Phosphors einsetzt. 
Die Schwellenspannung U, fiir die Leuchtsubstanz H,-weiß-blau mit PVA als 
Bindemittel wurde zu 300 V und für Wasserglas als Bindemittel zu 360 V ermittelt, 
Die Schwellenspannung ist 


demnach, wie zu erwarten - 

it Wossergli 
wesentlich von der Art des 
Bindemittels abhängig. 

R 107 A mit PYA 
x 
{av 
& 
2S/znSe/ag 
(Hywel 
hrsubstanz M 31 
10 
ZnS (M31) Ä 
S 
S 
& 
1201 0 -D-W% 10 kV 
celat. Leuchtdichte -Unterschiede Spannung U— 


Abb. 2. Relative Leuchtdichte-Unter- Abb.3. Leuchtdichte in Abhängigheit von der 

schiede von PVA- und Wasserglas- Strahlspannung bei konstantem Strahlstrom für 

Leuchtschirmen (bzw. auf Wasserglas- gleiche Leuchtsubstanz und verschiedene Binde- 
Vergleichsschirm) mittel 


Nach den obigen Ergebnissen hat der Exponent in Gl. (1) für Spannungen bis 
zu etwa 4 kV den Wert x = 2, d.h. die Leuchtdichte steigt quadratisch mit der 
wirksamen Spannung (U—U,) an. 

PVA-Leuchtschirme erweisen sich somit bei Strahlspannungen unter 7 kV 
den Wasserglasschirmen an Lichtausbeute überlegen. Das Auflösungsvermögen 
bei elektronischer Anregung wurde zwar nicht bestimmt, doch ist zu beobachten, 
daß vor allem bei höheren Strahlspannungen der PVA-Schirm kontrastreicher 
erscheint. Ein weiterer Vorteil bietet sich darin, daß der PVA mit der suspen- 
dierten Leuchtsubstanz zu einem gleichmäßigen und löcherfreien Leuchtstoff- 
film ausgegossen und eventuell gehärtet werden kann. Der PVA-Leuchtstoff- 
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film oder die PVA-Leuchtstoff-Folie ist so elastisch, daß sie weder durch Wärme 
noch im Vakuum reißt. Gehärtete Leuchtstoff-Folien sind ferner abwaschbar, 
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& Leuchtsubstanz 
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10 / 1 Wosserglas 
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I mit 
mossergias 


mıt Wosserglas 
u mit Al-Folie (sehr dick) 
u. ohne Zwischenhlie 


~ 


10 kV 


Spannung U ——~ 


Abb. 4. Abhängigkeit der Leuchtdichte von der 
bei Strahlstrom mit 


Strahlspannung 


konstantem 


und ohne Al-Folie 


Solche hochauflösenden 
Leuchtstoff-Folien können in 
die verschiedenartigsten Ge- 
rate eingebaut und auch aus- 
gewechselt werden. Als An- 
wendungsgebiet seien bei- 
spielsweise genannt: Leucht- 
stoff-Folien für Braunsche 
Röhren, Bildröhren, Bild- 
wandler-Röhren, Elektronen- 
Mikroskope, Röntgen-Verstär- 
ker-Folien und dergleichen ®). 

Die Vorteile der aufge- 
perlten Schirme können nun- 
mehr mit denen der sedimen- 
tierten Schirme vereint wer- 
den. Zu diesem Zweck wird 
zuerst ein dünnster PVA-Film 
gegossen und dann auf diesem 
Film der in Äther, Alkohol 
oder in sonst einer geeigneten 
Lösung suspendierte Leucht- 
stoff aufsedimentiert. Die or- 
ganische Suspensions-Flüssig- 
keit wird verdampft und der 
Leuchtstoff haftet in kollo- 
idaler Verteilung auf der da- 
nach gehärteten und damit 
wasserunlöslich gewordenen 
PVA-Folie. 


4) Bekanntlich werden ge- 
färbte und gestreckte PVA- 
Folien als Polarisationsfilter in 
der Optik verwendet. 


Berlin-Buch, Institut fur Festkörperforschung der Deutschen 
Akademie der Wissenschaften. 


(Bei der Redaktion eingegangen am 5. Juli 1952.) 
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Untersuchung von Geiger-Müller-Zählrohren 
mit Röntgenimpulsen 


Von P. Müller 
(Mit 12 Abbildungen) 


Inhaltsübersicht 


Die Leistungsfähigkeit einer neu entwickelten Röntgenimpulsapparatur wird 
mit einigen Daten charakterisiert. Ihre Verwendbarkeit zur Untersuchung des 
Mechanismus der Zählrohrentladung wird in verschiedenen Abbildungen gezeigt. 
Eine Diskussion der Wirkung der Raumladung führt zu einer Ergänzung der 
Vorstellungen vom Entladungsablauf. 


In den vergangenen Jahren haben Pepinsky u.a.!) eine Röntgenröhre mit 
steuerbarem Wehneltzylinder benutzt, um die Totzeit von Zählrohren zu be- 
stimmen. D. Willard und C. G. Montgomery?) führten Nachentladungs- 
messungen mit einer Röntgenanlage aus, die durch einen Radar-Modulator ge- 
steuert wurde. Für die vorliegende Arbeit wurde eine Röntgenimpulsapparatur 
mit einem Schaltventil entwickelt, die es gestattet, periodische Röntgenimpulse 
im Frequenzbereich von niedrigsten Frequenzen bis zu 10 kHz (und auch Einzel- 
impulse) mit einer einstellbaren Zeitdauer von etwa 10-5 sec bis zu 2,0-10-* sec 
zu erzeugen. (Natürlich kann die Impulsdauer nicht größer als die Periodendauer 
eingestellt werden.) Die Impulsspannung kann bis zu 60 kV betragen, der erreich- 
bare Impulsstrom ist größer als 100 mA. 

Läßt man die Röntgenimpulse auf ein Zählrohr einwirken, so kann man die 
in diesem ausgelösten Entladungen leicht mit einem Kathodenstrahloszillographen 
beobachten und auch photographieren. Die Abb. 1 bis 4 zeigen einige typische Bil- 
der, wie man sie erhält, wenn man einfache Röntgenimpulse benutzt., Die Zählrohr- 
entladungen werden zunächst über einen Breitbandverstärker verstärkt und dann 
einem System eines Zweistrahl-Oszillographen zugeführt. Mit Hilfe des anderen 
Systems wird der Strom durch die Röntgenröhre wiedergegeben. Die benutzte 
Impulsfrequenz betrug 1000 Hz, die Ablenkfrequenz am Oszillographen 500 Hz. 
Das Zählrohr wurde mit einem Außenwiderstand von 800 KQ betrieben, seine 
Eigenkapazität (+ Schaltkapazität usw.) betrug 20 pF. Je höher die Zählrohr- 
spannung über dem Einsatzpunkt liegt, umso mehr Nachentladungen beobachtet 
man. Die kurz nach dem echten Impuls auftretenden Nachentladungen sind zu- 
nächst klein und erreichen nach einer gewissen Zeit die normale Impulsgröße. 
Dies liegt nicht am Absinken der Spannung am Zähldraht, sondern an der Ab- 


si das Pepinsky, P. Jarmotz, H. M. Long and D. Sayre, Rev. Sci. Instrum. 19, 
48). 
*) D. Willard and C. G. Montgomery, Rev. Sci. Instrum. 21, 520 (1950). 
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nahme der Feldstärke im Innern des Zählrohres infolge der entstehenden positiven 
Raumladung. Es läßt sich also auch mit Hilfe von Nachentladungsmessungen 
die Erholungscharakteristik gewinnen (s. Abb. 3). Interessant ist, daß die 
Nachentladungen sich in einem breiten Zeitraum nach echten Entladungen bemerk- 
bar machen. Sie entstehen im wesentlichen durch das Auftreffen der positiven 


Abb. 1. U=11W0V Abb. 2. U = 1120V 


Abb. 3. U = 1160V Abb. 4. U =1200V , 


Abb. 1—4. Mit Röntgenimpulsen von 1000 Hz ausgelöste Entladungen eines Zählrohres 

und dadurch verursachte Nachentladungen bei verschiedenen Spannungen. Oberer 

Strahl: Strom durch die Röntgenröhre: Unterer Strahl: Zählrohrentladungen. Während 

der Belichtung durchliefen beide Strahlen das Bild 50 mal. Ein Impuls wird durch. den 

dunkelgesteuerten Rücklauf nicht sichtbar. Die in den Abb. links erkennbaren Nachent- 
ladungen gehören zu diesem Impuls 


Ionen auf der Kathode. Diese Ionen kommen also nicht gleichzeitig am Zylinder 
an. Die Ionenhülle, welche sich um den Zähldraht in seiner ganzen Länge aus- 
breitet, kann daher kein dünner Schlauch sein. Man vermag diesen Befund zu 
verstehen, wenn man die Feldstärke (in Abhängigkeit vom Radius) berechnet, wie 
sie durch eine positive Raumladung mit endlicher radialer Ausdehnung entsteht. 
Es ergibt sich aus dem Verlauf der Feldstärke kurz nach Bildung der Ionenhiille, 
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tiven daß die äußeren Ionen mit wesentlich größerer Geschwindigkeit zur Kathode 
ngen laufen, als die weiter innen dicht am Zähldraht gebildeten. Die Dicke der Ionen- 
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Abb. 5a—c. Berechnete Feldstärke im 
Zahlrohr. a) ohne Raumladung; b) nach 
Bildung der Ionenhülle und Sammlyng 
der Elektronen am Zähldraht. c) nach 
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beschriebene und wird für eine Näherungsbetrachtung vernachlässigt. In der Abb. 
5a ist der Feldverlauf im Zählrohr ohne Raumladung und in Abb. 5b kurz nach 
Bildung der Ionenhülle dargestellt. In Abb. 5c ist der aus einer Näherungsrechnung 
folgende Feldverlauf nach einer Zeit von t= 1,0 - 10-5sec wiedergegeben. Die La- 
dung pro Längeneinheit ist dabei zu 1,0 - 10-:°Coulomb/em angenommen. Es ist klar, 
daß der von Stever?) definierte „kritische Radius‘ hierdurch eine etwas andere 
Bedeutung bekommt. Gleichzeitig wird aber auch das Scheitern der Berechnung 


Abb. 6. U = 1120V; =1,35-10- s Abb. 7. U = 1120V; 4-2 =0,50-10-*s 


Abb. 6—8. Messung der Totzeit eines Zähl- 

rohres mit Doppelréntgenimpulsen von 

1000 Hz. t,_,: Zeitabstand der Impulse 1 und 

2. Zeitabstand der entsprechenden Doppel- 
impulse: 1,0-10-*s 


Abb. 8. 7 = 1120V; = 0,37 - 10-4 s 


der Totzeit mit Hilfe dieses ,,kritischen Radius“ verständlich. H. den Hartog 
und F. A. Muller‘) haben 1949 auf anderen Messungen beruhende ähnliche Über- 
legungen angestellt und eine genauere Beziehung für die Totzeit abgeleitet. Der 
Feldverlauf wurde von den genannten Verfassern jedoch nicht diskutiert. 

Mit der beschriebenen Röntgenimpulsanlage lassen sich nun auch Doppel- 
impulse mit veränderbarem Abstand erzeugen, die eine unmittelbare Messung der 
Totzeit auf dem Leuchtschirm ermöglichen. Die Abb. 6, 7 und 8 geben Beispiele 
für solche Messungen. Je dichter der zweite Röntgenimpuls dem ersten folgt, um 
so kleiner ist die entsprechende zweite Entladung im Zählrohr. Es liegt das wie- 
derum an der Schwächung der Feldstärke in der Nähe des Zähldrahtes infolge der 
Raumladungshülle. 


3) H. G. Stever, Physic. Rev. 61, 38 (1942). 
‘) H. den Hartog and F. A. Muller, Physica 15, 789 (1949). 
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Um die Abhängigkeit der Totzeit von der Größe der Impulse direkt sichtbar 
zu machen, wurde eine Methode ausgearbeitet, die es gestattet, einen dritten 
Röntgenimpuls mit wählbarem Abstand an den Doppelimpuls anzuhängen. Aus 


Abb. 9. U = 1120V; 4,_, = 1.08- 10-*s Abb. 10. U = 1120V;#,_, =1,07-10-‘s 

1,03;- 10" 4, = 0;37-10-*s 

Abb. 9—10. Messung der (von der Entladungsgröße) abhängigen Totzeit eines Zähl- 

rohres mit Dreifach-Röntgenimpulsen von 1000Hz. Mit t,_s und t,_, sind die Zeitabstände 

der Impulse 1 und 2, bzw. 2 und 3 bezeichnet. Der Abstand der entsprechenden Drei- 
fachimpulse ist 1,0.10-®s 


Abb. 9 ist zu ersehen, daß bei gleichen Impulsabständen die dritte Entladung im 
Zählrohr schon wieder eine wesentlich größere Amplitude hat als die zweite. 
Abb. 10 zeigt einen Zeitabstand des dritten Impulses, der etwa dem Abstand der 


Abb. 11. U = 1080V; t,_, = 0,86 - Abb. 12. U = 1080 V;#,_, = 0,4:10-'s 
= 0,86 + 10-'s te, = 0,4°10-* 8 


Abb. 11—12. Zählrohrbetrieb im Proportionalbereich 


Doppelimpulse in Abb. 8 gleicht. Während die dritte Entladung der Abb. 10 
noch gut aufgelöst wird, ist die mit nahezu gleichem Zeitabstand folgende zweite 
Entladung in Abb. 8 nicht nachweisbar. 


12% 


Abb. 
nach 
mung 
e La- 
klar, 
ndere 
: 
5 
ber- 
der 
um 
vie: 
der 
|| 


180 Annalen der Physik. 6. Folge. Band 11. 1952 


Zur Verdeutlichung der bekannten Erscheinung, daß die Totzeit im Propag 
tionalbereich wesentlich kleiner ist als im Auslösebereich, sind in den Abb. 11 umm 
12 zwei Aufnahmen wiedergegeben, die mit dem gleichen Zählrohr mit ei 
Spannung von 20 V unter dem Einsatzpunkt gewonnen wurden. Man sieht, da 
jetzt die drei Röntgenimpulse vom Zählrohr recht gut aufgelöst werden. 


Die beschriebenen Untersuchungen werden z.Z. noch fortgeführt. Herm 
Prof. Dr. A. Eckardt sei an dieser Stelle für die Anregung zu dieser Arbeit sow 
für wertvolle Förderung besonders gedankt. 


Jena, Technisch-Physikalisches Institut der Universität. 


(Bei der Redaktion eingegangen am 16. Juni 1952.) 
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